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ABSTRACT 

Using UV spectra obtained with FUSE, HST, and/or IUE together with higher resolution 
optical spectra, we determine interstellar column densities of 12 CO, 13 CO, and/or C2 for ten 
Galactic sight lines with E(B - V) ranging from 0.37 to 0.72. The Ar(CO)/JV(H 2 ) ratio varies 
over a factor of 100 in this sample, due primarily to differences in iV(CO). For a given A^H^), 
published models of diffuse and translucent clouds predict less CO than is observed. The J = 
1-3 rotational levels of 12 CO are sub-thermally populated in these sight lines, with T ex typically 
between 3 and 7 K. In general, there appears to be no significant difference between the excitation 
temperatures of 12 CO and 13 CO. Fits to the higher resolution CO line profiles suggest that CO 
(like CN) is concentrated in relatively cold, dense gas. 

We obtain C2 column densities from the F-X (0-0) band at 1341 A (three sight lines; J = 
0-12), the F-X (1-0) band at 1314 A (one sight line; J = 0-12), the D-X (0-0) band at 2313 A 
(four sight lines; J = 0-18), and the A-X (3-0) and (2-0) bands at 7719 and 8757 A (seven sight 
lines; J = 0-12). Comparisons among those column densities yield a set of mutually consistent 
band /-values for the UV and optical C2 bands, but also reveal some apparent anomalies within 
the F-X (0-0) band. Both the kinetic temperature Tk inferred from the C2 rotational populations 
and the excitation temperature Jo2(C2) are generally smaller than the corresponding Toi(H2) — 
suggesting that C2 is concentrated in colder, denser gas than H 2 . 

Incorporating additional column density data for K I, HD, CH, C2, C3, CN, and CO from 
the literature (for a total sample of 74 sight lines), we find that (1) CO is most tightly correlated 
with CN; (2) the ratios 12 CO/H 2 and 13 CO/H 2 both are fairly tightly correlated with the density 
indicator CN/CH (but C2/H2 is not); and (3) the ratio 12 CO/ 13 CO is somewhat anti-correlated 
with both CN/CH and A(CO). Sight lines with 12 CO/ 13 CO below the average Galactic value of 
12 C/ 13 C appear to sample colder, denser gas in which isotope exchange reactions have enhanced 
13 CO, relative to 12 CO. 
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1. Introduction 

Apart from H2 and (in some cases) HD, CO is the most abundant interstellar molecule. Determinations 
of the abundance and rotational excitation of the dominant 12 C 16 and its isotopomers thus can provide 
significant information about the physical conditions in Galactic and extragalactic molecular clouds. The 
CO/H 2 abundance ratio is often assumed to be roughly constant in dense molecular gas, with a "canonical" 
value of order 10 -4 derived primarily from millimetric and radio observations of CO and inferred abundances 
of H 2 in Galactic molecular clouds (Young & Scoville 1991). Under that assumption, measurements of CO are 
used to estimate the total H 2 content when no data are available for H 2 . In more diffuse gas, both H 2 and CO 
were detected directly via their absorption bands in the far-UV using spectra from the Copernicus satellite; 
longer wavelength UV bands of CO were observed with IUE. Analysis of the Copernicus spectra showed 
that the CO/H 2 ratio varied by two orders of magnitude (from about 5 x 1CP 8 to 5 x 1CP 6 ) (Federman 
ct al. 1980) — so that CO is not a good predictor for H 2 in more diffuse material. Theoretical models of 
cloud chemistry suggested that such large variations in the CO/H 2 ratio could be expected in diffuse (Ay 
< 1 mag) and translucent (1 mag < Ay < 5 mag) regions still permeated by dissociating photons (e.g., van 
Dishoeck & Black 1986a, 1988) as well as in regions perturbed by the recent passage of dissociative shocks 
(Neufeld & Dalgarno 1989). The relatively low sensitivity and resolution of Copernicus and IUE, however, 
allowed only limited investigation of more heavily reddened sight lines, where the CO/H 2 ratio might be 
expected to approach the dense cloud value. 

Determining observationally the population distribution in the individual rotational levels of any molec- 
ular species can yield constraints on the physical properties characterizing the molecular gas (e.g., van 
Dishoeck et al. 1991; Federman et al. 1997b; Wannier et al. 1997). Observations of the rotational excitation 
of C 2 , for example, have yielded estimates for the temperatures and densities in diffuse molecular clouds 
(e.g., van Dishoeck & de Zeeuw 1984); data for CN have been used both to constrain the temperature of the 
cosmic microwave background (e.g., Roth & Meyer 1995) and to estimate local hydrogen densities (e.g., Black 
& van Dishoeck 1991). The pervasive abundance of the CO and H 2 molecules in the interstellar gas make 
those species especially important tools for investigations of the gas properties. Unfortunately, however, the 
spectrographs on Copernicus and IUE could not resolve the rotational structure of the CO bands. 

Isotopic ratio measurements provide another means to investigate the chemical processes occurring 
in the interstellar medium. Carbon, whose isotopes 12 C and 13 C are products of primary and secondary 
nucleosynthesis processes (Langer & Penzias 1990), respectively, is well suited for such studies. Within 
about 1-2 kpc of the Sun, the current average 12 C/ 13 C ratio is about 70 (e.g., Wannier et al. 1982; Stahl 
& Wilson 1992; Lucas & Liszt 1998), with some variations; there also is an apparent systematic decrease 
in the ratio toward the Galactic center (e.g., Wilson 1999; Milam et al. 2005). When determined from 
molecules unaffected by chemical fractionation (e.g., CH + ), the 12 C/ 13 C ratio is an indicator of Galactic 
stellar evolution. The (optical) absorption lines from 13 CH + , however, are typically very weak. Because 
chemical fractionation can occur for CO, the 12 CO/ 13 CO ratio gives additional information on chemical 
processes in the molecular gas (e.g., Bally & Langer 1982; Sheffer et al. 1992). Again, however, the limited 
sensitivity of Copernicus and IUE made it difficult to obtain reliable measures of the generally weak 13 CO 
absorption. 

The higher far-UV sensitivity of the FUSE instrument (Moos et al. 2000), together with the higher 
sensitivity and resolution of the HST GHRS and STIS echelle modes in the near-UV, however, has allowed 
us to revisit these issues. The FUSE spectra contain numerous absorption lines from the Lyman and Werner 
bands of molecular hydrogen, from which accurate total H 2 column densities can be derived even toward 
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fainter, more heavily reddened stars (e.g., Rachford et al. 2002). While the B-X (0-0) (1150 A), C-X (0-0) 
(1087 A), and E-X (0-0) (1076 A) bands arising from the A X E + ground electronic state of CO are also present 
within the FUSE wavelength coverage, an accurate determination of the CO column densities from the FUSE 
data alone is often hindered by line blending, saturation, and unresolved structure in the absorption profiles. 
With GHRS and STIS, however, one could observe weaker bands in the CO A-X system between 1240 and 
1550 A at high enough resolution and S/N to discern and measure the individual rotational (J) levels of the 
bands — allowing accurate estimates of their individual column densities. Moreover, some of the weaker 
CO spin-forbidden (triplet-singlet) intersystem bands between 1300 and 1550 A can be detected once the 
permitted bands become saturated (Morton & Noreau 1994) — enabling more reliable determinations of 
higher CO column densities. Several recent studies, for example, have used UV spectra from HST and/or 
FUSE to determine accurate N(CO) toward the moderately reddened stars £ Oph (Lambert et al. 1994), p 
Oph and \ Oph (Federman et al. 2003), X Per (Sheffer, Federman, & Lambert 2002a; Sheffcr, Lambert, & 
Federman 2002b), and HD 203374A (Sheffer, Federman, & Andcrsson 2003); Pan et al. (2005) have reported 
A^H^) and JV(CO) toward a number of stars in Cep OB2 and Cep OB3. High-resolution GHRS spectra of 
the C2 D-X (Mulliken) (0-0) band toward ( Oph have enabled measurements of the C2 rotational levels up 
to J = 24 (Lambert, Sheffer, & Federman 1995) — higher J than had been obtainable from the weaker A-X 
(Phillips) bands in the optical and near-IR. 

In this paper, we present a study of the 12 CO, 13 CO, and/or C2 molecules toward ten moderately 
reddened [E(B - V) = 0.37 0.72] Galactic stars, based on UV spectra obtained with FUSE, IUE, and/or 
HST; on higher resolution optical spectra of CH, CN, and K I; and on very high S/N optical spectra of 
C2. This combination of UV and optical data has allowed us to perform detailed analyses of the CO and 
C2 profiles and to derive column densities for the individual rotational levels of the ground electronic and 
vibrational states — up to J = 3 for 12 CO, J = 2 for 13 CO, and J = 18 for C2. In § 2, we discuss the 
observational data and the data reduction procedures. In § 3, we describe the methods used to determine 
column densities from the spectra. In § 4, we combine our new results with existing literature data for 
H 2 , HD, CH, C2, C3, CN, and CO (for a total sample of 74 stars) and discuss the physical conditions 
characterizing the diffuse molecular gas using this expanded star sample. Column densities of H2, CH, C2, 
CN, 12 CO, and 13 CO for that total sample are tabulated in the Appendix. In § 5, we summarize the main 
results of this study. 



2. Observations and Data Reduction 

In order to investigate the behavior of the CO/H2 and 12 CO/ 13 CO ratios at higher overall extinctions 
and column densities, we selected sight lines according to the following general criteria: (1) inclusion in the 
FUSE survey of translucent sight lines (Rachford et al. 2002; for H2 and CO); (2) availability of higher 
resolution HST GHRS or STIS spectra of the CO A-X bands (for 12 CO, 13 CO, and the CO rotational 
structure); and (3) availability of high-resolution optical spectra of K I, CH, and CN (for information on 
component structure). For the sight line toward HD 73882 (Snow et al. 2000), which was not observed with 
HST, and for several other sight lines with HST data, we also have analyzed archival high-dispersion IUE 
spectra of the CO A-X and intersystem bands. The ten sight lines in the resulting sample (Table 1) have 
E(B — V) between 0.37 and 0.72 mag, visual extinction Ay between 0.83 and 2.36 mag, A(H2) from 3-13 
x 10 20 cm~ 2 , and molecular fractions /(H 2 ) between 0.14 and 0.76. The ten sight lines also probe a variety 
of regions in the Galactic interstellar medium — with two (HD 24534, HD 27778) in Taurus/Perseus, one 
(HD 147888) in Sco-Oph, three (HD 206267, HD 207198, HD 210839) in Cep OB2, and one (HD 210121) 
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through a high-latitude molecular cloud. 

2.1. FUSE Data 

FUSE observed our sample stars under programs P116 and X021 (PI: T. P. Snow), which were dedicated 
to the study of translucent sight lines (Rachford ct al. 2002). The exposures used for each sight line are 
listed in Table 2. The observations were obtained in time tag mode through the low- resolution aperture 
(LWRS) and were processed with version 2.0.5 of the CALFUSE pipeline (Kruk & Murphy 2001; Sahnow 
ct al. 2000). Eight detector segments cover sections of the total wavelength range from 912 to 1187 A, with 
most spectral regions covered by two to four segments. The resolution is ^0.062 A (FWHM), corresponding 
to about ^18 km s _1 or ~9 pixels. For each segment, the exposures were co-added after cross-correlating 
and shifting the spectra (by at most 8 pixels) with respect to the brightest exposure. The spectra were 
binned by 4 pixels (slightly less than one-half resolution element) in order to increase the signal-to-noise 
(S/N) ratio without degrading the resolution delivered by the spectrograph optics. The co-added spectra 
typically have S/N of about 10 per resolution element in the LiFlA spectrum and up to 15 per resolution 
element in the LiF2A spectrum. The LiFIB spectrum was excluded from the data analysis because of the 
presence of a well-known detector artifact which causes an artificial flux deficiency in this segment (Kruk & 
Murphy 2001). Figure 1 presents the profiles of the CO E-X (0-0) (1076 A), C-X (0-0) (1087 A), and B-X 
(0-0) (1150 A) bands, as functions of heliocentric velocity, detected in the far-UV spectra of seven of the ten 
sight lines. 

2.2. HST Data 

The calibrated GHRS or STIS spectra listed in Table 2 were retrieved from the MAST archive located 
at the Space Telescope Science Institute. All the STIS data were automatically reduced with the latest 
version of the CALSTIS pipeline. The STIS spectra were obtained using either the FUV MAMA detector 
and E140H cchelle mode or the NUV detector and E230H echelle mode, at resolutions of either 2.75 km s _1 
(for the 0."2 x 0."2 or 0."2 x 0."09 slits) or 1.5 km s" 1 (for the 0."1 x 0."03 slit). The E140H spectra 
generally cover the wavelength range from about 1150 to 1370 A — including the A-X (7-0) through (12-0) 
bands of 12 CO and 13 CO (between 1240 and 1350 A) and the F-X (0-0) and (1-0) bands of C 2 at 1341 A 
and 1314 A (as well as many atomic absorption lines). Unfortunately, the E140H spectra generally do not 
include the strongest CO intersystem bands. The more extensive E140H spectra of HD 24534, which do 
include those intersystem bands, are described by Sheffer et al. (2002a, 2002b). The E230H spectra cover 
the range from about 2130 to 2400 A — including the D-X (0-0) band of C 2 at 2313 A (and various atomic 
lines). The extracted spectral orders were corrected for scattered light contamination using the procedures 
described by Howk & Scmbach (2000) or Valenti et al. (2002). For the spectra of HD 206267 and HD 210839, 
which were obtained through the narrowest slit, additional processing was performed by E. Jenkins (private 
communication) in order to account for unbalancing of the even and odd pixels on the detector (see Jenkins 
& Tripp 2001 for details). The continua were normalized to unity using low-order polynomials. In most 
cases, the S/N ratios (estimated from the fluctuations in the continuum regions) range from 16 to 36 per 
resolution element; the C 2 F-X band spectra for HD 24534, however, have S/N ~ 130, and the C 2 D-X band 
spectra (for three stars) have S/N ~ 85-120. Figures 2-4 show the profiles of some of the 12 CO A-X bands 
toward six of the stars; Figure 5 shows the 13 CO A-X (7-0) (1347 A) band profiles toward four of the stars 
and the 13 CO A-X (4-0) (1421 A) band toward HD 24534. The individual rotational levels (J = 0, 1, 2, 3) 
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of the ground vibrational state resolved in each CO band are numbered above the spectra; for CO, the R, 
Q, and P branches are left-to-right in each case. 

The GHRS ECH-A spectra of HD 210121 cover the 12 CO A-X (2-0) (1477 A) and (4-0) (1419 A) bands, 
the corresponding 13 CO bands (1478 and 1421 A), and the relatively strong 12 CO aT4 intersystem band 
(1419 A), at a resolution of about 3.5 km s _1 . The GHRS ECH-B spectra of HD 24534 cover the C 2 D-X 
(0-0) band at 2313 A. Because the spectra were obtained at several fp-split offsets, the STSDAS routines 
POFFSETS, DOPOFF, and SPECALIGN were used to align and combine the individual spectra. The S/N 
values range from 15 to 20 per resolution element for HD 210121 and are about 100 per resolution element 
for HD 24534. 



2.3. IUE Data 

High-dispersion (FWHM - 25 km s" 1 ) IUE spectra of HD 27778, HD 73882, and HD 207198 were 
retrieved from the archive (Table 2). For all three stars, multiple spectra had been obtained at several 
different offsets within the large aperture, in order to reduce the effects of fixed-pattern noise when the 
spectra are combined. For these three stars, the S/Ns range from about 15-45 per resolution element in the 
aligned, co-added spectra, for wavelengths between about 1320 and 1550 A [covering the CO A-X bands (0-0) 
through (8-0)]. Figure 6 shows some of the stronger 12 CO A-X bands toward HD 27778 and HD 73882 (see 
also Figure 7 of Morton & Noreau 1994); the corresponding 13 CO A-X bands and some of the intersystem 
bands of 12 CO can also be discerned. While the permitted A-X bands are saturated toward both of those 
stars, the stronger intersystem bands toward HD 73882 indicate a higher iV(CO) in that line of sight. The 
IUE data for HD 210121 were described by Welty & Fowler (1992). 

2.4. Optical Data 

High-resolution (FWHM ~ 0.6-3.6 km s" 1 ) spectra of K I (7698 A), CH (4300 A), and CN (3874 A) were 
obtained with the Kitt Peak coude feed telescope and the McDonald 2.7m telescope. A detailed discussion 
of the reduction and analysis of the K I spectra has been given by Welty & Hobbs (2001); the CH and CN 
spectra, processed in similar fashion, will be discussed by Welty, Snow & Morton (in preparation). Profiles 
of the optical lines for six of the sight lines discussed in this paper are included in Figures 2-4; some have 
also been shown by Pan et al. (2004). 

Medium resolution (FWHM ~ 8 km s _1 ) but very high S/N (> 1000) optical spectra were acquired for 
nine of the stars with the 3.6m telescope and the ARC echelle spectrograph at Apache Point Observatory as 
part of an observing program designed to investigate the nature of the diffuse interstellar bands. The details 
on the reduction and analysis of these spectra and, in particular, of the C2 A-X (Phillips) bands detected in 
seven of the sight lines studied here, are described in Thorburn et al. (2003). 

3. Data Analysis 

Several factors have made it difficult to derive accurate CO column densities for more heavily reddened 
sight lines. The lower resolution UV spectra from Copernicus, IUE, and FUSE yield equivalent widths 
for a number of the CO bands, but do not resolve either the (often complex) component structure or the 
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rotational structure of the bands. Blends with stellar lines can also complicate the interpretation of the 
spectra, especially when the S/N and/or vsini are relatively low. In addition, for these higher column 
density sight lines, many of the permitted CO bands that can be detected are significantly saturated. 

Fortunately, the higher resolution (and often higher S/N) spectra obtained with the HST GHRS and 
STIS echelle modes at least partially resolve the CO rotational structure and enable reliable measurement 
of weaker CO A-X and intersystem bands even toward fainter, more heavily reddened targets. Agreement 
between several recent observational and theoretical studies suggests that the /-values now available for 
the intersystem bands are reliable (Sheffer et al. 2002a; Eidelsberg & Rostas 2003). Finally, still higher 
resolution optical spectra of K I, CH, and (especially) CN can provide useful information as to the velocity 
component structure — enabling more accurate modeling of the observed CO line profiles. 

Three methods were used to estimate the interstellar CO column densities. First, empirical curves 
of growth were constructed using the measured equivalent widths of the CO bands. The empirical curves 
were then compared with theoretical curves generated for single components with a range in the Doppler 
broadening parameter (6), for two representative CO rotational population distributions characterized by 
different excitation temperatures. Second, the Apparent Optical Depth method (Hobbs 1971; Savage & 
Scmbach 1991) was applied to the observed CO line profiles. Third, the absorption-line profiles were fitted 
using multi-component models derived from higher resolution optical spectra of K I, CH, and CN. Profile 
fitting was used to determine column densities from the UV and optical lines of C2. Because the individual 
gas components along these sight lines are at best partially resolved in the UV spectra, determination of 
the column densities for individual components is often uncertain with STIS data and essentially impossible 
with the lower resolution FUSE data. We will therefore only discuss total sight line column densities in the 
present work. 



3.1. Curves of Growth 

In order to construct the empirical curves of growth (COG), total equivalent widths of the 12 CO and 
13 CO bands seen in the normalized UV spectra were measured over velocity ranges which include all the 
detected rotational lines. Table 3 lists the CO equivalent widths found for the A-X, B-X, C-X, and E-X bands 
from FUSE and STIS data (eight sight lines); Table 4 gives the values for the A-X and intersystem bands 
found from GHRS, STIS, and/or IUE data (six sight lines). The CO equivalent widths for the well-studied 
sight lines toward X Per (HD 24534) and ( Oph (HD 149757) are included in Table 4 for comparison. In 
all cases, the listed \-a uncertainties include contributions from both photon noise and continuum-fitting 
uncertainty, added in quadrature; upper limits are 3-cr. Because of the wide velocity range covered by the 
ensemble of rotational lines, the continuum-fitting uncertainties generally are dominant. Inspection of the 
equivalent widths for HD 27778 and HD 207198 indicates that there is generally good agreement between 
the values determined from IUE and STIS spectra for the CO A-X (7-0) and (8-0) bands at 1322 and 1344 
A. The equivalent widths for HD 185418 and HD 192639 are consistent with those listed by Sonnentrucker 
et al. (2002, 2003). Because the CO C-X (0-0) band at 1087 A is often blended with an adjacent CI I 
line (see Fig. 1) and, in the B stars in our sample, with a stellar absorption feature, it was not used in the 
construction of the curves of growth. An estimate of its equivalent width is, nevertheless, given for the two 
stars toward which the blending was less severe. The E-X (0-0) band at 1076 A is in the wing of the very 
strong H 2 (2-0) R(0) line at 1077 A, making the appropriate continuum difficult to estimate in some cases. 

The rest wavelengths and oscillator strengths listed in Tables 3 and 4 were taken from Morton & Noreau 
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(1994) for the CO A-X bands and from Federman et al. (2001) for the far-UV B-X (0-0) (1150 A), C-X 
(0-0) (1087 A), and E-X (0-0) (1076 A) bands. For the CO intcrsystem bands, we have adopted the rest 
wavelengths of Eidelsberg & Rostas (2003) and (where possible) the empirical /-values of Sheffer et al. 
(2002a), which yield smoother curves of growth for the equivalent widths measured toward both HD 24534 
(X Per) and HD 149757 (( Oph) than do the corresponding /-values of Eidelsberg & Rostas. For the 12 CO 
bands with both FUSE and STIS data, the range in /A can span nearly three orders of magnitude from the 
weakest A-X (12-0) (1246 A) to the strongest E-X (0-0) (1076 A) band; a similar range in /A is covered for 
sight lines with both STIS and IUE data. 

In principle, the CO curve of growth depends on (1) the rotational structure of the individual bands, (2) 
the excitation temperatures characterizing the various rotational levels, and (3) the interstellar component 
structure (e.g., Black & van Dishoeck 1988). Theoretical curves of growth were generated for single interstel- 
lar components with a range of b- values (0.3-5.0 km s _1 — as an approximation to more complex component 
structures), using the CO rotational structure of the A-X (6-0) band (taken as representative of all the A-X 
bands). Separate curves were generated for the far-UV E-X, C-X, and B-X bands, for which the rotational 
structure is somewhat different. For 12 CO, curves were generated for two sets of relative populations in the 
J = 0-3 rotational levels. The first had level populations 0/1/2/3 = 1.0/1.0/0.2/0.05, which are similar to 
the average ratios found for the sight lines in our small sample. The relative populations in the J = 0-2 
levels, for example, correspond to an excitation temperature of about 5 K. The second had level populations 
0/1/2/3 = 1.0/0.5/0.02/0.0001, corresponding to an excitation temperature of about 3 K (typical of some 
other sight lines reported in the literature). For 13 CO, equal populations were assumed for J = and 1 
(the only levels generally detected in the STIS spectra) . The sums of the /-values for the various individual 
transitions (R, Q, and P branches) for each J are all equal to the /-value of the band as a whole, and the 
COG analysis of the integrated band equivalent widths yields the total column density over all J. 

Columns (2) and (7) in Table 5 list the total 12 CO and 13 CO column densities (respectively) derived 
by comparing the empirical curves of growth with the theoretical single-component curves for T cx ~ 3 K 
and/or for T cx ~ 5K; column (3) lists the best-fit "effective" b- value. In each case, the uncertainty given for 
the column density corresponds to the range in AT(CO) giving rms deviations less than 1.5 times that of the 
best (N,b) pair, for that best-fit b- value. Figure 7 shows the 12 CO curves of growth for the sight lines toward 
HD 27778 and HD 207198, with the adopted values for TV and b. In the figure, the solid lines show the T ox ~ 
5 K theoretical curves for the A-X bands (which in most cases should be adequate for the intersystem, B-X, 
and C-X bands); the dotted lines show the theoretical curves for the E-X (0-0) band at 1076 A, displaced 
downward by 1.0 for clarity. The data point for the observed E-X (0-0) band is plotted against both sets of 
theoretical curves to show that the curves for the E-X (0-0) band lie slightly below those for the A-X bands. 

For most of the sight lines, at least two CO bands lie on or near the linear part of the curve of growth. 
Accurate total column densities can be derived from such optically thin bands with minimal assumptions 
regarding the gas velocity distribution. When additional, stronger bands are measured, the COG analysis 
can (in principle) also yield some information regarding the overall gas distribution, as long as (1) the /-value 
range for the stronger bands is large enough that different degrees of saturation are observed and (2) the CO 
excitation temperature is reasonably well determined. For eight of the nine sight lines analyzed, the effective 
6-values for the empirical curves, obtained by comparison with the theoretical curves, range from about 0.5 
to 2.0 km s _1 . The very small effective 6-values for HD 147888 and HD 210839 (0.5 km s^ 1 ) likely reflect 
the widths of the single dominant narrow components seen in CN in those two sight lines. The higher values 
for HD 27778 and HD 206267 (both - 1 km s" 1 ) and for HD 207198 (~ 2.0 km s" 1 ), however, presumably 
correspond to the multiple components seen in CN in those sight lines (rather than the widths of single 
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components). Slight differences between the empirical and theoretical curves for several sight lines may be 
due to the more complex component structure (seen in the higher resolution optical spectra) in those cases. 
Comparisons of the observed equivalent widths with the theoretical curves for T C x ~ 3 K (which fall slightly 
below the corresponding curves for T ox ~ 5 K) generally exhibited the best agreement for somewhat higher 
column densities at the same effective b or else for similar column densities and slightly higher effective b. 
While the column densities derived via the COG thus depend on knowledge of T cx (obtained via detailed 
fits to the high-resolution line profiles), they do provide useful confirmation of the total N(CO) derived in 
the profile fits. 

3.2. Apparent Optical Depth 

The Apparent Optical Depth (AOD) analysis was applied to all detected 12 CO and 13 CO bands, over 
the same velocity ranges used for the equivalent width measurements. In the absence of unresolved saturated 
structure(s) within the line profile, the AOD method yields the "apparent" (instrumentally smeared) column 
density as a function of velocity, with some improvement over the simplest assumption of a linear relationship 
between equivalent width and column density. If the line is somewhat saturated, the AOD analysis yields a 
lower limit to the actual column density, especially if the spectral resolution is not very high. In principle, 
the degree of saturation in the profiles of a given species can be assessed if multiple lines, with /-values 
differing by at least a factor of 2, are available (Jenkins 1996). For CO, the distribution of the total column 
density over multiple rotational levels (some blended) renders the specific N(v) not very meaningful, but it 
does delay somewhat the onset of saturation. As for the COG analysis, the total column density (over the 
whole profile and for all rotational levels) can still be estimated. 

For most of the sight lines, smaller total column densities are obtained from the stronger CO bands 
(relative to those obtained from the weakest observed bands) — indicative of more severe saturation in the 
stronger bands. Toward HD 207198, for example, the AOD analysis of the strongest CO band observed with 
STIS [A-X (7-0) at 1344 A] yields a total N(CO) smaller by a factor of about 1.5 than that found for the 
weaker CO bands. Toward HD 210839 and HD 206267, saturation effects are already apparent in the A-X 
(9-0) (1301 A) band, and the total column densities inferred from the A-X (7-0) band at 1344 A are factors 
of 3-4 lower than those obtained from the weakest bands. In such cases, however, accurate column densities 
still can be determined by using only the weaker A-X (12-0), (11-0), and (10-0) bands (at 1246, 1263, and 
1281 A), which yield concordant results. The AOD column densities inferred from the weakest CO bands are 
listed in Columns (3) and (7) of Table 5; the associated error bars include contributions from both photon 
noise and continuum-fitting uncertainty. In most cases, the total column densities derived from the AOD 
analysis are consistent with those obtained with the COG method. 

3.3. Profile Fitting 

While the STIS resolution is high enough to at least partially resolve the individual rotational levels 
in the CO and C2 band structures (Figs. 2-5 and 8-10), it is not sufficient to separate or resolve all of 
the multiple, narrow gas components that produce the observed absorption lines. The blending between 
components — both interstellar and rotational — is particularly severe for the shorter wavelength R-branch 
lines of the CO A-X bands (Figs. 2-4) and the C2 F-X bands (Fig. 10). In order to determine accurate (total) 
column densities for the individual rotational levels contributing to the molecular bands, the absorption-line 
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profiles were fitted with multi-component models based on fits to higher resolution optical spectra. The 
profile fitting program "Owens" (Lemoine et al. 2002) and a variant of the program FITS6P (e.g., Welty, 
Hobbs, & Morton 2003) — both of which perform iterative least-squares fits to the observed line profiles — 
were used for the detailed component analyses of the optical and UV spectra. 

While lines of molecular species such as CH, C2 , and CN are the most direct optical tracers of diffuse 
molecular gas, the good correlations observed among the column densities of H 2 , CH, and K I — and the 
often striking correspondence between the line profiles of CH and K I — suggest that component information 
from K I may be relevant as well (Welty & Hobbs 2001). More detailed structure is often discernible in the 
K I line profiles because (1) the potassium atom is heavier than the CH molecule, (2) the hyperfine splitting 
of the K I A7698 line (~ 0.3 km s _1 ) is smaller than the separation between the lambda-doubled components 
of the CH A4300 line (1.43 km s -1 ), and (3) the K I line is typically stronger. The more detailed K I profile 
structures (i.e., the relative velocities and 6- values for each discernible component) were therefore used to 
model the corresponding CH profiles. For the ten sight lines in this study, between three and thirteen 
components were needed to fit the K I profiles; not all were detected in CH, however. Because CN is likely 
to be more concentrated in the colder, denser regions of each sight line, independent fits were performed 
to the CN profiles; one to four components (typically corresponding to the main components seen in K I 
and CH) were required. The structures found from the high-resolution optical spectra of K I, CH, and/or 
CN (Appendix A; Table 14) were then taken to represent the component distributions for the UV CO and 
C2 bands in each sight line. Knowledge of the component structure can be crucial for determining column 
densities from the stronger CO lines, but is generally less important for the corresponding analyses of the 
much weaker C2 lines. 



3.3.1. CO Column Densities 

The component structures found for K I and CH were used for the initial fits to the CO band profiles 
observed with GHRS and STIS. In the Owens fits, the relative velocity and the b- value of each component 
were fixed, while the continuum placement and individual component column densities were allowed to vary. 
The FITS6P fits were performed on the normalized spectra, fixing the relative velocities (usually) and the 
6-values but allowing the component column densities to vary. In both cases, all the CO A-X band profiles 
were fitted simultaneously. For sight lines where both Owens and FITS6P were used to analyze the CO 
bands, the two methods generally yielded consistent results (within the respective uncertainties) both for 
the column densities of the individual rotational levels (J = 0-3) and for the total A^(CO). 

Toward HD 192639, HD 185418, and HD 207198, all the CO bands were adequately fitted using the 
K I component structures. Toward HD 27778, HD 206267, and HD 210839, however, fits using the K I 
structures tended to under-produce the weakest bands and overproduce the strongest bands. In an attempt 
to obtain mutually consistent column densities for all the bands in those sight lines, two additional scries 
of fits were undertaken. In the first additional series, the relative velocities of the various gas components 
were again fixed to those of K I, but both the column densities and the 6-values were allowed to vary. For all 
three sight lines, the best fits to the CO bands were realized with 6-values between 0.3 and 0.7 km s _1 for 
the main components — i.e., somewhat smaller than the 0.6-1.1 km s _1 values found for K I. The second 
additional series of fits, using the (simpler) CN velocity structures (but with slightly smaller 6-values), also 
yielded more consistent fits to the ensemble of observed profiles. Inspection of the observed line profiles 
suggests that the profiles of the individual CO rotational lines seen in the higher resolution STIS spectra 
of HD 206267 and HD 210839 generally are more similar to the profiles of CN than to those of K I and 
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CH. A close correspondence between the component structures for CO and CN has also been found by Pan 
et al. (2005) toward a number of stars in the Cep OB2 and Cep OB3 associations, including HD 206267, 
HD 207198, and HD 210839. For those three lines of sight, our total CO column densities are consistent 
with those determined by Pan et al. 

Examples of the fits to some of the CO A-X system profiles toward HD 27778 and HD 210839 are shown 
in Figures 8 and 9, respectively, together with the corresponding higher resolution optical profiles of K I, 
CH, and CN. The fits to the profiles for HD 27778 utilize two components separated by 2.3 km s _1 , with 
b <~ 0.5-0.6 km s _1 and relative column densities ~ 3:1 — i.e., similar to the structure found for CN, but 
with slightly smaller b- values and somewhat different relative strengths (Table 14). That structure (with the 
rotational level column densities Nj obtained from fitting the ensemble of lines) provides acceptable fits to 
both weak and strong CO bands, including the stronger A-X bands observed at lower resolution with IUE. 
The best fits to the profiles for HD 210839 employ two components separated by 2.1 km s _1 , with b ~ 0.4 
km s _1 (slightly smaller than the 6-values adopted for CN) and relative strengths roughly 3:1. 

Because the lower resolution IUE spectra of HD 27778, HD 73882, HD 207198, and HD 210121 do not 
resolve the CO rotational structure, the line profiles and total equivalent widths of the various permitted 
A-X and intersystem CO bands listed in Table 4 were fitted by assuming rotational excitation temperatures 
of 3 and 5 K, using the component structures determined for K I and CN. The resulting total sight line 
CO column densities toward HD 27778 and HD 207198 are consistent with those derived from the higher 
resolution STIS spectra, but differ somewhat from the values determined by Joseph et al. (1986) (and 
Federman et al. 1994a) from curve of growth analyses of more limited IUE spectra of the permitted A- 
X bands (only). While no high-resolution CN spectra are available for HD 73882, three components are 
seen in both K I absorption and mm- wave CO emission (van Dishoeck et al. 1991). Three-component fits 
to the intersystem and weakest permitted A-X bands of 12 CO, assuming T ox = 5 K and several choices 
for the relative column densities in each component, yield column densities consistent both with those 
determined from the COG analysis and with the value estimated from the CO emission lines (van Dishoeck 
et al. 1991). Because the higher resolution GHRS spectra of HD 210121 allow better recognition of narrow 
stellar absorption features blended with several of the interstellar CO lines, the column density of 13 CO is 
significantly smaller than that estimated by Welty & Fowler (1992) from the IUE spectra alone. Adoption 
of a smaller CO 6-value (0.5 km s _1 vs. 1.0 km s _1 ) for that sight line, however, yields a somewhat higher 
column density for 12 CO. 

The results of the CO profile fitting (FIT) analyses are summarized in Table 5, which lists the total 
12 CO and 13 CO column densities (for all gas components and all rotational levels), and in Table 6, which 
lists the column densities Nj for the individual rotational levels (J = 0-3) for the eight sight lines with high- 
resolution UV spectra. Where only Nq has been measured for 13 CO, we have assumed that No/Ni (or the 
excitation temperature Tqi) is the same as for 12 CO. The uncertainties in the Nj include contributions from 
photon noise, continuum uncertainty, and uncertainties in the component 6-values (± 0.05-0.10 km s _1 ). 
The continuum uncertainties are most significant for the weaker lines; the 6-value uncertainties generally 
are most significant for the stronger lines (especially for b < 0.5 km s _1 ). While the column densities 
derived from the AOD, COG, and FIT methods agree within the mutual uncertainties, we have adopted the 
N(CO) determined in the profile fits, as they account explicitly for saturation (using the detailed component 
structures), they use information from all the observed transitions, and they yield values for the individual 
rotational levels. The better fits to the CO bands obtained both with smaller 6-values and with the CN 
component structure suggest that CO (like CN) may preferentially sample colder, denser gas than K I and 
CH (see also Pan et al. 2005). Discussions of the physical conditions in the gas in the following sections 
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seem consistent with that interpretation. 



3.3.2. C2 Column Densities 

Nearly all of the existing studies of interstellar C2 have been based on optical and/or near-IR spectra 
of the A-X (Phillips) bands [(1-0) at 10144 A, (2-0) at 8757 A, and (3-0) at 7719 A] (e.g., Souza & Lutz 
1977; Hobbs & Campbell 1982; van Dishoeck & de Zeeuw 1984; van Dishoeck & Black 1986b; Federman 
ct al. 1994a; Gredel 1999). In one exception, high-resolution GHRS ECH-B spectra of the D-X (Mulliken) 
(0-0) band at 2313 A and lower resolution G160M spectra (FWHM - 16 km s" 1 ) of the F-X (0-0) band at 
1341 A toward C Oph were used to obtain the C2 rotational level populations up to J = 24 and to estimate 
the relative /-values for the C2 A-X, D-X, and F-X systems (Lambert, Shcffer, & Federman 1995). Many of 
the lines from the individual rotational levels in the F-X (0-0) band are inextricably blended in the G160M 
spectrum, however. Kaczmarczyk (2000) has discussed GHRS echelle spectra of the F-X (0-0) and D-X (0-0) 
bands toward HD 24534. 

In this paper, we present high-resolution profiles of the UV C2 F-X (0-0) and (1-0) bands (as seen in the 
STIS E140H spectra of HD 24534, HD 27778, and HD 206267; Figures 10 and 11), along with high-resolution 
GHRS ECH-B and STIS E230H spectra of the D-X (0-0) band toward HD 24534, HD 27778, HD 147888, 
and HD 207198 (Figure 12) and high S/N optical spectra of the A-X (3-0) and (2-0) bands toward seven 
stars (Figures 13 and 14). All these UV and optical spectra at least partially resolve the rotational structure 
of the C2 bands — enabling detailed line profile analyses of the C2 absorption. The D-X (0-0) band is 
particularly useful, as the individual rotational lines from both the R and P branches are well separated 
and relatively strong. Equivalent widths for the individual D-X (0-0) lines are given in Table 7; equivalent 
widths for lines in the A-X (2-0) and (3-0) bands are listed in Table 8. Wavelengths and relative strengths 
for the individual rotational lines in the various C2 bands — as well as the adopted total band /-values — 
are listed in Table 15 (Appendix B). 

By comparing column densities derived from the various UV and optical bands, a set of mutually consis- 
tent band oscillator strengths may be determined. Lambert et al. (1995) adopted a band oscillator strength 
/ = 0.0545 for the D-X (0-0) band, based on two independent (and closely concordant) theoretical estimates. 
Comparison of the overall strengths of the D-X (0-0) and F-X (0-0) bands toward Q Oph then yielded / 
= 0.10±0.01 for the latter band — in excellent agreement with the theoretical value 0.098 determined by 
Bruna & Grcin (2001). With those /-values, fits to the profiles of the D-X (0-0) and F-X (0-0) bands to- 
ward HD 24534 and HD 27778 yield reasonably consistent C2 column densities for rotational levels J = 0-6 
(Table 9), though there appear to be some systematic discrepancies for J — 8 and 10 (see below). Similar 
fits to the F-X (1-0) band at 1314 A toward HD 24534 then suggest a total oscillator strength / ~ 0.06 for 
that band. 

As for CO, the K I velocity structure was initially adopted to determine the column densities of the 
individual C2 rotational levels present in the spectra of the three UV bands. Additional fits were performed 
using either a single component (with an effective b) or the structure that best accounted for the saturation in 
the stronger CO bands; a separate three-component fit was made to the D-X (0-0) band toward HD 207198, 
where more complex structure (similar to that seen in K I and CH) is evident (Figs. 4 and 12). In each 
case, the various fits led to similar results (within the uncertainties), indicating that the total C2 column 
densities derived from the C2 D-X (0-0) and F-X (0-0) and (1-0) bands generally are not very sensitive to 
the details of the gas velocity distribution along these lines of sight. Column densities for individual C2 
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rotational levels up to J — 12 were thus obtained from the F-X band(s) toward HD 24534, HD 27778, and 
HD 206267 and up to J = 18 from the D-X (0-0) band toward HD 24534, HD 27778, HD 147888, and 
HD 207198 (Table 9). The listed uncertainties include contributions from photon noise, continuum- fitting 
uncertainty, and uncertainty in the effective b; the uncertainties for the D-X (0-0) band are similar to those 
that would be inferred from the uncertainties in the equivalent widths of the individual rotational lines. The 
relative rotational level populations Nj/N 2 appear to be fairly similar for these five lines of sight (but see 
the caveats below). Toward HD 210839, the C2 F-X (0-0) band absorption is significantly weaker, so that 
detailed rotational distributions could not be determined from the UV spectra. 

Detailed examination of the fits to the various C2 band profiles and comparisons of the C2 column 
densities obtained from those bands, however, indicate that the strength and structure of the UV C2 bands 
— particularly the F-X (0-0) band — are not completely understood. On the one hand, the fits to the D-X 
(0-0) and F-X (1-0) bands toward HD 24534 are quite good (Figs. 12 and 10) and yield mutually consistent 
column densities for the individual rotational levels J = 0-10. Moreover, while there are some "discrepant" 
Nj for individual sight lines, there do not appear to be any systematic differences in the Nj determined from 
the D-X (0-0), A-X (3-0), and A-X (2-0) bands (Table 9). The fits to the F-X (0-0) band toward HD 24534, 
HD 27778, and HD 206267, however, all are rather poor near the "pile-up" of the R(6)-R(12) lines; some 
additional broad (?) absorption appears to be present there. In addition, the calculated wavelengths of the 
clearly detected F-X (0-0) lines (J < 10) seem to be too small by about 36 mA (8 km s _1 ). Finally, the 
lines from J = 6-10 seem to be weaker in the F-X (0-0) band than in the intrinsically weaker F-X (1-0) 
band; compare, for example, the relatively unblended Q(8), Q(10), and P(8) lines of the two bands toward 
HD 24534 in the upper half of Figure 10. 

The discrepancies between the C2 D-X (0-0) and F-X (0-0) bands previously noted in GHRS spectra of 
£ Oph (Lambert et al. 1995) and HD 24534 (Kaczmarczyk 2000) may reflect similar effects. Lambert et al. 
might not have recognized either the "extra" absorption near the F-X (0-0) band R-branch "pile-up" or the 
8 km s _1 velocity offset in the lower resolution G160M spectrum. Their fits to the whole F-X (0-0) band, 
which uniformly scaled the rotational level column densities determined from the D-X (0-0) band at 2313 A, 
would thus have tended to predict too much absorption over the rest of the band — as seen, for example, 
for the Q(10)-P(4), Q(12)-P(6), and Q(14)-P(8) blends in their Figure 3. Kaczmarczyk found that the Nj 
derived from the D-X (0-0) band did not yield a good fit to the Q(12)-P(6) and Q(14)-P(8) blends in the 
F-X (0-0) band and noted (unspecified) problems with the F-X (0-0) band wavelengths. While the "extra" 
absorption near the F-X (0-0) band R-branch "pile-up" is not evident in the lower S/N GHRS spectrum 
analyzed by Kaczmarczyk, he did mention a "lower than expected" continuum near several of the stronger 
lines. 

The causes of the discrepant behavior observed for the F-X (0-0) band are not understood. In principle, 
the "extra" broad absorption near the R-branch pile-up could be a stellar feature, but the three stars 
observed (HD 24534, HD 27778, and HD 206267) have somewhat different spectral types, radial velocities, 
and projected rotational velocities. Moreover, the "extra" absorption toward HD 24534 is very similar in 
three separate sets of STIS spectra taken at epochs separated by about 3 and 33 days, respectively. The 
wavelengths adopted for both of the F-X bands (Table 15) were calculated from molecular constants derived 
from flash discharge spectra analyzed by Herzberg, Lagerqvist, & Malmberg (1969). The 8 km s _1 velocity 
offset for the clearly detected (and unblended) lines in the F-X (0-0) band (all with J < 10) does not appear 
to be due to errors in the STIS wavelengths, as the (0-0) band appears in two adjacent orders, one of which 
also contains the 12 CO A-X (7-0) band at 1344 A and the CI I line at 1347 A. While the velocity offset for 
those lower J F-X (0-0) band lines would be consistent with a slightly smaller value for the upper level i/q 
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(74530.9 cm -1 instead of 74532.9 cm -1 ), the wave numbers for the higher J lines measured by Herzberg 
ct al. (1969) generally are consistent (within ± 0.2 cm -1 ) with the values predicted using their tabulated 
constants. Unfortunately, Herzberg et al. apparently could not measure the positions of any of the lower 
J lines in the F-X (0-0) band, and there are no clear detections of any of the unblended higher J lines in 
the spectra shown in Figure 10. There are, however, tantalizing hints of weak absorption at the predicted 
(i.e., unshifted) positions of several of the higher J Q-branch lines in the spectrum of HD 24534 [which was 
taken from data set o66p01020, the longest single STIS exposure of HD 24534 covering the F-X (0-0) band]. 
When STIS spectra of HD 24534 obtained at two other epochs are combined with the spectrum shown in 
Figure 10 (more than doubling the total exposure time), those weak features persist. Fits to that combined 
spectrum with only the lines from J < 10 shifted yield low values of Nj for J = 8-12, but quite reasonable 
values of Nj for J = 14-18 (Figure 11, Table 9). 

Column densities for the individual C2 rotational levels were also obtained from the A-X (2-0) and 
(3-0) bands at 8757 and 7719 A (Thorburn et al. 2003), using both the measured equivalent widths of 
the individual rotational lines in the medium-resolution ARCES spectra and fits to the ensemble of lines 
(including blended features). Figure 13 shows normalized spectra of the optical C2 A-X (2-0) band toward 
seven stars; Figure 14 shows corresponding spectra of the A-X (3-0) band toward five stars. The profiles 
observed toward HD 204827, which has high column densities of both C2 and C3 (Thorburn et al. 2003; Oka 
et al. 2003; Adamkovics et al. 2003) are included to show the band structures more clearly. Table 8 lists 
the wavelengths, log(/A) values, and equivalent widths for the A-X (2-0) and (3-0) band transitions toward 
seven of our stars (plus HD 204827); C2 was not detected toward HD 185418 or HD 192639 and was not 
observed toward HD 73882 (though see Gredel et al. 1993). The equivalent widths of A-X (2-0) R-branch 
lines located in the steep wing or deep core of the stellar Paschen 12 line of H I can be rather uncertain; in 
some cases (e.g., toward HD 24534), some of those R-branch lines could not be reliably measured. There is 
also a weak diffuse interstellar band nearly coincident with the A-X (3-0) Q(2) line in several cases (Herbig 
& Leka 1991). The equivalent widths in Table 8 are consistent with most of the (generally lower precision) 
values previously reported for HD 27778 (Federman et al. 1994a); HD 24534, HD 206267, HD 207198, and 
HD 210839 (Federman & Lambert 1988); and HD 210121 (Gredel et al. 1992). 

Unfortunately, the band oscillator strengths for the C2 A-X bands have been somewhat uncertain, with 
persistent differences between theoretical and experimental estimates (e.g., van Dishoeck & Black 1989; 
Langhoff et al. 1990; Lambert et al. 1995; though see Erman & Iwamae 1995). Lambert ct al. (1995) 
adopted / = 0.00123±0.00016 for the A-X (2-0) band, based on comparisons with the UV D-X (0-0) and 
F-X (0-0) bands toward ( Oph and the D-X (0-0) band /-value noted above. That empirical value lies 
roughly halfway between the theoretical estimate / = 0.00144±0. 00007 of Langhoff ct al. (1990) and the 
experimental value / = 0.00100 adopted by van Dishoeck & Black (1989). Similar comparisons based on 
our fits to the profiles of the F-X, D-X, and A-X bands observed toward HD 24534, HD 27778, HD 147888, 
HD 206267, and HD 207198 suggest a somewhat higher /-value for the A-X (2-0) band, however. We have 
therefore adopted the value / = 0.00140 — an average of the / = 0.00136±0. 00015 measured by Erman & 
Iwamae (1995) and the value calculated by Langhoff et al. (1990). On average (toward HD 24534, HD 27778, 
HD 204827, HD 206267, HD 210121, and several other stars with data for these two Phillips bands), the 
strength of the C2 A-X (3-0) band then is consistent with / ~ 0.00065 — i.e., weaker than the (2-0) band 
by the expected factor <~ 2.2 (van Dishoeck & Black 1982). The C2 column densities derived from fits to 
the A-X (2-0) and (3-0) bands (with the adopted /-values) are listed in Table 9. The values for ( Oph 
given in Table 9 were derived from published equivalent widths (Hobbs & Campbell 1982; Danks & Lambert 
1983; van Dishoeck & Black 1986b; Lambert et al. 1995). The adopted C2 band /-values yield quite good 
agreement for the column densities derived from the D-X (0-0), A-X (3-0), and A-X (2-0) bands toward 
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£ Oph as well. 

For most of the sight lines in this study, observations of the UV and/or optical C 2 bands, together with 
the corresponding adopted /-values, appear to yield fairly reliable and mutually consistent column densities 
Nj for the various rotational levels. We have therefore adopted the weighted average Nj listed in Table 9 for 
sight lines with data for multiple bands. Figure 15 compares the average C2 rotational populations toward 
HD 24534 with those derived from the individual UV (upper panel) and optical (lower panel) bands. In 
general, there is good agreement; the main exceptions are the low values found for J = 8 and 10 from the 
F-X (0-0) band. [In estimating the average C 2 in levels J = 8 and J = 10 toward HD 24534, HD 27778, 
and HD 206267, we have doubled the contributions from the F-X (0-0) band, in view of that apparent (and 
unexplained) deficiency, compared to the values obtained for those levels in the other C2 bands.] 

While we have C2 data for J = 0-8 for most of the sight lines in our sample, examination of the C 2 
rotational population distributions toward HD 24534, HD 27778, and HD 207198 (for J = 0-18) and ( Oph 
(for J = 0-24; Lambert et al. 1995) — and extrapolation of those distributions to even higher J — suggests 
that a non- negligible fraction of the total C2 is in rotational levels higher than J = 8. Toward HD 24534 
and C Oph, for example, roughly 20% and 40% of the total C 2 , respectively, appear to be in rotational levels 
above J = 8. In order to estimate the total C 2 column density in all rotational levels for each line of sight, 
we have multiplied the observed total for J = 0-6 (-/Vo-6> which we have for all sight lines) by a factor 
No-20/ 'No— 6 derived from the theoretical rotational distribution that best fits the observed Nj for that line 
of sight (see below, § 4.5.2). Although the theoretical distributions are calculated only to J — 20, No-20 is 
within 10% of AT to t for all the sight lines in our sample; the factor N0-20/N0-6 ranges between about 1.1 
and 2.3 for those sight lines. The estimated total -/V(C 2 ) (listed in the last columns of Tables 7 and 9) are 
generally somewhat smaller than the corresponding values given in Thorburn et al. (2003) — due mostly 
to differences in the adopted A-X band /-values, but also to differences in the extrapolation to all J, in the 
estimation of saturation effects, and (in some cases) in the measured equivalent widths. 

3.3.3. Adopted Column Densities 

Table 10 summarizes the total column densities of H 2 , 12 CO, 13 CO, and C 2 found for the ten lines of 
sight in this study, together with the resulting 12 CO/H 2 , 12 CO/ 13 CO, and C 2 /H 2 ratios. The H 2 column 
densities are from Rachford et al. (2002), except for the value toward HD 147888 (Cartledge et al. 2004). 
The column densities for CO and C 2 are from this paper, except for the 12 CO and 13 CO values toward 
HD 24534 (Sheffer et al. 2002a). Because the STIS and FUSE data do not resolve the individual gas 
components seen in the higher resolution optical spectra, a component-by-component analysis is rendered 
uncertain (especially for H 2 ). The discussions below therefore consider only the total column densities - 
which should in any case represent the primarily molecular material in each line of sight. 

4. CO, C 2 , and Physical Conditions in Translucent Sight lines 

The so-called "translucent" clouds are thought to be intermediate or transitional objects between diffuse, 
primarily atomic clouds and dense molecular clouds (e.g., van Dishoeck & Black 1988, 1989). It has been 
difficult, however, to isolate and characterize specific examples of "true" translucent clouds using optical/UV 
absorption lines, even with the more sensitive spectroscopic capabilities provided by FUSE and HST. In a 
recent FUSE survey of interstellar H 2 in translucent sight lines, Rachford et al. (2002) obtained accurate 
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measurements of the H2 column densities in the J = 0-1 rotational levels of the ground electronic and 
vibrational state for 23 lines of sight with Ay between 0.6 and 3.4 mag. While all the sight lines exhibit 
strong H2 absorption, with 7V(H2) > 3 x 10 20 cm~ 2 , none is characterized by the very high molecular 
fractions /(H2) > 0.8 expected for "true" translucent clouds. Rachford et al. concluded that most of the 
observed sight lines contained multiple components consisting of mixtures of diffuse and denser gas, with 
"little evidence for the presence of individual translucent clouds". That assessment has been confirmed 
toward HD 185418 and HD 192639 via detailed studies of those sight lines (Sonnentrucker et al. 2002, 2003). 
Similarly, Snow, Rachford, & Figoski (2002) did not find significantly enhanced depletion of iron (as would 
have been expected for translucent clouds) , in essentially the same sample of sight lines (though it is possible 
that individual components with very severe depictions could be masked by other components with less 
severe depletions in the low-resolution FUSE spectra). Nine of the ten stars considered here were included 
in the FUSE H 2 and Fe II surveys. 



4.1. Behavior of CO 

Because of the abundance of CO and its importance as a tracer and diagnostic of molecular gas, a number 
of studies have been undertaken to understand and predict its behavior in diffuse, translucent, and denser 
clouds. Semi-analytic models (e.g., Federman & Huntress 1989) identify the most significant channels for the 
formation and destruction of each molecular species and their dependences on the local physical conditions. 
Comparisons of the predictions of such simplified chemical models with observed column densities, for an 
ensemble of sight lines, can be used to constrain (sometimes poorly known) rate coefficients and other 
parameters (e.g., the cosmic ray ionization rate) of the models. More detailed numerical models (e.g., van 
Dishoeck & Black 1988, 1989; Viala et al. 1988; Warin, Bcnayoun, & Viala 1996; Le Petit et al. 2006) 
typically include more extensive reaction networks and radiative transfer effects, and attempt to determine 
the abundances of various species (including the individual rotational levels for H2 and CO) as functions of 
depth within the model clouds. All those detailed models, however, are for single, isolated plane-parallel 
clouds (often of uniform density). 

In diffuse and translucent clouds with Ay < 3 mag, the primary route for the production of CO begins 
with the reaction C + + OH — > CO + + H. The CO + then reacts with H2 to form HCO + , which subsequently 
forms CO via dissociative recombination. In thicker, denser clouds, reactions involving O, CH, C2, and 
other species may also contribute. The destruction of CO is via line absorption into predissociated bound 
states; the relevant transitions are at far-UV wavelengths between 912 and 1100 A (e.g., van Dishoeck & 
Black 1988; Viala et al. 1988; Warin et al. 1996). At high enough column densities [V( 12 CO) > 10 14 
cm -2 ], those transitions become saturated, self-shielding begins to reduce the photodissociation rate, and 
the column density of 12 CO increases rapidly. Recent observations of some of the far-UV CO bands suggest 
that the self-shielding by those bands may be more effective than previously thought (Sheffer, Federman, 
& Andersson 2003). In addition, some of the 12 CO transitions coincide with strong absorption lines of H I 
and/or H2, which can dominate the shielding of 12 CO at lower column densities. The other isotopomers 
(e.g., 13 CO, C 18 0) are much less abundant, and so do not self shield (nor are they effectively shielded by 
cither H 2 or 12 CO). At much higher column densities, reactions between CO and other species become the 
dominant destruction mechanism for CO. 

In the sections below, the observed column densities of CO, H2, and several other species will be 
compared with predictions from both the models of van Dishoeck & Black (1988, 1989), Warin et al. (1996), 
and Le Petit et al. (2006). The T, H, and I model sequences of van Dishoeck & Black cover the range 
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of extinctions and H 2 column densities relevant to most current observations of CO absorption (Ay < 2.5 
mag) . The T models are characterized by the average interstellar radiation field, total hydrogen densities % 
= n(H I) + 2n(H 2 ) of 300-1000 cm" 3 , temperatures of 15-60 K, H 2 column densities of 0.5-4.0 x 10 21 cm" 2 , 
and visual extinctions of 0.7-5.2 mag. Several of the models allow for some variation in the temperature 
and density with depth. The H models have half the average interstellar field, nn = 500 cm" 3 , T = 40 K, 
A(H 2 ) = 0.3-3.0 x 10 21 cm' 2 , and Ay = 0.4-3.8 mag. The I models have ten times the average interstellar 
field, n H = 2000 cm' 3 , T = 30 K, iV(H 2 ) = 0.5-10.0 x 10 21 cm' 2 , and Ay = 0.8 12.8 mag. All three series 
assume depletions of carbon, nitrogen, and oxygen of 0.1-0.4, 0.5, and 0.6, respectively (with respect to 
the solar or "cosmic" abundances adopted at that time); all three use the average Galactic extinction curve 
[though van Dishoeck & Black (1989) also considered the effects of curves that are shallower (HD 147889) 
and steeper (HD 204827) in the UV]. 

Warin et al. (1996) computed models for diffuse clouds [Ay = 1.09 mag, A(H) = 2 x 10 21 cm" 2 , nn 
= 500 cm" 3 ], translucent clouds [Ay = 4.34 mag, iV(H) = 8.1 x 10 21 cm" 2 , n H = 1000 cm" 3 ], and dark 
clouds [Ay = 10.86 mag, iV(H) = 20 x 10 21 cm" 2 , n H = 10 4 cm" 3 ] — both at constant T and in thermal 
equilibrium. All the Warin et al. models are for uniform density slabs; all assume depletions of 0.1 for both 
carbon and oxygen. Unlike the van Dishoeck & Black models, the Warin et al. models allow for differences 
in rotational excitation temperature for different J — resulting from the coupling between level populations, 
self-shielding, and photodissociation rates. The total unshielded photodissociation rate adopted by Warin 
et al. (1996), based on more recent CO spectroscopic data, is of order 20% smaller than that used by van 
Dishoeck & Black (1988). 

Le Petit et al. (2006) give an overview of the current Meudon photon-dominated region (PDR) code (e.g., 
Le Bourlot et al. 1993), which computes the radiative transfer, thermal balance, and chemistry (including 
the excitation of H 2 and CO) for a plane-parallel cloud in steady state. In principle, both temperature and 
density can vary with depth in the cloud. While Le Petit et al. do not give extensive numerical results, 
they do provide a coarsely sampled table of column densities of H 2 , CO, and several other species (for nn = 
100 cm -3 , the standard interstellar radiation field, and total Ay ranging from 0.2 to 7.0) and several figures 
showing the column densities of CH, CN, and CO, as functions of nn/x, for a set of diffuse cloud models 
[Ay = 1) 1h ranging from 20 to 2000 cm -3 , and \ (the scaling factor for the radiation field) from 0.1 to 10]. 
All the models use the current observed gas-phase interstellar abundances of carbon, nitrogen, and oxygen. 

The various models of cloud chemistry suggest that the column density of CO generally will increase 
with increasing local density, H 2 column density, and optical depth. For diffuse clouds, Federman et al. 
(1980) concluded that the local density of CO would be proportional to [n(H 2 )] 2 (from formation) and to 
exp[r(OH)+r(CO)] (from destruction; the optical depths r are for the wavelength ranges appropriate for 
photodissociation of OH and CO), but that the density dependence would be only linear with n(H 2 ) once 
reactions with other species dominate the destruction of CO. Isobaric models with mean densities of 150 and 
2500 cm" 3 appeared able to account for most of the range in N(CO) observed at any given A^(H 2 ). (The 
role of line absorption and self-shielding in the dissociation of CO was not yet appreciated, however.) For 
clouds with Ay < 2.5, the more detailed models of van Dishoeck & Black (1988), which include self-shielding, 
predict that: (1) increasing the gas-phase carbon abundance by a factor of four [for constant density and 
A(H 2 )] increases the total CO column density by a factor of 1.5-3.0 (models T1-T4); (2) decreasing the 
radiation field by a factor of 2 [for constant depletion, density, and A(H 2 )] increases iV(CO) by a factor of 
3-4 (models H2-H4 vs. models T1-T3); and (3) increasing A(H 2 ) by a factor of 2 (for constant depletion and 
radiation field) increases N(CO) by more than an order of magnitude (models T1-T4, H2-H5, and 11-14). 
The Warin et al. (1996) translucent cloud model has twice the density, four times the H 2 column density, 
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and roughly 3000 times the CO column density of the diffuse cloud model. The earlier models suggest that 
CO becomes the dominant carbon-containing species for clouds with Ay > 3, with the transition between 
C + (and C°) and CO typically occurring at optical depths ry between 0.5 and 1.0 (van Dishocck & Black 
1988; Warin et al. 1996) — where the hydrogen is already almost completely in molecular form. The very 
high predicted molecular fractions /(H 2 ) > 0.9 — even for the total (integrated) column densities in the 
diffuse cloud models — have not yet been observed, however (Rachford et al 2002). In the Meudon PDR 
models with nn =100 cm~ 3 and x = 1, C + remains the dominant form of carbon even for Ay = 7; in the 
diffuse cloud models, iV(CO) does not show a very strong dependence on density, for rin/x between 10 and 
1000 cm~ 3 . Additional aspects of the behavior of CO (e.g., rotational excitation, isotopic abundances), as 
predicted by the various models, will be discussed further below (§ 4.5.1 and § 4.6, respectively). 

4.2. Behavior of C 2 

While C 2 is not as abundant as CO, it does play a role in the networks of chemical reactions occurring 
in diffuse and translucent clouds, and analyses of its absorption lines can provide constraints on the density, 
temperature, and radiation field in those clouds. Federman & Huntress (1989) have described the key 
reactions controlling the abundance of C2; van Dishoeck & Black (1982) have provided a detailed discussion 
of its excitation in diffuse and translucent clouds. A number of papers (e.g., van Dishoeck & de Zeeuw 1984; 
Federman & Lambert 1988; van Dishoeck & Black 1989; Federman et al. 1994a) have discussed observations 
of C2 absorption in the context of models for its abundance and/or excitation. 

In diffuse and translucent clouds, the formation of C2 begins with the reaction C + + CH — > C^ + H; 
subsequent hydrogen abstraction reactions and dissociative recombinations yield C2 via several channels. 
Destruction of C2 can occur by photodissociation (more important at relatively low optical depths) or by 
reactions with atomic oxygen (yielding CO) or nitrogen (yielding CN), which become dominant at higher r. 
The detailed models T1-T6 of van Dishoeck & Black (1989), with total Ay = 0.7-5.2 mag, predict a roughly 
constant C2/H2 ratio ~ 3.3-5.7 x 10~ 8 over that range. The C2 abundance is somewhat enhanced if the 
oxygen depletion is more severe than the value of 0.5 usually assumed in the models, is essentially unaffected 
by the assumed nitrogen depletion, and is somewhat affected by the shape of the far-UV extinction. As 
a "second generation" molecule dependent on the prior existence of H2 and CH, C2 is likely to be more 
concentrated in somewhat colder, denser regions. The models of van Dishoeck et al. (1991) suggest that 
the average density traced by C2 may be either higher or lower than that traced by CO, depending on the 
structure of the cloud and how it is illuminated by the ambient radiation field. The rotational excitation of 
C2 will be discussed below (§ 4.5.2). 

4.3. Column Density Comparisons 

In order to compare the column densities of 12 CO and 13 CO obtained in this paper with those found 
for other Galactic sight lines, we have collected CO column densities derived from UV absorption- line data 
from the literature. Table 16 (Appendix C) lists total column densities for H2, CH, C 2 , CN, and CO for 74 
sight lines for which CO has been measured (53 detections, 21 upper limits) — a reasonably large sample 
which includes recently improved iV(CO) values for both diffuse and more heavily reddened sight lines. For 
a number of those sight lines, column densities also have been recently reported for HD (Lacour et al. 2004) 
and for C 3 (Maier et al. 2001; Roucff et al. 2002; Galazutdinov et al. 2002; Oka et al. 2003; Adamkovics, 
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Blake, & McCall 2003). Tabic 17 (Appendix) lists various derived quantities (e.g., excitation temperatures, 
column density ratios) for the sight lines with data for CO, H 2 , CN, and CH. There are now 19 sight lines 
with relatively well determined N(CO) > 10 15 cm" 2 derived from UV absorption lines [plus a few others 
with values from mm- wave CO absorption (Liszt & Lucas 1998, 2000)], compared to one such sight line (( 
Oph) in the early survey of Federman et al. (1980). This enlarged higher column density sample allows 
more meaningful comparisons with the predictions of theoretical models of translucent clouds. Comparisons 
among these molecular species (and of various ratios of these species) may provide some clues as to their 
behavior in different environments and to the physical conditions characterizing the clouds in which they 
reside (e.g., Federman et al. 1994a). While there are likely to be some differences in the spatial distribution 
of these species, they all trace the main components in each sight line, where the bulk of the molecular 
material is located. 



4.3.1. N(CO) andN(C 2 ) vs. N(H 2 ) 

Figure 16 shows plots of the column densities of 12 CO and 13 CO versus A(H 2 ). In the plots, the different 
symbols refer to the sources of the CO data (Copernicus, IUE, FUSE, or HST); the filled points represent 
sight lines for which 13 CO has been detected. The dotted lines in both figures indicate the relationships 
predicted by the models of van Dishoeck & Black (1988); series H, T, and I (increasing radiation field; all 
with carbon depletion 5c = 0.4) are left-to-right in each case. The letters "W" indicate the values predicted 
by the diffuse, translucent, and dense cloud models of Warin et al. (1996). The dashed lines show the 
predictions of the Mcudon PDR models with nn = 100 cm -3 , \ = 1) an d ranging from 0.2 to 7 (Le Petit 
et al. 2006). The observed column densities confirm the steep increase of A(CO) with A(H 2 ), for A(H 2 ) 
> 3 x 10 20 cm~ 2 , as predicted by the various models and as noted by Federman et al. (1980). The much 
more extensive data for A( 12 CO) > 10 15 cvnT 2 , however, appear to indicate a somewhat steeper increase 
for the higher column density sight lines than that predicted by Federman et al. — presumably due to the 
effects of self-shielding. The models of van Dishoeck & Black (1988) appear to reproduce the trends in the 
observed data reasonably well, but predict too little 12 CO by a factor of order 2-3 [for a given A(H 2 )]. For 
the sight lines with detected 13 CO, the observed A( 13 CO) appear to agree somewhat better with the model 
predictions than do the observed A( 12 CO). The models of Warin et al. (1996) and Le Petit et al. (2006) 
predict much smaller CO column densities than the observed values. While the density in these models 
from Le Petit et al. is relatively low, the Meudon diffuse cloud models (Ay = 1) do not show an increase in 
N(CO) until 71h is greater than about 1500 cm~ 3 (for \ = 1). More extensive exploration of the parameter 
space available in the Meudon code (though beyond the scope of this paper) might prove instructive. 

As suggested by the three series of models (H, T, I), some of the scatter in the relationship between 
A(CO) and A(H2) may be due to differences in the strength of the ambient radiation field. For sight lines 
with column density data for the lowest few rotational levels of H 2 , the relative amount in J = 4 (which 
is assumed to be populated primarily via photon pumping) may be used to estimate the strength of the 
radiation field (Jura 1975; Lee et al. 2002; Welty et al. 2006). The four open squares just to the left of the 
high radiation field (I) model sequence correspond to four sight lines in Cep OB3 (Pan et al. 2005). For 
three of those sight lines, the H 2 rotational populations determined by Pan et al. do suggest fairly strong 
ambient fields. The molecular material seen toward HD 37903 [open triangle near A(H 2 ) = 8 x 10 20 cm~ 2 
and A(CO) = 10 14 cm -2 )] also is likely subject to an enhanced field (Lee et al. 2002). However, there 
are other sight lines with similarly high inferred fields whose CO column densities lie closer to the model 
sequences characterized by average or low fields. In addition, Lacour et al. (2005b) obtained progressively 
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higher b- values for the higher rotational levels of H2 in several moderately reddened lines of sight (including 
HD 73882, HD 192639, and HD 206267) and argued that the column densities of those higher J levels could 
be significantly overestimated under the usual assumption of constant b for all J. If A^(H2) is too high, then 
the radiation field will be similarly overestimated. While differences in the radiation field probably do affect 
the abundances of CO observed in different sight lines, other factors must also play a role. 

Increases in the gas-phase carbon and oxygen abundances and/or in the self-shielding of 12 CO (Sheffer et 
al. 2003) in the models could yield somewhat better agreement between the observed 12 CO column densities 
and those predicted by the models. While interstellar C II column density data are still rather limited, there 
does appear to be a fairly well-defined average value C/H = 161±17 ppm for both diffuse and translucent 
sight lines (Sofia et al. 2004). For a total (gas + dust) carbon abundance of 245 ppm (Lodders 2003), the 
carbon depiction 8c is thus about 0.66 [or a logarithmic depiction D(C) ~ —0.18 dcx]. Only the sight line 
toward HD 27778 shows evidence for a gas-phase carbon abundance significantly below that average value 
[< 108 ppm (3-er)] 2 ; the carbon abundance toward X Per (155±45 ppm), with higher Ay and /(H 2 ) and 
similar molecular abundances, is quite consistent with the average (Sofia, Fitzpatrick, & Meyer 1998). The 
contribution of CO to the overall gas-phase carbon abundance in those two sight lines is only ~ 6 ppm. The 
Warin et al. (1996) models assumed a total carbon abundance of 360 ppm and 8c = 0.1. The resulting 
gas-phase carbon abundance (~ 36 ppm) is a factor of ~ 4.5 lower than the observed average interstellar 
value. Increasing the gas-phase carbon abundance would move the predicted iV(CO) closer to the observed 
values, but would generally not be enough to yield good agreement. The van Dishoeck & Black (1988, 
1989) models assumed a total carbon abundance of 467 ppm and <5c = 0.1-0.4. For <5c = 0.4, the resulting 
gas-phase carbon abundance (~ 190 ppm) is slightly higher than the average interstellar value, so that other 
adjustments (e.g., increasing the 12 CO self-shielding) would be required to obtain better agreement with the 
observed N(CO). 

Several other suggested modifications to the (relatively) simple early models may be able to increase 
the CO abundances predicted for diffuse and translucent clouds. (1) Both van Dishoeck & Black (1989) and 
Warin et al. (1996) noted that one way would be to reduce the (relatively poorly known) intensity of the 
radiation field below 1200 A. The field would need to be weaker by an order of magnitude (or more), however, 
for the Warin et al. diffuse cloud model to reproduce the observed abundances of 12 CO and 13 CO (even with 
an increased carbon abundance). (2) Liszt & Lucas (1994) concluded that the low predicted values of N (CO) 
toward £ Oph should be ascribed to insufficient production of CO (rather than excessive photodissociation), in 
view of the similar underprediction of its immediate precursor HCO + . An increase in the rate of the initiating 
reaction C+ + OH -> CO+ + H (e.g., Dubcrnct, Gargaud, & McCarroll 1992; Federman et al. 2003) would 
yield corresponding increases in both HCO + and CO. (3) In view of observational evidence for complex, 
inhomogeneous structure in molecular clouds, Spaans (1996) explored the consequences of having higher 
density clumps embedded in a lower density interclump medium. On the one hand, dissociating radiation 
may penetrate further into the interior of such a clumpy structure than is the case for the previously modeled 
homogeneous slabs. Both theoretical considerations (§ 4.1) and empirical evidence (§ 4.4), however, suggest 
that CO will be enhanced at higher densities. Spaans found that CO (which is located mostly in the denser 
clumps) could actually be enhanced by factors of 5-6 in the clumpy models, compared to the abundances 
found for the corresponding homogeneous models. (4) Because the photodissociation of CO is dominated by 
line absorption, it can depend on details of the velocity distribution in the molecular gas. Consideration of the 
effects of a turbulent velocity field has suggested that resulting changes to the CO absorption-line profiles 



2 Whilc the CO column density toward HD 27778 estimated by Federman ct al. (1994a) seemed suggestive of an unusually 
high CO abundance in that line of sight (Sofia ct al. 2004), our value for 7V(CO) is a factor of 2 smaller. 
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could significantly increase the self-shielding of CO — thus increasing the CO abundance in moderately 
dense molecular clouds (e.g., Rollig, Hegmann, & Kegel 2002). Much smaller effects are seen, however, 
when the full CO rotational structure is incorporated (instead of the single-line approximation used in the 
initial calculations). Finally, however, if many sight lines contain a significant amount of H2 without much 
associated CO (e.g., in more diffuse molecular gas), then the agreement with the models would not be as 
good as Figure 16 would otherwise seem to suggest. Detailed studies of individual lines of sight (to separate 
the contributions from different molecular components) and comparisons with more sophisticated models 
(incorporating these additional physical effects) may yield better agreement between observed and predicted 
CO abundances. 

Figure 17 shows the relationship between the column densities of C2 and H2. The solid lines are 
linear least-squares fits (weighted and unweighted) to the points representing detections of both quantities, 
taking into account that both quantities have associated uncertainties. 3 Table 11 lists the linear correlation 
coefficients (r), the slopes of the best-fit lines, and the residual scatter about the best-fit lines found for this 
comparison (and for others discussed below). The dotted line, which shows the roughly linear relationship 
predicted by models T1-T6 of van Dishoeck & Black (1988, 1989), falls slightly below many of the observed 
data points. Much of the systematic difference between the models and the observed data may be ascribed 
to the depletion of carbon (5c = 0.1-0.15) assumed in the models — as the corresponding gas-phase carbon 
abundances are factors of 2-3 lower than the values observed for nearly all diffuse and translucent clouds. 
Apart from that systematic offset, the observed data seem to follow a somewhat steeper relationship (slope ~ 
1.8) than the roughly linear predicted trend (though the ranges in both column densities are rather limited); 
the much steeper slope found for the unweighted fit is due largely to the point for 40 Per. While Federman 
& Huntress (1989) suggested progressively more severe depletion of carbon (to a factor of 15 by Ay = 4 
mag.) was required to explain the C2 abundances in translucent sight lines, none of the sight lines in Fig. 17 
has Ay higher than 2.7 mag. 



4.3.2. N(CO) vs. Column Densities of CN, CH, C 2 , C 3 , HD, and KI 

Figures 18 and 19 show plots of the column density of 12 CO versus the column densities of CN, CH, C2, 
and K I; similar comparisons were also made with the column densities of C3 and HD. Again, the different 
symbols denote the sources of the CO data, with the filled points representing sight lines for which 13 CO 
has been detected, and the solid lines are the linear least-squares fits to the data. For the full sample, there 
are clear, significant correlations (r from 0.733 to 0.905) with CN, CH, HD, and K I (Table 11). The smaller 
correlation coefficients for C2 and C3 appear to be due to several sight lines (40 Per, HD 200775) with 
very high CO/C2 ratios and to the small range in column density for the few sight lines with detections 
of C3, respectively. The strongest correlations are with CN (r = 0.900; slope ~ 1.50) and with HD (r = 
0.905; slope ~ 1.75); the relationship with CN shows significantly smaller scatter, however. The slopes of 
the relationships with CH, C2, and K I are much steeper (between about 3.1 and 4.0), though consideration 
of the upper limits for some sight lines suggests that the relationships may be somewhat shallower at lower 
column densities. In view of the scatter in those relationships, it would be difficult to accurately predict 
7V(CO) from the column densities of CH, C2, or K I for any individual sight line. 

The dotted lines in Figures 18 and 19 give the relationships between the column density of 12 CO and 



3 A slightly modified version of the subroutine regrwt.f, obtained from the statistical software archive maintained at Penn 
State (http://www.astro.psu.edu/statcodes), was used to perform the fits; see Welty & Hobbs (2001). 
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those of CN, CH, and C 2 predicted by models T1-T6 of van Dishoeck & Black (1988, 1989). In all three 
cases, the average slopes of the predicted relationships are very similar to those determined from the fits 
to the observed data, but the predicted N(CO) are lower than the observed (fitted) values by factors of 
2-4. To first order, the low values assumed for the carbon abundance in the models presumably should not 
significantly affect the predicted ratios of CO to CN or CH. The gas-phase abundances assumed in the van 
Dishoeck & Black (1988, 1989) models for oxygen and nitrogen are closer to the current observed values 
(Meyer, Cardelli, & Sofia 1997; Cartledge et al. 2004). The dashed lines in Figure 18 show the predictions 
of the Meudon diffuse cloud models (Le Petit et al. 2006), for x — 1 and nn ranging from 20 to 2000 cm~ 3 . 
The predicted N(CO) would be consistent with the observed trends for densities of about 200-500 cm -3 - 
though there are sight lines with Ay < 1 and iV(CO) much higher than those predicted values. 

As was found from fitting the observed CO absorption-line profiles, these column density comparisons 
suggest that the distribution of CO is most similar to that of CN. The correlation coefficient with CN is 
larger than those for CH, C2, C3, and K I, and the slope of the relationship with CN is closer to unity. 
Weak components in K I and CH, for example, are probably much weaker in CO, though they may still be 
detectable in the strongest CO bands (e.g., Sheffer et al. 2003; see also the profiles for HD 206267 in Fig. 3). 
While the abundances of both CO and CN appear to increase strongly with overall column density, however, 
the reasons for the increase may be different for the two species. Because the production of CN depends on 
the prior existence of C2, CH, and/or NH, the "third-generation" molecule CN is expected to trace colder, 
denser, more shielded regions where those precursor molecules have already become abundant (Federman, 
Danks, & Lambert 1984; Federman et al. 1994). On the other hand, the strong increase in N( 12 CO) is likely 
to be due both to increased formation at higher densities and to reduced photodissociation resulting from 
the onset of self-shielding, once sufficient 12 CO is present [i.e., iV( 12 CO) > 10 14 cm -2 ]. 



4.4. Variation of the CO/H 2 Ratio 

Both theoretical modeling and observations of the CO abundance indicate that the CO/H2 column 
density ratio can vary by three orders of magnitude, from < 10~ 7 in diffuse clouds to about 10~ 4 in dense 
molecular clouds (van Dishoeck & Black 1988; van Dishoeck & Black 1990; Warin et al. 1996); the transition 
between those values occurs in the translucent regime. The analogous ratio of N(H2) to CO emission 
strength [the so-called "conversion factor" X C o = ^(H 2 )/I( 12 CO), of order 2 x 10 20 cm~ 2 (K km s -1 ) -1 
in our Galaxy] is often used to derive the H2 column density in dense molecular clouds, since H2 is often 
difficult to measure directly in such regions (Lacy et al. 1994). That conversion factor reflects the relationship 
between the width of the (saturated) 12 CO J = 1-0 emission line and the mass of a gravitationally bound 
dense cloud. A unique conversion factor cannot be defined, however, in diffuse and translucent clouds (e.g., 
Magnani & Oncllo 1995) — which are not gravitationally bound and where iV(CO) is not high enough for 
the emission lines to be saturated. 

Figure 20 shows the same data as in Figure 16, but plotted as A^(CO)/A^(H 2 ) versus 7V(H 2 ). Again, 
the dotted lines show the relationships predicted by models H, T, and I of van Dishoeck & Black (1988), 
the letters "W" show the values predicted by the models of Warin et al. (1996), and the dashed line (for 
12 CO) shows the predictions of the Meudon code (Le Petit et al. 2006). In our sample often sight lines, the 
12 CO/H 2 ratio varies by a factor of about 100 — from 2.8 x 10~ 7 toward HD 192639 to 2.8 x 10~ 5 toward 
HD 73882 (Table 10); values as low as 2.5 x 10~ 8 are seen in the larger sample. For our sight lines, the 
variations in the overall column densities of both H2 and total hydrogen (H I + 2H2) are much smaller; for 
example, seven of the ten have N(B. 2 ) in the range 4.9-7.2 x 10 20 cm~ 2 (Tabic 10). The differences in the 
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CO/H2 ratios in our sight lines thus are due primarily to variations in the CO column densities, which in 
turn depend on the physical conditions in the molecular gas in those particular sight lines. Since the total 
hydrogen content is similar for those sight lines, variations in the shielding of the CO bands by H and H2, 
lines cannot explain the differences in the CO column density. It seems most likely that the range in iV(CO) 
observed for our ten sight lines is due to differences in the CO photodissociation rate (ambient radiation 
field, self-shielding) and/or in the CO formation rate (local hydrogen density). 

In all ten of our sight lines, CO contains only a minor fraction of the total carbon abundance. For the 
current solar carbon abundance C/H ~ 2.5 x 10~ 4 (e.g., Asplund, Grevesse, & Sauval 2005) and typical 
interstellar depletions (5c ~ 0.66), the total gas-phase ratio C/H2 would be about 3.2 x 10~ 4 (if hydrogen 
were fully molecular). The highest values of the CO/H2 ratio in our sample — ~ 2-3 x 10~ 5 toward HD 
27778 and HD 73882 — are more than an order of magnitude smaller. Most of the gas-phase carbon is still 
in the form of C + (Sofia et al. 2004; see also Liszt & Lucas 1998). 

In order to explore observationally the relationship between the CO/H2 ratios and local physical con- 
ditions (e.g., T, njj, radiation field), empirical diagnostics of those parameters are needed. The simplified 
chemical model of Federman & Huntress (1989) implies, for example, that the CN/CH ratio in diffuse molec- 
ular gas should be roughly proportional to the square of the local hydrogen density (n 2 H ) (Cardelli, Federman, 
& Smith 1991) — suggesting that the CN/CH column density ratio may be used as a density indicator. The 
more detailed models of van Dishoeck & Black (1986a, 1989), however, predict a somewhat weaker increase 
of JV(CN)/AT(CH) with nu- While the densities estimated by Pan et al. (2005) for individual components 
toward stars in Ccp OB2 suggest that the CN/CH ratio may be roughly proportional to n][ 2 (with much 
scatter), the densities inferred below from C2 (§ 4.5.2) appear to indicate a somewhat steeper relationship. 
The CN/CH ratio thus does appear to be correlated with density, but other factors (e.g., Ay and the steep- 
ness of the extinction curve in the far-UV) may contribute to the scatter (see also Cardelli 1988; Welty & 
Fowler 1992). Some CH may also be formed together with CH+ in lower density gas (Federman, Welty, & 
Cardelli 1997a; Zsargo & Federman 2003), but such contributions generally appear to be fairly minor for 
sight lines in which CN has been detected. 

Figure 21 shows the observed ratios A( 12 CO)/A(H 2 ) and A( 13 CO)/A(H 2 ) versus A(CN)/A(CH) for 
sight lines with data for those species (Appendix C; Table 17). There are clear correlations in both cases 
(r = 0.815 and 0.868); the solid lines show the linear least-squares fits to the data (with slopes of order 
1.55-1.75; Table 11). The dotted line for 12 CO/H2, showing the relationship predicted by models T1-T6 of 
van Dishoeck & Black (1989) [for 5c — 0.15, 5^ — 0.6, So — 0.5, and grain model 2], is very close to the best 
fit to the observed data. The letter "M" , corresponding to the values predicted by the Meudon code for Ay 
= 1, x = 1, and 7Jh = 100 cm -3 , also lies very close to the observed trend. These correlations are perhaps 
not surprising — in view of the good correlation between N( 12 CO) and A(CN) and the known roughly 
linear relationship between -/V(H 2 ) and iV(CH) (e.g., Danks, Federman, & Lambert 1984). Nonetheless, 
the correlations do suggest that differences in the local density «h [as indicated by N (CN) / N (CK)} may 
play a significant role in explaining the variations in the CO/H 2 ratios observed in diffuse and translucent 
molecular gas. As the dispersions relative to the best fits to the data (~ 0.2 dex) are larger than the typical 
measurement uncertainties, it is likely that other factors (e.g., radiation field, UV extinction, cloud structure; 
see § 4.3.1) also affect the CO abundances, however. In contrast, the ratio A'(C 2 )/A^(H 2 ) does not show any 
significant dependence on N (CN) / N (CU) . 

The column density ratios observed for our ten sight lines are quite consistent with the trends observed 
for the larger sample. The largest ratios of both CO/H 2 and CN/CH are observed toward HD 73882 and 
HD 27778 (Table 17); the smallest ratios are found toward HD 185418 and HD 192639 (where CN is weak). 
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For the three targets in the Cep OB2 association, the CN/CH ratio toward HD 206267 is higher than those 
toward HD 207198 and HD 210839 by factors of 1.4-2.1, while the CO/H 2 ratio is about 2.6-3.0 times larger. 
Simple chemical modeling of the CH and CN absorption by Pan et al. (2005) yields densities of 525-1150 
cnr 3 toward HD 206267, 225-425 cm" 3 toward HD 210839, and 150 cm~ 3 toward HD 207198 for the main 
components seen in CN. Analysis of the C2 rotational excitation (see below) implies slightly higher densities 
toward HD 206267 than toward HD 207198. For those three Cep OB2 sight lines, the CN/CH and CO/H 2 
ratios thus both appear to increase with density. 

4.5. Rotational Excitation of CO and C 2 

Because the populations of the ground and excited rotational levels of molecules such as H 2 , CO, and 
C 2 can depend on both collisional and radiative processes, comparison of the observed rotational excitation 
of those molecules with predictions from theoretical models can provide information on the density, tem- 
perature, and radiation environment characterizing the molecular gas. In general, measurements of both 
lower and higher rotational states are needed to disentangle the collisional and radiative effects. The rel- 
ative rotational level populations are often expressed in terms of rotational excitation temperatures = 
[(Ej/k) — (Ei/k)] I [\n(Ni/gi) - \n(Nj/gj)] (usually for adjacent levels i < j), where Ei, Ni, and gi are the 
energy, column density, and statistical weight of level i. Because H 2 , CO, and C 2 each respond somewhat 
differently to the local physical conditions and may well trace somewhat different volumes of gas, deter- 
minations of their respective rotational excitations can provide both complementary constraints on those 
conditions and some information on the structure of the clouds. Tables 12 and 13 list various excitation tem- 
peratures derived from the rotational level populations of CO, H 2 , and C 2 , which will be used to constrain 
the physical conditions in the molecular gas. 

4-5.1. CO Excitation 

Because of the abundance and ubiquity of CO in interstellar space, observations of its emission lines 
have often been used to estimate the physical conditions in dense molecular gas (tih > 10 4 cm~ 3 , T ~ 10 
K). Although CO is a fairly good density and temperature diagnostic in such dense clouds, where collisional 
processes dominate, the use of CO as a diagnostic of the physical conditions is more problematic in diffuse and 
translucent clouds still permeated by UV photons (Wannier et al. 1997). While the detailed models of van 
Dishocck & Black (1988) incorporated both chemical reactions and various radiative and collisional processes 
in order to predict both the total CO abundance and its rotational excitation, the excitation temperature 
was assumed to be the same for all rotational levels. Warin et al. (1996) removed that constraint — and 
were thus able to study the variations in the excitation of the various CO rotational levels for "typical" 
diffuse, translucent, and dense clouds. 

For diffuse clouds, the Warin et al. (1996) models predict that the photodissociation rate remains both 
large throughout the cloud and globally independent of the J level. While CO in the excited rotational levels 
is easily photodissociated by UV photons, the relatively low gas densities do not allow efficient collisional 
population of those excited levels (achieved primarily via collisions with H 2 but also with H). The rotational 
population distribution of 12 CO (and its isotopes) is therefore sub-thermal — i.e., the excitation temperatures 
Tij of the various excited rotational levels (< 10 K for J = 1-4) are well below the assumed temperature of 
the gas (50 K) — even though the energies of those levels are less than or comparable to the thermal energy 
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(kT) at that temperature (Table 6). The models predict that the relative population of the J = 1 level is ~ 
10-50% higher than that of the J = level and that the populations of the higher levels decrease by nearly 
a factor 10 for each level above J = 1. For the typical column densities measured in diffuse clouds (and the 
sensitivity of the instrumentation available to this point), the J = 0, 1, and (perhaps) 2 rotational levels of 
12 CO should be detectable. Corresponding lines from the much less abundant isotopomers are expected to 
be weak in diffuse clouds. 

For translucent clouds, the Warin et al. (1996) models predict a strong dependence of the photodis- 
sociation rate on the rotational level, as self-shielding becomes effective for the lowest-lying levels of 12 CO. 
Because of this differential photodissociation effect, the excitation temperatures are predicted to vary both 
with rotational level and with depth in the cloud. Thermalization is still not achieved, however, and all exci- 
tation temperatures remain sub-thermal throughout the cloud. For example, the isothermal model predicts 
that the relative populations in J = 0/1/2/3 at optical depth (due to dust) tv = 0.5 are ~ 1.0/2.0/0.5/0.07, 
while the relative populations at the cloud core (tv = 2.0) are ~ 1.0/2.3/2.3/1.5. The corresponding pre- 
dicted excitation temperatures for J = 1-3 are 6-12 K at tv = 0.5 and 19-21 K at the cloud core 
(compared to the assumed constant temperature T = 25 K). The J = and 1 levels of the isotopomers 
13 CO and C 18 also should be detectable, with the 13 CO excitation temperature T i ~9K throughout the 
cloud. 

Figure 22 shows the 12 CO normalized excitation diagrams (ln[Nj/((2J+l) No)] vs. Ej/k) for the eight 
sight lines whose rotational level column densities were derived from our fits to high-resolution UV spectra 
obtained with GHRS or STIS (Table 6). The dotted lines show the extrapolated straight-line fit to the points 
for J = and 1, the slope of which is the negative of the inverse of the excitation temperature Toi. Table 12 
lists the corresponding 12 CO (and 13 CO) excitation temperatures for those eight sight lines, together with 
the CO excitation temperatures reported in the literature for eight other sight lines and the average values 
derived from the CO rotational populations in the diffuse and translucent cloud models listed in Table 2 of 
Warin ct al. (1996). For 12 CO, the previously reported Toi are all between 2.7 and 3.6 K. The Toi derived 
from our profile fits range from ~ 2.6 K toward HD 192639 and 3.3 K toward HD 185418 to ~ 11 K toward 
HD 147888, with the values for the other five sight lines all between 3.8 and 6.3 K; similar values (generally 
between 3 and 8 K) have been derived from mm-wave CO absorption observations toward extragalactic radio 
continuum sources (Liszt & Lucas 1998). The 12 CO Toi values for all 16 observed sight lines are significantly 
lower than those predicted by the Warin et al. (1996) translucent cloud models — which is perhaps not 
surprising, since none of the sight lines has Ay or iV(H 2 ) as high as the values assumed for the models. A 
number of the values determined in this paper are similar to those predicted for the thermal equilibrium 
diffuse cloud models. Even for those models, however, the predicted dependence of excitation temperature 
on rotational level (Ti 2 < T i < T23) is generally not seen for the eleven sight lines with data for J > 1, many 
of which have T i < T 12 < T23. For most of the sight lines in our sample, those three observed excitation 
temperatures are equal, within the uncertainties. 

Wannicr ct al. (1997) pointed out that CO line emission from dense molecular gas adjacent or near to 
that seen in absorption could in principle increase the populations of levels J > above the values due to 
collisional and radiative processes in the observed gas alone. Analyses of the CO rotational excitation which 
neglect that emission would thus overestimate the density in the gas seen in absorption. For most of the 
sight lines in the upper part of Table 12, Wannier et al. concluded that the proposed radiative excitation 
could account entirely for the observed rotational excitation — which thus could not be used to infer «h in 
the gas. One consequence (and test) of the proposed radiative excitation is that the excitation temperatures 
for 13 CO should be lower than the corresponding values for 12 CO, as the 13 CO emission lines should be 
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much weaker. That does not appear to be the case, however, for the seven sight lines in Table 12 with values 
for Xbi( 13 CO). In six of the sight lines, 7oi( 13 CO) is nominally greater than T i( 12 CO); in all cases, the 
two temperatures are equal, within the mutual uncertainties. And while two of the six components detected 
in both 12 CO and 13 CO absorption by Liszt & Lucas (1998) have T i( 12 CO) > T i( 13 CO), the other four 
components have T i( 12 CO) ~ T i( 13 CO). The radiative excitation mechanism proposed by Wannier et al. 
thus may not be significant, in many cases. 

It is not clear whether the higher values of Ty (~ 8-11 K) for 12 CO observed toward HD 147888 can 
be associated with a higher density in the molecular gas in that sight line. While a relatively high density 
could be suggested by the observed relatively strong absorption from the excited fine-structure states of C I 
and O I, the density inferred from the C2 rotational excitation (§ 4.5.2) is similar to the values found for 
many other sight lines. Unfortunately, 13 CO was not detected toward HD 147888. 

4-5.2. C'2 Excitation 

In diffuse and translucent clouds, most C2 is in various rotational levels of the ground vibrational state 
of the ground electronic state. Levels in upper electronic states are populated by absorption of either 

near-IR photons [to A 1 H U (Phillips system)] or (less frequently) UV photons [to D (Mulliken system) 
or F 1 II It ]. Quadrupole and intersystem (singlet-triplet) transitions then return the molecule to rotational 
levels in the ground vibrational state, whose populations are also affected by collisions with both H and H2 
(van Dishoeck & Black 1982; van Dishoeck & de Zeeuw 1984). 

Because both collisional and radiative processes contribute to the population of the C2 rotational levels in 
diffuse and translucent clouds, care must be taken when deriving the gas physical conditions from observations 
of C2. The effect of the radiative excitation (and subsequent cascades back to the ground electronic state) 
is to increase the populations of the excited rotational levels above their thermal equilibrium values - 
thus yielding excitation temperatures for those levels that are greater than the actual kinetic temperature. 
(While the energies of the excited C2 rotational levels are similar to those of CO for the same J, the relative 
populations in those excited levels are much higher for C2, primarily because C2 has no dipole moment.) 
Because the radiation-induced departures from thermal equilibrium increase with J, it is important (1) to get 
accurate measures of the lowest rotational levels (which yield the best estimates for the kinetic temperature) 
and (2) to measure as many of the higher J levels (typically J > 8) as possible (to determine the extent to 
which radiative processes affect the rotational population distribution and to get tighter constraints on the 
density). 

There are several ways of using the observed lower J populations of C2 to characterize the kinetic 
temperature T^. In principle, the best estimate for Tk comes from T02. Unfortunately, however, T02 can be 
somewhat uncertain, as N is determined from the single, generally relatively weak R(0) line in each of the 
observed C2 bands. For the sight lines analyzed in this study, the relatively well determined X02 range from 
about 20 K (HD 207198) to about 40 K (HD 206267, HD 210121, C Oph); the nominal values for HD 24534, 
HD 27778, and HD 147888 are all somewhat larger (but are also very uncertain). Within the uncertainties in 
the Nj, it is often possible to fit a single excitation temperature to the lowest 3-4 levels (i.e., to J = 0-4 or 
J = 0-6), generally at a value somewhat larger than T02 (e.g., for HD 24534 in Fig. 15). Toward HD 206267, 
for example, the excitation temperatures for the lowest two, three, and four levels are T02 = 36±9 K, T04 
= 51±5 K, and T 6 = 53±3 K, respectively. The populations of the higher J levels, however, generally are 
increasingly higher than would be expected from thermal equilibrium at a local kinetic temperature of T 2 
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— especially for J > 8. As discussed above, such elevated populations of the higher J levels can be produced 
by radiative excitation. 

The set of observed C2 rotational level populations may also be compared with models for the excitation 
in order to estimate both the actual kinetic temperature and the density (for an assumed ambient radiation 
field). Figure 23 shows the normalized excitation diagrams (ln[(5 Nj)/((2J+l) N2)] vs. Ej/k) for six sight 
lines whose rotational level column densities were derived from our fits to the C 2 F-X, D-X, and/or A-X 
bands; the corresponding column densities were given in Table 9. For each sight line, the observed relative 
C2 rotational populations have been compared with those predicted by the models of van Dishoeck & Black 
(1982; see also van Dishoeck & de Zeeuw 1984), for T ranging from 10 to 150 K (in steps of 5 K) and the 
density n ranging from 50 to 1000 cm" 3 (in steps of 50 cm" 3 ) 4 . The models may be characterized by 
the parameter na/I, which is a measure of the relative importance of collisional de-excitation and radiative 
excitation. Here n is the number density of collisional partners (H and H2); a is the (assumed constant) 
collisional de-excitation cross-section for transitions from J +2 to J; and J is a scaling factor for the strength 
of the interstellar radiation field in the near-IR (i.e., near the A-X bands), which is most important for 
the radiative excitation of C2. The cross-section a has not been measured for conditions appropriate for 
diffuse or translucent clouds, but has been estimated to be about 2xl0 -16 cm 2 from a detailed analysis of 
the C Per sight line (van Dishoeck & Black 1982). Calculations by Lavendy et al. (1991), however, suggest 
a somewhat higher value for the de-excitation cross-section for collisions with H 2 (Federman ct al. 1994). 
Higher values of na/I — i.e., higher n and/or lower I — yield predicted rotational populations closer to the 
thermal equilibrium values. The solid curves in Figure 23 show the populations predicted for the best-fit Tk 
and n (assuming 1=1). For that best-fit temperature, the long-dashed straight lines show the rotational 
populations predicted for thermal equilibrium. For the adopted cross-section a, the A-X band /-values noted 
above, and 1=1, the two dotted curves in each panel of Figure 23 correspond to n = 1000 (near the thermal 
equilibrium line) and 100 cm" 3 . The two short-dashed curves correspond to the rotational populations 
predicted for slight differences from the best-fit Tk and n — for {/Ty + 10 , n + 50) and (2\ — 10 , n — 50). 

The values for Tk obtained in this way usually are lower than those for T04, and are more comparable 
to those for T 2- Estimates for T 02 , T 4, and Tk (including values derived from equivalent widths or column 
densities reported in the literature for some additional sight lines) are listed in Table 13. It can sometimes 
be difficult, however, to find a unique, tightly constrained combination of T and na/I for which the observed 
and predicted rotational populations agree within the (often considerable) observational uncertainties - 
especially if data are not available for rotational levels above J = 8. In some cases, "acceptable" fits arc 
possible for a fairly wide range of values, while in others, none of the models seem able to reproduce all 
the observed Nj. Toward HD 206267, for example, the best fit to most of the Nj (for these particular 
models) is for T = 35 K (very similar to the derived T02) and n = n(H) + n(H 2 ) ~ 250 cm" 3 , assuming the 
average interstellar radiation field. Within the uncertainties, however, acceptable agreement between the 
observed and predicted populations is also possible for slightly lower temperatures and densities and also 
for somewhat higher temperatures and densities. [Note that n (which counts collision partners) is smaller 
than «h (which counts hydrogen nuclei) by a factor of (1 — /(H 2 )/2)" 1 ; estimates for both are given in 
columns (7) and (8) of Table 13.] It seems somewhat unlikely, however, that njj toward HD 206267 could 



4 The comparisons were performed using C2 rotational populations calculated via a web-based tool (available at 
http://dibdata.org) developed by B. McCall (sec Adamkovics et al. 2003), which reproduces the relative Nj in Table 7 of 
van Dishoeck & Black (1982) for densities smaller by a factor of 1.35 (to account for differences in the adopted C2 band /- 
values; see van Dishoeck & de Zeeuw 1984). Our derived densities were reduced by an additional factor of 1.2 to account for 
the C2 /-values used in this study. 
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be as high as 600-1200 cm~ 3 , as estimated from the CN/CH ratio (sec § 4.4 above); CN may well trace 
denser gas, on average, than C2. If the (near-IR) radiation field is stronger than average, however, then the 
estimated densities would be proportionally higher. The gas toward HD 206267 traced by H2 probably has 
a somewhat higher temperature ~ 65 K poi(H 2 )] and somewhat lower average density njj = 10-100 cm~ 3 
(obtained from an analysis of O I fine-structure excitation, to be presented elsewhere). Much of the H 2 may 
thus reside in more diffuse gas, with the contributions from the denser clumps masked by blending in the 
lower resolution FUSE spectra. 

Of the ten sight lines with detected C2 (mostly from this paper) listed in the bottom section of Table 13, 
seven (HD 24534, HD 27778, HD 147888, C Oph, HD 204827, HD 206267, and HD 210121) have relative 
rotational populations roughly consistent with T ~ 35-55 K and n ~ 150-300 cm~ 3 (or nn ~ 200-400 cm -3 ). 
On average, the C 2 -containing gas toward HD 207198 appears to be slightly warmer (Tk <~ 60 K) but of 
comparable density; the gas toward HD 73882 appears to be somewhat colder (Tk ~ 20 K) and denser. In 
many of the cases, somewhat lower Tk and slightly lower n and/or somewhat higher Tk and higher n would 
also be possible, within the uncertainties in the Nj. (And recall that since the theoretical curves are for 
values of na/I, changes to the adopted a or I will imply corresponding changes in n.) The values for Tk 
and density in Table 13 are consistent with those estimated by van Dishoeck et al. (1991), for nine sight 
lines in common, except that Tk is higher toward £ Oph and HD 210121 (based on more recent C2 data) 
and lower toward HD 62542. We note that relatively few sight lines reported thus far in the literature have 
sufficiently accurate Nj for high enough J for very strong constraints to be placed on Tk and n; a number of 
these are included in the top part of Table 13. More accurate column densities — especially for J = and 
for the higher J levels (> 8) — would be needed to better constrain the temperature and density in those 
sight lines. 

Comparing the excitation temperatures derived for C 2 with the corresponding values for H 2 (Tables 13 
and 17) indicates that T 2(C 2 ) and Tk(C2) are almost invariably less than T i(H 2 ) — sometimes significantly 
so. For the H2 column density regime sampled by the sight lines with detected C 2 , Toi(H2) is thought to 
be an accurate measure of the local kinetic temperature, while (as noted above) To2(C2) may represent 
only an upper limit to Tk if radiative excitation is significant. The systematically lower excitation and 
kinetic temperatures estimated for C2 suggest that the C2 is preferentially concentrated in somewhat colder, 
denser gas than that traced (on average) by H2 — consistent with predictions of the various models of cloud 
chemistry (e.g., van Dishoeck & Black 1986a). 

Given the indications that many (or most) sight lines are comprised of mixtures of gas with different 
characteristics — and given the differences in physical conditions (e.g., density) often inferred from different 
tracers — we may ask whether such mixtures can be detected and/or characterized via analyses of C2 
rotational excitation. Cccchi-Pcstellini & Dalgarno (2002) have examined the behavior of the integrated 
sight line C2 rotational populations when some fraction of the C2 in the sight line is in dense clumps 
characterized by thermal equilibrium rotational populations. Not surprisingly, the presence of such dense 
clumps is most noticeable when both the temperature and density of the more abundant diffuse gas are low 
- e.g., for Tk < 20 K and n < 100 cm~ 3 . In such cases, the low-J populations (J < 8) for the sight line 
as a whole remain closer to the thermal equilibrium values, and the higher J levels are somewhat enhanced 
as well. Cecchi-Pestellini & Dalgarno conclude, however, that such dense clumps will be very difficult to 
discern via analysis of the C2 excitation (given the typical observational uncertainties) when the bulk of the 
more diffuse gas is warmer than Tk ~ 20 K and/or denser than n ~ 100 cm~ 3 . All of the sight lines in 
our primary sample (with the possible exception of that toward HD 210839) appear to fall in that latter 
category. Several of the sight lines listed in the upper part of Table 13 have Tk < 20 K, but n > 150 cm~ 3 ; 
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the sight line toward HD 110432 has n <~ 100 cm 3 , but uncertain temperature. 

4.6. Isotopic Ratio: 12 C/ 13 C 

The carbon isotopic ratio 12 C/ 13 C in the interstellar medium is often used as an indicator of the degree 
of stellar processing and, hence, of chemical evolution in the Galaxy. It is thought that as the Galaxy evolves 
more 13 C should be produced relative to the dominant isotope 12 C. While the current average local Galactic 
interstellar value of 12 C/ 13 C is about 70 ± 3 (Stahl & Wilson 1992), variations in that ratio of factors of 2-3 
have been observed (Wilson 1999; Casassus, Stahl, & Wilson 2005). Estimates of this ratio have been based 
on observations of the isotopomers of species such as CH + (e.g. Hawkins et al. 1993), CO (Lambert et al. 
1994), H 2 CO (Henkel et al. 1985), and CN (Savage et al. 2002) — each of which presents its own difficulties. 
For example, while CH+ is thought to be unaffected by chemical fractionation, measurement of 13 CH+ is 
often difficult because the (optical) lines can be quite weak. CO, on the other hand, is subject to chemical 
fractionation — either from isotope exchange reactions or from selective photodissociation. Observations of 
mm- wave absorption spectra of the isotopomers of HCO + , HCN, and HNC (which should not be significantly 
fractionated) have yielded 12 C/ 13 C = 59 ± 2 (Lucas & Liszt 1998). 

The models of Warin et al. (1996) predict a somewhat complicated dependence of the 12 CO/ 13 CO 
ratio on density, temperature, and Ay for clouds with 71h between 10 2 and 10 4 cm~ 3 and Ay < 5 (see 
also van Dishoeck & Black 1988). The isotope exchange reaction 13 C+ + 12 CO 13 CO + 12 C+ [which is 
exothermic by only E/k ~ 35 K (Watson, Huntress, & Anicich 1979)] can enhance 13 CO in low temperature 
gas at intermediate densities and Ay < 3 by as much as a factor of 5. At higher temperatures, lower 
densities, and low Ay, however, 13 CO is more readily photodissociated than 12 CO because of its much lower 
abundance, which does not allow effective self-shielding. At higher densities (> 10 5 cm" 3 ) and at higher Ay, 
where chemical reactions dominate the destruction of CO and there is little 13 C + present, the 12 CO/ 13 CO 
ratio is restored to the average value corresponding to the carbon isotopic abundances. Comparison of 
the 12 CO/ 13 CO ratio with the corresponding average 12 C/ 13 C ratio therefore yields indirect information 
concerning the physical conditions in the gas containing CO. 

Reasonably reliable column densities of 13 CO from UV absorption are now available for twelve sight 
lines (six from the literature and six from this paper; Table 16). Of those twelve sight lines, three (p Oph, 
X Oph, C Oph) have 12 CO/ 13 CO ratios significantly above the average value of 70, five (X Per, HD 27778, 
HD 207198, HD 210121, HD 210839) have 12 CO/ 13 CO ratios consistent with the average value, and four 
(HD 73882, 20 Aql, HD 200775, HD 206267) have 12 CO/ 13 CO ratios significantly below the average value. 
Federman et al. (2003) ascribed the higher than average 12 CO/ 13 CO ratios toward p Oph, \ Oph, and ( 
Oph to selective photodissociation and proposed that the ratio A^( 12 CO)//uv, where I\jy gives the relative 
strength of the UV radiation field, could be viewed as a measure of the degree of fractionation. Liszt & Lucas 
(1998) obtained 12 CO/ 13 CO ratios between 15 and 54 (smaller than the average 12 C/ 13 C ratio derived from 
other species) for a small set of mm- wave absorption components, with smaller ratios for higher A^(CO). 

Figure 24 plots the 12 CO/ 13 CO ratios versus the corresponding 12 CO column densities (left) and CN/CH 
ratios (right) for the twelve sight lines in Table 16 with measured N( 13 CO). For comparison, the asterisks in 
the left-hand panel show the total sight line values derived from mm- wave absorption by Liszt & Lucas (1998). 
For this relatively small sample of sight lines, the 12 CO/ 13 CO ratio appears to be anti-correlated with both 
N( 12 CO) (as noted by Liszt & Lucas 1998) and the CN/CH ratio, with correlation coefficients r = —0.647 
and —0.549, respectively (Table 11). If the CN/CH ratio is a density indicator, that anti-correlation would 
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suggest that (in general) the 12 CO/ 13 CO ratio decreases with increasing density — presumably because the 
isotope exchange reaction has enhanced 13 CO in denser, colder gas — consistent with the behavior predicted 
by Warin et al. (1996) for densities in the range 100-1000 cm~ 3 . The excitation and kinetic temperatures 
derived from observations of C2 (Table 17) suggest that T might in fact be low enough in the clouds with 
low values of the 12 CO/ 13 CO ratio for the isotope exchange process to be effective — particularly if the 
effective temperatures are lower for the more centrally concentrated, density-sensitive CN and CO than for 
C 2 . Three of the four sight lines with lower than average 12 CO/ 13 CO ratios (toward HD 73882, HD 200775, 
and HD 206267) have T 2(C2) and/or Tk(C2) < 35 K; the four sight lines with average or above average 
12 CO/ 13 CO ratios and well-determined Tk all have Tk > 45 K. We note, however, that three of the four 
low 12 CO/ 13 CO ratios (toward HD 73882, 20 Aql, and HD 200775) were derived from IUE spectra alone, 
and may be less reliable than the values obtained from the higher resolution, higher S/N HST spectra. As 
there are also indications that the interstellar 12 C/ 13 C ratio may be somewhat variable (e.g., Casassus et al. 
2005), additional data on the 12 CO/ 13 CO ratio and corresponding physical conditions are needed for a more 
secure understanding of the fractionation of CO in diffuse and translucent clouds. Several of the other sight 
lines listed in the upper part of Table 13 with Tk(C2) < 30 K might be interesting candidates for further 
study. 



5. Summary 

Archival UV spectra from HST (GHRS and STIS), FUSE, and/or IUE, together with high and medium 
resolution optical spectra, were analyzed to study the molecular abundances and physical conditions in the 
molecular gas along ten translucent sight lines. Absorption from 12 CO was detected in the far-UV E-X (0-0) 
(1076 A), C-X (0-0) (1087 A), and B-X (0-0) (1150 A) bands, in a number of bands in the A-X system 
(between 1240 and 1550 A), and/or in several of the stronger intersystem bands (between 1360 and 1550 
A); corresponding A-X bands of 13 CO were detected in six of the sight lines. Accurate column densities 
were determined for 12 CO and 13 CO using a combination of curve of growth, apparent optical depth, and 
profile fitting methods. High-resolution optical spectra of K I, CH, and CN were used to determine the 
velocity structure of the gas along the sight lines in order to account for potential saturation effects in the 
UV spectra. 

The UV C 2 F-X (0-0) absorption band at 1341 A was detected in high-resolution STIS echelle spectra of 
HD 24534, HD 27778, and HD 206267; the weaker F-X (1-0) band at 1314 A was detected toward HD 24534; 
and the stronger D-X (0-0) band at 2313 A was detected toward HD 24534, HD 27778, HD 147888, and 
HD 207198. These constitute nearly all of the reported high-resolution spectra of the UV bands of C2. 
Absorption from the optical C2 A-X (3-0) and/or (2-0) bands at 7719 and 8757 A was detected toward 
seven stars. Comparisons among the various C2 bands have yielded a (reasonably) mutually consistent set 
of band /-values; the value adopted for the A-X (2-0) band is consistent with several recent theoretical and 
experimental determinations. There appear to be some discrepancies between the observed and predicted 
wavelengths and/or /-values for (at least) some individual lines in the F-X (0-0) band. 

Reliable column densities (or limits) are now available for 12 CO, CN, CH, C2, and H2 for at least 30 
diffuse or translucent sight lines (11 with detections of 13 CO; 13 with detections of HD; 7 with detections 
of C3; 15 with Ay between 1.0 and 2.7 mag) — enabling investigations of correlations between the various 
species and providing more stringent tests for models of the structure and chemistry of translucent clouds. 
Of those molecular species, CN shows the best correlation with CO — both in overall column density 
and in component structure. Strong correlations between the ratios 12 CO/H 2 and 13 CO/H 2 and the density 
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indicator CN/CH suggest that the abundances of both CO isotopomers increase with increasing local density. 
Differences in the strength of the ambient radiation field (as inferred from H2 rotational excitation) also are 
likely to affect the CO abundances. For our sample of ten translucent sight lines, the iV(CO)/-/V(H2) ratio 
varies over a factor of 100, due primarily to differences in iV(CO). For the total sample, the 12 CO/H2 ratio 
ranges from 2.5 x 10~ 8 to 2.8 x 10~ 5 ; CO is not yet the dominant carbon-containing species in any of the 
sight lines. The models of van Dishoeck & Black (1988, 1989) predict relationships with H 2 , CN, CH, and 
C2 that have slopes similar to those obtained in fits to the observed data, but the predicted iV(CO) generally 
are lower than the observed values by factors of 2-4. The models of Warin et al. (1996) and Le Petit et al. 
(2006) predict much less CO than is observed, for a given iV(H 2 ) — though more of the parameter space 
available in the Meudon PDR code should be explored. Inclusion of the effects of inhomogeneous/clumpy 
structure (e.g., Spaans 1996) and increasing the CO self-shielding and (for the Warin et al. models) the 
gas-phase carbon abundance should yield better agreement with the observed CO abundances. 

Rotational level populations and corresponding excitation temperatures were determined up to J = 3 
for 12 CO (nine sight lines) and up to J = 2 for 13 CO (four sight lines). The 12 CO excitation temperatures 
T01 typically are between 3 and 7 K (well below the likely kinetic temperature of the gas); the corresponding 
values for the higher rotational levels are generally equal to X01, within the uncertainties. Some of these 
excitation temperatures are slightly higher than the values previously reported for eight other comparably 
reddened sight lines, but they are similar to the values found via mm-wave CO absorption by Liszt & 
Lucas (1998) and to the values predicted for diffuse clouds with Ay ~ 1.1 mag by Warin et al. (1996). 
The excitation temperatures Xbi for 13 CO are equal (within the mutual uncertainties) to the corresponding 
values for 12 CO — suggesting that the radiative excitation mechanism proposed by Wannier et al. (1997) 
may not be significant, in many cases. 

Rotational level populations and excitation temperatures were also determined for C2 — up to J = 12 
from the UV F-X (0-0) and (1-0) bands (three sight lines), up to J = 18 from the near-UV D-X (0-0) band 
(four sight lines), and up to J = 12 from the optical A-X (2-0) and (3-0) bands (seven sight lines). The 
excitation temperature for the lowest rotational level (X02) ranges from about 20 to 45 K for those sight lines 
where it is reasonably well determined; corresponding values for the higher rotational levels increase with 
J (especially for J > 8), due to radiative excitation of the upper states. In most cases, the C2 rotational 
excitation seems consistent with model predictions for kinetic temperatures of 35-55 K and local total 
hydrogen densities of 200-400 cm~ 3 , assuming the typical ncar-IR interstellar radiation field. 

Of the twelve sight lines with measured 7V( 13 CO), three have 12 CO/ 13 CO ratios significantly above the 
average Galactic 12 C/ 13 C value of 70, five have 12 CO/ 13 CO ratios consistent with the average value, and 
four have 12 CO/ 13 CO ratios significantly below the average value. There appears to be an anti-correlation 
between the 12 CO/ 13 CO ratio and the CN/CH ratio — suggestive of a decrease in 12 CO/ 13 CO with increasing 
local density, as predicted by the models of Warin et al. (1996) for densities in the range 100-1000 cm~ 3 . 
The 12 CO/ 13 CO ratio is also anti-correlated with the column density of CO, as found (for a smaller sample) 
by Liszt & Lucas (1998). Sight lines with below average 12 CO/ 13 CO ratios also appear to have kinetic 
temperatures (from analysis of the C2 rotational excitation) less than about 35 K; sight lines with average 
or above average 12 CO/ 13 CO ratios appear to be somewhat warmer. While 12 CO/ 13 CO ratios above the 
average value likely result from selective photodissociation (Federman et al. 2003), which favors the self- 
shielded 12 CO, the below-average ratios probably are due to isotope exchange reactions, which favor 13 CO 
in cold gas (T < 35 K). 

While the H 2 excitation temperature T i ranges from about 50 K to 110 K for our ten sight lines, there 
are several indications that C2 and (especially) CO trace colder, denser gas in many cases. The excitation 
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temperature T02 of C2 (which provides an upper limit to the local kinetic temperature for the gas containing 
the C2) and the Tk inferred from comparisons with models of the excitation are generally lower than T i 
for H2 (which should equal the kinetic temperature of the gas traced by H2). Fits to the CO line profiles 
indicate that smaller 6-values (0.3-0.7 km s _1 ) are needed than for fits to the profiles of K I, CH, and CN; 
the CO structure also appears more similar to that of CN (which is concentrated in denser regions). The 
column densities of both 12 CO and 13 CO are most tightly correlated with that of CN, and their abundances, 
relative to H2, are correlated with the density indicator 7V(CN)/7V(CH). 

These results illustrate the utility of combining information from multiple species sensitive to different 
combinations of physical parameters in order to probe regions characterized by different properties (along 
a sight line or within a single "cloud"). They also emphasize the need for developing models capable of 
incorporating ranges and/or mixtures of physical conditions in order to accurately assess and predict the 
abundances of various atomic and molecular species in the diffuse ISM. 

This work is based in part on data obtained for the Guaranteed Time Team by the NASA-CNES- 
CSA FUSE mission operated by the Johns Hopkins University. Financial support to U. S. participants has 
been provided by NASA contract NAS5-32985. P. Sonnentrucker acknowledges partial support from NASA 
LTSA grant NAG5-13114 to the Johns Hopkins University. D. E. Welty acknowledges support from NASA 
LSTA grants NAG5-3228 and NAG5-11413 to the University of Chicago. We thank E. Jenkins for providing 
processed STIS spectra of HD 206267 and HD 210839, S. Federman and K. Pan for providing results from 
their study of the Cep OB2 and OB3 regions in advance of publication, S. Federman for comments on an 
earlier version of the paper, B. McCall and T. Oka for advice on the analysis of C2, and the anonymous 
referee for clarifying several points in the discussion. This work has used the profile fitting procedure Owens. f 
developed by M. Lemoine and the FUSE French Team. 

Facilities: FUSE, HST (GHRS, STIS), IUE, ARC, KPNO:coude feed, McD:2.7m 

A. Component Structures 

Table 14 lists the component structures derived from high-resolution optical spectra of K I, CH, and CN 
that have been used to model the profiles of C2 and CO in this paper. The K I structures for HD 206267, 
HD 207198, and HD 210839 were originally listed by Welty & Hobbs (2001); they differ slightly from those 
adopted by Pan ct al. (2004). The K I and CH structures for HD 73882, HD 185418, and HD 192639 were 
used by Snow et al. (2000) and by Sonnentrucker et al. (2002, 2003) in studies of those individual lines of 
sight. Because the total sight line column densities of K I and CH are generally very tightly correlated, and 
because the absorption-line profiles of the two species often are strikingly similar in detail (Welty & Hobbs 
2001), we have used the more complex structure discernible in the K I profiles to model the CH profiles. 
Because CN is likely to be more concentrated in the cooler, denser regions of each sight line, however, we have 
performed independent fits to the CN profiles. The CN structures for HD 73882 and HD 210121 are based 
on analyses of spectra originally reported by Grcdel ct al. (1991) and Welty & Fowler (1992), respectively. 
A full description of the optical data and profile fits — for these and other translucent sight lines — will be 
given by Welty, Snow, & Morton (in preparation). 
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B. C2 Wavelengths and Line Strengths 

The energy levels of the various rotational levels in the ground vibrational state of the ground X 1 H£ 
electronic state of C2 were calculated via the usual formula E* =0 (J) = Bq J {J + 1) — J 2 (J + l) 2 , using 
values for the constants and D* from Chauville, Maillard, & Mantz (1977). Because C2 is homonuclear, 
with zero nuclear spin, only even J levels are possible in that ground electronic state. The energy levels 
of the rotational levels of the v = and v = 1 vibrational states of the F 1 Ii u and A 1 !^ upper electronic 
states were calculated using the same formula, using constants from Herzberg et al. (1969) and Chauville et 
al. (1977), respectively. The expected wavelengths are then just the reciprocals of the differences between 
energy levels (in appropriate units); the values for the individual rotational lines in the F-X (1-0) and (0-0) 
(1314 and 1341 A) bands and in the A-X (3-0) and (2-0) (7719 and 8757 A) bands arc listed in columns 
(3), (4), (8), and (9) of Table 15. In the same way, the energy levels of the rotational levels in the v = 
vibrational state of the upper electronic state were calculated using constants from Sorkhabi et al. 

(1998), with i/Q = 43227.2 cm -1 estimated by comparing predicted wave numbers with those measured from 
their laser-induced fluorescence spectrum. The corresponding wavelengths for the rotational transitions in 
the D-X (0-0) band at 2313 A are listed in column (6) of Table 15; these predicted values are in better 
agreement with the D-X (0-0) lines observed toward HD 24534 than are the measured wave numbers listed 
by Sorkhabi et al. The relative intensities of the individual F-X and A-X rotational transitions are given by 
the usual formulas >Sr(j) = 2(27+1) ' ^Q(J) = I' an< ^ ^ p ( J ) = 2(27+1) ' ^ nc corres P on ding relative strengths 
for the lines in the D-X (0-0) band (which has no Q branch) are given by Sr(j) = 2J+1 an ^ &p(J) = 27+1 )■ 



C. Compilation of CO and Other Molecular Data 

In order to compare the column densities of 12 CO and 13 CO obtained in this paper with those found 
for other Galactic sight lines, we have collected CO column densities derived from UV absorption- line data 
from the literature. Some of the values for sight lines in the Copernicus survey of Federman et al. (1980) 
[most with relatively low jV(CO)] have recently been revised by Crenny & Federman (2004). In the past 
few years, results from studies of more heavily reddened sight lines [generally with higher AT(CO)], based on 
spectra obtained with IUE, HST, and/or FUSE, have become available (e.g., Lambert et al. 1994; Knauth 
et al. 2001; Shcffcr et al. 2002a; Federman et al. 2003, 2004; Pan et al. 2005). In this paper, we have 
analyzed IUE, FUSE, and/or HST spectra of additional translucent sight lines in an attempt to enlarge 
the sample at higher iV(CO). Table 16 lists total column densities for H2, CH, C2, CN, and CO for sight 
lines where CO has been measured. 5 This compilation is thus similar to that of Federman et al. (1994a), 
except that sight lines with no information on CO have not been included. Moreover, we have generally not 
included sight lines whose CO column densities were derived from COG analyses of equivalent widths from 
single IUE spectra (in view of the low resolution and S/N of the data), particularly if no constraints from 
higher resolution optical or UV spectra were used or if the CO band structure was not used in computing 
the theoretical COG. Of the sight lines in common with the earlier compilation of Federman et al. (1994a), 
many of those with CO column densities below 10 14 cm~ 2 have revised values that are significantly smaller; 
some with column densities above 10 14 cm~ 2 have revised values different by factors > 2. 

The C2 column densities in Table 16 have been adjusted for differences in the adopted /-values — e.g., 



5 The column densities in Table 16 are included in a compilation of atomic and molecular column densities maintained at 
http://astro.uchicago.edu/~welty/coldens.html. 
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by -0.15 dex for references using /[A-X(2-0)] = 0.001 (Federman et al. 1994a; Thorburn et al. 2003). Where 
the C2 reference is in parentheses, the total column density is based on comparing the Nj derived from the 
equivalent widths in that reference with theoretical C2 rotational distributions (see § 3.3.2 and § 4.5.2). 

Table 17 lists color excesses, visual extinctions, total hydrogen column densities, molecular fractions, 
excitation temperatures for H 2 and C2, as well as several ratios derived from the column densities listed in 
Table 16. 
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Fig. 2.- STIS spectra of 12 CO A-X bands toward HD 27778 and HD 147888. At the top arc higher 
resolution optical spectra of K I, CH, and CN, which show more detailed component structure; the CN 
profile toward HD 147888 was calculated from the component structure in Pan et al. 2004. CO rotational 
levels are marked above each spectrum (R, Q, and P branches are left-to-right). 
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Fig. 3.— STIS spectra of 12 CO A-X bands toward HD 185418 and HD 206267. At the top are higher 
resolution optical spectra of K I, CH, and CN, which show more detailed component structure. CO rotational 
levels are marked above each spectrum (R, Q, and P branches are left-to-right). 
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Fig. 4.— STIS spectra of 12 CO A-X bands toward HD 207198 and HD 210839. At the top are higher 
resolution optical spectra of K I, CH, and CN, which show more detailed component structure. CO rotational 
levels are marked above each spectrum (R, Q, and P branches are left-to-right). 
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Fig. 5.— STIS spectra of the 13 C0 A-X (7-0) (1347 A) band toward four stars and the 13 C0 A-X (4-0) 
(1421 A) band toward HD 24534. Rotational levels are marked above each spectrum; lines from J = 2 are 
evident in the A-X (4-0) band. Additional weak absorption red-ward of the A-X (4-0) band is from 12 C 18 
(Sheffer et al. 2002b). 
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Fig. 6.— WE spectra of 12 CO A-X bands toward HD 27778 and HD 73882. CO rotational levels are marked 
above each spectrum (R, Q, and P branches are left- to-right). Additional absorption from 13 CO A-X bands 
and strongest 12 CO intersystem bands is also marked. Note the larger velocity range covered in this plot, 
compared to plots of STIS spectra. 
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Fig. 7.— Curves of growth for 12 CO toward HD 27778 and HD 207198. Solid lines are theoretical curves 
[single component; using structure for A-X (6-0) band] for b = 0.3, 0.5, 1.0, 3.0, and 5.0 km s _1 (bottom to 
top). Dotted lines are corresponding theoretical curves for the E-X (0-0) band, displaced downward by 1.0 
for clarity. Filled circles represent values for permitted A-X bands; open circles represent intersystem bands, 
triangles represent B-X, C-X, and E-X (0-0) bands; the point for the E-X (0-0) band appears on both sets 
of theoretical curves. Column density and (approximate) effective &-value are given for each sight line. 
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Fig. 8.— Fits to UV line profiles toward HD 27778. Observed spectra (from KPNO coude feed, STIS, and 
IUE) are given by solid histograms; fits to the UV profiles are given by dotted lines. Note the 12 CO a'14 
intersystem band centered at v ~ 120 km s^ 1 near the permitted A-X (4-0) band. Individual components are 
marked for higher resolution spectra of K I, CH, and CN R(0) lines; note also the CN R(l) and P(l) lines. 
CO rotational lines (R, Q, and P branches; J = 0-3) are marked at the mean velocities of the unresolved 
component blends; two components of roughly equal column density, separated by 2.3 km s _1 , were used in 
fitting the CO profiles. 
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Fig. 9.— Fits to UV line profiles toward HD 210839. Observed spectra (from KPNO coude feed and STIS) 
are given by solid histograms; fits to the UV profiles are given by dotted lines. Individual components are 
marked for higher resolution K I, CH, and CN lines. The CH profile is based on the component structure in 
Crane et al. (1995). CO rotational lines (R, Q, and P branches; J = 0-3) are marked at the mean velocities 
of the unresolved component blends; two components, separated by 2.1 km s _1 , were used in fitting the CO 
profiles. 
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Fig. 10.— UV C 2 F-X (1-0) and (0-0) (1314, 1341 A) bands toward HD 24534, HD 27778 and HD 206267. 
Spectra have been shifted so that R(0) lines are at v = km s _1 ; note the expanded vertical scale for the top 
three profiles. Observed spectra (from STIS) are given by solid histograms; fits to the profiles are given by 
dotted lines. C2 rotational lines (R, Q, and P branches; J = 0, 2, 4, ...) are marked at the mean velocities 
of the unresolved component blends. Note the poor fit to the F-X (0-0) R-branch "pile-up" for J = 6-12 
near —50 km s -1 (toward all three stars) and the relative weakness of the lines from J = 8 and 10 in the F-X 
(0-0) band [relative to those in the F-X (1-0) band] toward HD 24534. 
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Fig. 11— C 2 F-X (0-0) (1341 A) band toward HD 24534 (combined spectra from three epochs). The 
spectrum has been shifted so that the R(0) line is at v — km s -1 ; note the expanded vertical scale. The 
observed spectrum (from STIS) is given by the solid histogram; a fit to the spectrum is given by the dotted 
line. C2 rotational lines (R, Q, and P branches; J = 0, 2, 4, ...) are marked at the mean velocities of 
the unresolved component blends. The solid tick marks have all been shifted by 8 km s _1 , relative to the 
predicted wavelengths (see Sec. 3.3.2); the dotted tick marks for J > 10 show the unshiftcd locations. In 
the fit, only the lower J lines (J < 10) have been shifted. Note the somewhat better fit to the R-branch 
"pile-up" and the weak absorption features at the unshifted positions of the higher J Q-branch lines. 
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Fig. 12.— UV C 2 D-X (0-0) (2313 A) band toward HD 24534, HD 27778, HD 147888, and HD 207198. 
Spectra have been shifted so that R(0) lines are at v = km s _1 ; note the expanded vertical scale. Observed 
spectra (from GHRS or STIS) are given by the solid histograms; fits to the profiles are given by the dotted 
lines. C 2 rotational lines (R and P branches; J = 0, 2, 4, ...) arc marked at the mean velocities of 
the unresolved component blends; note the three components (or component groups) discernible toward 
HD 207198. 
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Fig. 13. — Optical C2 A-X (2-0) band at 8757 A toward seven stars. Spectra have been shifted so that the 
R(0) line is at v = km s _1 ; note the expanded vertical scale for all profiles. Sections of the spectrum of 
HD 24534 within the stellar H Paschen 12 (emission) feature were not normalized (several R-branch lines 
are thus missing here). C2 rotational lines (R, Q, and P branches; J = 0, 2, 4, ...) are marked at the mean 
velocities of the unresolved component blends. The spectrum of HD 204827, with much stronger C2 lines, 
is included to show the band structure more clearly. Observed spectra (from ARCES) are given by solid 
histograms; fits to the profiles are given by dotted lines. 
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Fig. 14. — Optical C2 A-X (3-0) band at 7719 A toward five stars. Spectra have been shifted so that the 
R(0) line is at v = km s _1 ; note the expanded vertical scale for all profiles. C2 rotational lines (R, Q, and 
P branches; J — 0, 2, 4, ...) arc marked at the mean velocities of the unresolved component blends. The 
spectrum of HD 204827, with much stronger C2 lines, is included to show the band structure more clearly. 
Observed spectra (from ARCES) are given by solid histograms; fits to the profiles are given by dotted lines; 
note the additional broad absorption feature (a weak diffuse interstellar band) nearly coincident with the 
Q(2) line in several cases. 
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Fig. 15. — Rotational populations of C2 (J = 0-14) toward HD 24534. The upper panel shows the normalized 
Nj/{2J+l) obtained from the UV bands — F-X (1-0) at 1314 A (upright triangles), F-X (0-0) at 1341 A 
(squares), and D-X (0-0) at 2313 A (inverted triangles) — compared with the average values derived from 
all the UV and optical bands (filled circles). Note the apparent deficits for J = 8 and 10 in the F-X (0-0) 
band; similar differences are seen for those levels in the F-X (0-0) bands toward HD 27778 and HD 206267. 
The lower panel shows the corresponding Nj/(2J+1) obtained from the optical bands — A-X (3-0) at 7719 
A (triangles) and A-X (2-0) at 8757 A (squares) — again compared with the average values (filled circles). 
In both panels, the rotational levels J are noted at the bottom. The dotted lines show the extrapolated 
straight line fit to the average Nj for the lowest three levels; the corresponding excitation temperature, XJh 
~ 74 K, is the negative of the inverse of the slope of that line. The higher rotational levels are characterized 
by progressively higher excitation temperatures. 
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Fig. 16.— iV( 12 CO) {left) and A^( 13 CO) {right) vs. N(tt 2 ). Circles denote CO data from HST, squares from 
FUSE, triangles from IUE, and crosses from Copernicus. Filled symbols are for sight lines with measured 
13 CO. The dotted lines show the predictions from models H1-H6, T1-T6, and 11-18 of van Dishoeck & Black 
(1988), for 0.5, 1.0, and 10 times the average Galactic radiation field (left to right). Letters "W" show the 
results the for "standard" diffuse, translucent, and dark cloud models of Warin et al. (1996). The dashed 
line (for 12 CO only) shows the predictions from the Meudon PDR code (Le Petit et al. 2006) for iih = 
100 cm -3 and the average Galactic field. 
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Fig. 17. — iV(C2) vs. ./V(H 2 ). The solid lines are linear least-squares fits (weighted and/or unweighted) to 
the data. The dotted line shows the predictions from models T1-T6 of van Dishoeck & Black (1988). 
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Fig. 18.— iV( 12 CO) vs. 7V(CN) (/e/t) and vs. iV(CH) (n^Af). Circles denote CO data from HST, squares 
from FUSE, triangles from IUE, and crosses from Copernicus. Filled symbols are for sight lines with measured 
13 CO. The solid lines are linear least-squares fits (weighted and/or unweighted) to the data. The dotted 
lines are the predictions from models T1-T6 of van Dishoeck & Black (1989). The dashed lines show the 
predictions from the Mcudon PDR code (Lc Petit et al. 2006) for diffuse clouds [Ay = 1.0, jih = 20- 
2000 cm -3 (left to right), and the average Galactic field]. 
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Fig. 19.— iV( 12 CO) vs. N(C 2 ) (left) and vs. N(K I) {right). Circles denote CO data from HST, squares from 
FUSE, triangles from IUE, and crosses from Copernicus. Filled symbols are for sight lines with measured 
13 CO. The solid lines are linear least-squares fits (weighted and/or unweighted) to the data. The dotted line 
for C2 shows the predictions of models T1-T6 of van Dishoeck & Black (1989). 
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Fig. 20.— A^( 12 CO)/7V(H 2 ) (left) and N ( 13 CO) / N (R 2 ) (right) vs. iV(H 2 ). Circles denote CO data from 
HST, squares from FUSE, triangles from IUE, and crosses from Copernicus. Filled symbols are for sight lines 
with measured 13 CO. Dotted lines show predictions from models H1-H6, T1-T6, and 11-18 of van Dishoeck 
& Black (1988), for 0.5, 1.0, and 10 times the average Galactic radiation field (left to right). Letters "W" 
show the results for the "standard" diffuse, translucent, and dark cloud models of Warin et al. (1996). The 
dashed line (for 12 CO only) shows the predictions from the Meudon PDR code (Le Petit et al. 2006) for n H 
= 100 cm~ 3 and the average Galactic field. 




log[N(CN)/N(CH)] log[N(CN)/N(CH)] 

Fig. 21.— Column density ratios 12 CO/H 2 (left) and 13 CO/H 2 (right) vs. density indicator CN/CH. The 
CO data are from HST (circles), FUSE (squares), IUE (triangles), and Copernicus (crosses); filled symbols 
are for sight lines with measured 13 CO. Solid lines arc linear least-squares fits (weighted and/or unweighted) 
to the observed data. Dotted line for 12 CO/H 2 shows the predictions from models T1-T6 of van Dishoeck 
& Black (1989). The letter "M" (for 12 CO only) shows the prediction from the Meudon PDR code (Le Petit 
et al. 2006) for a diffuse cloud with Ay = 1.0, nn = 100 cm~ 3 , and the average Galactic field. 
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Fig. 22. — Rotational excitation of 12 C0 toward eight stars. Open squares show the normalized 12 C0 
rotational level column densities, for J — 0-3. Dotted lines show the extrapolated straight line fits to J = 
0-1. In each case, the excitation temperature Xbi is the negative of the inverse of the slope of that line. 
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Fig. 23. — Relative C2 rotational population distributions toward six stars, as functions of excitation energy, 
compared to predicted values. In each case, the solid line represents the best-fit theoretical model (to the 
nearest 5 K in Tk and 50 cm~ 3 in n), based on the analysis of van Dishoeck & Black (1982) but adjusted 
for our adopted C2 /-values. The straight long-dashed line shows the thermal equilibrium distribution for 
the best-fit T<; the two dotted curves show the distributions for that Tk and n = 1000 cm~ 3 (nearer to the 
equilibrium line) and 100 cm~ 3 . The two short-dashed curves show the distributions for (Tk — 10, n — 50) 
and (Tk + 10, n + 50). The C2 rotational distributions yield well determined Tk and n in all six cases. 
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Fig. 24.— Column density ratio 12 CO/ 13 CO vs. N(CO) (left) and vs. CN/CH (right). The UV absorption- 
line data are from HST (circles), FUSE (squares), and IUE (triangles); mm- wave absorption-line data (aster- 
isks) are from Liszt & Lucas (1998). The solid lines are linear least-squares fits (weighted and unweighted) 
to the observed data. The dotted horizontal line indicates the average Galactic value (70) for the 12 C/ 13 C 
ratio. 
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Table 1. Sight lines 



Star 


Type 


V 


E(B - V) 


Av 


JV(H I) 


JV(H 2 ) 


I 


b 


d 


v v 


vsini 












(cm- 2 ) 


(cm" 2 ) 






(PC) 


(km s- 1 ) 


(km s- 1 ) 


HD 24534 


09.5pe 


6.1 


0.59 


2.05 


20.72 


20.92 


163.08 


-17.14 


590 


15.0 


150 


HD 27778 


B3 V 


6.36 


0.37 


1.01 


20.98 


20.79 


172.76 


-17.39 


223 


15.7 


105 


HD 73882 


08 V 


7.22 


0.72 


2.36 


21.12 


21.11 


260.18 


0.64 


1100 


22.0 




HD 147888 


B5 V 


6.74 


0.47 


1.91 




20.48 


353.65 


17.71 


136 


-4.0 


180 


HD 185418 


B0.5 V 


7.45 


0.51 


2.03 


21.12 


20.76 


53.60 


-2.17 


950 


-2.0 




HD 192639 


08c 


7.12 


0.66 


1.87 


21.32 


20.69 


74.90 


1.48 


1100 


-6.0 


120 


HD 206267 


06f 


5.62 


0.53 


1.37 


21.31 


20.86 


99.29 


3.74 


1000 


-7.8 


80 


HD 207198 


09 He 


5.95 


0.62 


1.36 


21.34 


20.83 


103.14 


6.99 


1000 


-18.4 




HD 210121 


B3 V 


7.67 


0.38 


0.83 


20.63 


20.75 


56.88 


-44.46 


210 


-8.0 


<30 


HD 210839 


06 If 


5.04 


0.57 


1.57 


21.16 


20.84 


103.83 


2.61 


505 


-74.0 


200 



Note. — Type, V, E(B - V), and d from Thorburn et al. 2003; A v , N(H I), and JV(H 2 ) from Rachford et al. 2002 (and 
references therein) or Cartlcdgc ct al. 2004 (HD 147888); v r and vsini from SIMBAD or Welty & Fowler 1992 (HD 210121). 
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Table 2. UV Data Sets 



Star 


Mission 


Instrument 


Data Sets 


HD 24534 


HST 


STIS 


c-64812010-030, o64813010-020, o66p01020, o66p02010 




HST 


GHRS 


z2do020at 


HD 27778 


IV E 


SWP 


21264-66, 21284-85, 21289-92, 21297-98, 24095 




HST 


STIS 


o59s01020, o61b01010-040 




FUSE 




pll60301 


HD 73882 


IV E 


SWP 


25437, 30064-65, 30075-77, 38299, 38311 




FVSE 




pll61301, pll61302, X0211301 


HD 147888 


HST 


STIS 


o59s05010, o61b04010-040 




FVSE 




pll61501 


HD 185418 


HST 


STIS 


o5c01q010 




FVSE 




pll62301 


HD 192639 


HST 


STIS 


o5c08t010 




FVSE 




pll62401 


HD 206267 


HST 


STIS 


o51h09010-040 




FVSE 




pll62701 


HD 207198 


IV E 


SWP 


17472, 25981-87, 26001, 26017 




HST 


STIS 


o59s06010-020, o61b06010-040 




FVSE 




pll62801 


HD 210121 


IV E 


SWP 


38860, 38872, 38886, 38949 




HST 


GHRS 


z2m2010ht, z2m22010kt 




FVSE 




pll63001 


HD 210839 


HST 


STIS 


O54304010 




FVSE 




pll63101 



Note. — All the FVSE data were obtained for programs P116 or X021 (T. Snow, P.I.)- The STIS 
data were obtained for HST GO programs 8273 (HD 27778, HD 147888, HD 207198; D. Meyer, 
P.I.), 8622 (HD 24534; D. Lambert, P.I.), 8623 (HD 24534; T. Snow, P.I.), 9465 (HD 27778, 
HD 147888, HD 207198: U. J. Sofia, P.I.), GTO programs 8043 and 8484 (HD 206267, HD 210839; 
E. Jenkins, P.I.), or SNAP program 8241 (HD 185418, HD 192639; J. Lauroesch, P.I.). The GHRS 
data were obtained for HST GO programs 5487 (HD 24534; T. Snow, P.I.) and 5739 (HD 210121; E. 
van Dishoeck, P.I.). The IVE data were obtained for programs IDHTS, IEETS, IMFTS, IMGTS, 
IMITS, and IMLTS (HD 27778, HD 73882, HD 207198; T. Snow, P.I.) and program IMLLH 
(HD 210121; L. Hobbs, P.I.). 
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Table 3. CO Equivalent Widths (FUSE and STIS) 



Species 


Band 


A a 


log /A a 


27778 


147888 


185418 


192639 


206267 


207198 


210839 


12 CO 


E-X (0-0) 


1076.034 


1.864 


83.9±7.3 




64.6±4.7 


38.9±5.5 


119±18 


124±15 


84.5±9.0 




C-X (0-0) 


1087.868 


2.126 


blend 




96±20 


59±10 


blend 


blend 


blend 




B-X (0-0) 


1150.483 


0.887 


69.9±5.0 








65.0±4.7 




39.6±2.6 




A-X (12-0) 


1246.059 


-0.950 


11.9±3.2 








11.2±4.3 




6.5±1.9 




A-X (11-0) 


1263.433 


-0.643 


21.1±3.4 


4.9±2.8 






26.3±6.5 


9.1±2.5 


6.9±2.0 




A-X (10-0) 


1281.866 


-0.279 


40.2±5.0 


11.3±2.6 


3.5±0.7 




36.0±6.3 


17.3±2.7 


14.3±2.0 




A-X (9-0) 


1301.403 


0.092 


52.5±5.1 


14.2±5.2 


8.1±1.9 




62.0±7.7 


29.9±2.5 


23.3±1.9 




A-X (8-0) 


1322.150 


0.423 


74.3±4.3 


21.4±4.6 


13.9±2.0 


3.3±2.4 


81.0±6.2 


46.5±2.2 


30.5±3.1 




A-X (7-0) 


1344.186 


0.745 


97.0±4.9 


33.9±4.8 


20.1±2.0 


6.8±1.9 


103.0±6.0 


71.8±1.9 


44.5±3.5 


13 co 


A-X (8-0) 


1325.793 


0.428 


3.7±1.4 








7.3±2.0 








A-X (7-0) 


1347.575 


0.747 


9.7±1.9 








14.5±3.0 


3.9±1.6 


1.5±0.5 



Note. — Values are equivalent width (mA) and 1-er uncertainty. Equivalent width for the B-X (0-0) band toward HD 73882 is 
63.9±5.5 mA. 

a Wavelengths and /-values are from Morton & Noreau (1994) for the A-X bands and from Federman ct al. (2001) for the E-X, C-X, 
and B-X bands. 



-64- 



Table 4. CO A-X Band Equivalent Widths (IUE and HST) 



Species 


Band 


A a 


log /V 


149757 b 


24534° 


27778 d 


73882 s 


207198 f 


210121 s 


12 CO 


(12-0) 


1246.059 


-0.950 




11.9±1.1 


11.9±3.2 










(11-0) 


1263.433 


-0.643 


5.4±0.7 


18.6±1.5 


21.1±3.4 




9.1±2.5 






(10-0) 


1281.866 


-0.279 


10.2±0.9 


28.9±1.6 


40.2±5.0 




17.3±2.7 






(9-0) 


1301.403 


0.092 


15.7±1.2 


40.8±1.6 


52.5±5.1 




29.9±2.5 






(8-0) 


1322.150 


0.423 


24.6±1.3 


55.1±1.4 


73± 6 


110±15 


50± 6 
















74.3±4.3 




46.5±2.2 






(7-0) 


1344.186 


0.745 


33.9±1.6 


72.3±2.1 


98± 9 


90±14 


65± 8 
















97.0±4.9 




71.8±1.9 






(6-0) 


1367.623 


1.036 


42.5±1.3 


88. Oil. 4 


>103 


134±20 


78± 8 






(5-0) 


1392.525 


1.305 


48.8±1.2 


107.7±2.3 


124± 8 


141±17 


119± 8 


105±15 




(4-0) 


1419.044 


1.520 


54± 7 


122.8±2.0 


152± 7 




149± 6 


115±13 




















97±14 




(3-0) 


1447.352 


1.700 


70.4±1.1 


139.6±2.0 


180± 9 


174±12 


170± 8 






(2-0) 


1477.565 


1.773 


74.4±1.2 




183±12 


189±12 


181± 6 


119±15 




















119±12 




(1-0) 


1509.748 


1.643 






158± 8 


176±10 


176± 7 


142±10 




el2 


1323.1 


-1.697 




3.09±0.24 












a'17 


1366.214 


-1.383 


1.06±0.21 


6.38±0.29 




18± 6 








dl2 


1366.443 


-1.143 


1.82±0.22 


9.62±0.23 




21± 7 








a' 14 


1419.535 


0.072 


17.46±0.57 


44.67±0.40 


54± 8 




(13± 3) 


46±10 




















37±11 




e5 


1449.5 


-0.952 


2.36±0.17 


15.99±0.23 


<11 


22± 5 


<12 






a'll 


1480.784 


-1.347 


1.43±0.24 


6.52±0.39 


<14 


17± 6 


<12 






el 


1543.1 


-0.168 






46± 6 


65±10 


35±20 


(64±11) 


13 CO 


(8-0) 


1325.793 


0.428 




5.0±0.8 


3.7±1.4 


26± 8 








(7-0) 


1347.575 


0.747 


<1.0 


8.6±0.7 


14± 5 


35± 8 


8± 4 
















9.7±1.9 




3.9±1.6 






(6-0) 


1370.617 


0.951 


<1.8 


14.3±0.8 


17± 6 


49±12 


<12 






(5-0) 


1395.229 


1.306 


3.6±0.3 


17.1±0.9 


27± 6 


61±12 


10± 4 






(4-0) 


1421.340 


1.537 


6.3±2.1 


24.8±1.0 


38± 6 


67± 8 


8± 5 


(51±12) 




















17± 6 




(3-0) 


1449.145 


1.701 


10.3±0.8 




33± 6 


68± 9 


17± 4 






(2-0) 


1478.793 


1.774 


11.9±0.6 




46± 7 


64±10 


16± 4 






















16± 5 




dl2 


1370.928 


0.375 




4.59±0.19 











Note. — Values for each sight line are equivalent width (mA) and 1-a uncertainty (or 3-tr limit). Integer values in 
the last four columns are from IUE or GHR.S data; decimal values are from STIS data. 

a Morton & Norcau 1994 (permitted bands); Shcffcr ct al. 2002a or Eidclsbcrg & Rostas 2003 (intcrsystcm bands) 

b Federman ct al. 1994b; Lambert ct al. 1994 

C STIS data (permitted bands); Sheffcr ct al. 2002a (intcrsystcm bands) 

d IUE data (A > 1320A); STIS data (A < 1360A); values for stronger 13 CO bands may contain some contribution 
from 12 C ls O 

e IUE data; values for stronger 13 CO bands may contain some contribution from 12 C ls O 
f IUE data (A > 1320A); STIS data (A < 1360A) 

S IUE data (Welty & Fowler 1992; first line); GHR.S data (second line) 
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Tabic 5. Total 12 CO and 13 CO Column Densities. 



Star 






N( 12 CO) 








N( 


13 CO) 






COG a 




AOD 


FIT 


Adopt b 


COG* 


AOD 


FIT 


Adopt b 


HD 24534 


16.15±0.15 


0.8 


>16.02 


16.01±0.09 


16.15±0.07 


14.25±0.05 


14.18±0.04 


14.29±0.03 


14.29±0.02 


HD 27778 


16.10±0.20 


1.0 


>16.03 


16.08±0.03 


16.08±0.03 


14.20±0.10 


14.16±0.09 


14.19±0.10 


14.19±0.10 


HD 73882 


16.6±0.2 


0.8 




16.55±0.20 


16.55±0.20 


15.2±0.2 




15.15±0.20 


15.15±0.20 


HD 147888 


15.30±0.15 


0.5 


15.30±0.11 


15.30±0.08 


15.30±0.08 






<13.5 


<13.5 


HD 185418 


14.70±0.15 


1.0 


14.71±0.06 


14.75±0.03 


14.75±0.03 










HD 192639 


14.15±0.15 


1.0 


14.04±0.12 


14.13±0.13 


14.13±0.13 










HD 206267 


16.10±0.25 


1.0 


16.06±0.12 


16.04±0.04 


16.04±0.04 


14.40±0.10 


14.40±0.10 


14.39±0.03 


14.39±0.03 


HD 207198 


15.45±0.25 


2.0 


15.52±0.08 


15.52±0.04 


15.52±0.04 


13.70±0.10 


13.78±0.17 


[13.75±0.10] 


[13.75±0.10] 


HD 210121 


16.00±0.20 


0.5 


>15.51 


15.83±0.08 


15.83±0.08 


13.70±0.15 


>13.67 


13.82±0.10 


13.82±0.10 


HD 210839 


15.75±0.25 


0.5 


15.44±0.08 


15.46±0.06 


15.46±0.06 


(13.30±0.15) 


(13.36±0.12) 


[13.59±0.08] 


[13.59±0.08] 



Note. — Logarithmic column densities and 1-cr uncertainties (or 3-cr limits); units arc cm -2 . Methods used to derive the column densities 
arc curve of growth (COG), apparent optical depth (AOD), and profile fitting (FIT). Values in parentheses are for TVo only. Values for 13 CO in 
square braces are for No + JVi, assuming that Ni/Nq, is the same for 13 CO as for 12 CO. 

a Curvcs of growth assume T cx — 3 K (for HD 185418, HD 192639, and HD 207198) or T cx — 5 K (for all others). 

b Valucs from Shcffer ct al. (2002a, 2002b) arc adopted for HD 24534. FIT values arc adopted for all other sight lines. 



Table 6. CO Rotational Level Column Densities 



Star 




Af( 12 CO) 


x 10~ 14 




JV( 


13 CO) x 10 


-13 




J = 


J = 1 


J = 2 


.7 = 3 


J = 


J = 1 


J = 2 


HD 24534 


48.8±17.8 


43.8±9.3 


9.2±1.4 


0.67±0.30 


11.9±1.3 


6.7±0.4 


0.70±0.08 


HD 27778 


52.1±5.9 


53.5±4.9 


12.5±1.4 


1.17±0.26 


8.2±1.2 


6.8±1.3 




HD 147888 


5.4±2.1 


9.8±3.2 


4.0±0.9 


0.88±0.20 








HD 185418 


3.5±0.3 


2.0±0.2 


0.21±0.14 










HD 192639 


1.00±0.27 


0.35±0.26 












HD 206267 


43.3±6.5 


53.7±6.3 


12.3±1.4 


1.40±0.34 


10.2±1.2 


14.1±1.4 




HD 207198 


18.4±2.6 


13.1±1.6 


1.32±0.41 


<0.3 


3.3±0.7 


(2.3) a 




HD 210121 


25±8 


31±9 


12±4 


0.31±0.08 


2.8±1.0 


3.8±1.1 




HD 210839 


16.4±3.6 


11.3±1.5 


1.33±0.13 


<0.1 


2.3±0.4 


(1.9) a 




E(./)/k (K) 


0.00 


5.53 


16.60 


33.19 


0.00 


5.53 


16.60 



a Assumes same excitation temperature as for 12 CO. 
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Tabic 7. C 2 D-X (0-0) Band Equivalent Widths 



Line 


A 


ga 


HD 24534 


HD 27778 


HD 147888 


HD 207198 


R(0) 


2313.163 


1.000 


3.7±0.8 


3.7±0.5 


2.6±0.5 


4.8±0.9 


R(2) 


2312.769 


0.600 


9.5±0.8 


8.4±0.6 


3.9±0.5 


12.1±0.9 


P(2) 


2313.745 


0.400 


7.1±0.8 


5.1±0.5 


3.0±0.5 


7.8±0.8 


R(4) 


2312.371 


0.556 


8.8±0.8 


9.2±0.6 


4.7±0.6 


11.5±0.9 


P(4) 


2314.126 


0.444 


8.0±0.8 


6.0±0.5 


3.6±0.5 


10.1±0.9 


R(6) 


2311.967 


0.538 


6.9±0.8 


6.2±0.6 


3.0±0.5 


9.4±0.9 


P(6) 


2314.503 


0.462 


6.4±0.8 


4.9±0.5 


2.8±0.6 


8.7±0.9 


R(8) 


2311.560 


0.529 


4.9±0.7 


4.0±0.6 


2.1±0.5 


5.9±0.8 


P(8) 


2314.874 


0.471 


4.4±0.8 


3.7±0.6 


1.9±0.6 


6.6±0.9 


R(10) 


2311.147 


0.524 


4.1±0.8 


2.6±0.5 


2.2±0.5 


4.3±0.8 


P(10) 


2315.240 


0.476 


3.0±0.9 


2.0±0.5 


1.6±0.6 


4.6±0.8 


R(12) 


2310.730 


0.520 


3.2±0.8 


2.0±0.5 


1.1±0.5 


3.8±0.7 


P(12) 


2315.601 


0.480 


2.5±0.7 


1.8±0.6 


0.7±0.5 


2.4±0.7 


R(14) 


2310.309 


0.517 


3.2±0.9 


1.4±0.5 


<1.5 


2.5±0.8 


P(14) 


2315.956 


0.483 


1.2±0.6 


1.2±0.5 


<1.5 


1.5±0.9 


R(16) 


2309.884 


0.515 


<2.4 


1.1±0.4 


<1.5 


1.7±0.8 


P(16) 


2316.307 


0.485 


<2.4 


1.2±0.5 


<1.5 


1.7±0.8 


R(18) 


2309.455 


0.514 


<2.4 


<1.5 


<1.5 


2.6±0.9 


P(18) 


2316.652 


0.486 


<2.4 


<1.5 


<1.5 


<2.4 


Note. 


— Equivalent widths 


are in mA. 


Uncertainties 


are 1(t; limits 


arc 3a". 



a Rclative /-value for individual line. Adopted band /-value is 0.0545 (Lambert et 
al. 1995). 
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Table 8. C 2 A-X Band Equivalent Widths 



Line A log /A HD 24534 HD 27778 HD 147888 HD 204827 HD 206267 HD 207198 HD 210121 HD 210839 



(2-0) band 



R(0) 


8757.686 


1.089 


1.2±0.3 


1.2±0.3 


1.0±0.3 


16.1±0.5 


3.5±0.4 


2.6±0.5 


3.8±0.4 


2.0±0.4 


R(2) 


8753.949 


0.690 


2.9±0.5 


2.0±0.3 


1.2±0.3 


20.1±0.6 


3.9±0.5 


<5.0 


3.6±0.4 




Q(2) 


8761.194 


0.788 


3.2±0.4 


3.2±0.4 


1.7±0.4 


28.2±0.5 


7.3±0.5 


3.6±0.6 


6.0±0.4 


1.9±0.5 


P(2) 


8766.031 


0.089 




1.2±0.3 




4.7±0.4 


2.2±0.4 


<1.5 






R(4) 


8751.685 


0.611 




(3.1±0.4) 


(0.9±0.3) 


(20.3±0.6) 


(6.3±0.7) 


(3.1±0.8) 


(5.4±0.8) 


(1.6±0.6) 


Q(4) 


8763.751 


0.788 


3.4±0.4 


3.5±0.4 


1.2±0.4 


25.3±0.6 


7.5±0.6 


5.3±0.5 


6.0±0.4 


3.3±0.6 


R(6) 


8750.848 


0.577 




2.7±0.4 


0.7±0.3 


9.7±0.6 


1.9±0.5 


1.5±0.6 


3.3±0.6 


1.9±0.4 


Q(6) 


8767.759 


0.788 


1.6±0.4 


3.4±0.3 


2.6±0.4 


18.6±0.6 


4.4±0.6 


3.4±0.6 


3.6±0.5 


1.1±0.5 


P(6) 


8782.308 


0.373 








5.8±0.5 


1.2±0.3 




1.3±0.5 




Q(8) 


8773.221 


0.788 


(2.1±0.4) 


(3.3±0.4) 


(1.5±0.5) 


(17.5±0.6) 


(4.7±0.6) 


(4.2±0.6) 


(6.9±0.8) 


(1.0±0.4) 


Q(io) 


8780.141 


0.789 




<2.0 


<2.0 




<2.0 








R(12) 


8757.127 


0.536 








1.0±0.3 










Q(12) 


8788.558 


0.789 




<2.0 


<2.0 


3.6±0.5 


0.8±0.5 


0.9±0.5 


0.7±0.3 




(3-0) band 




















R(0) 


7719.329 


0.700 


1.0±0.3 


1.0±0.3 




6.5±0.2 


1.2±0.2 




1.2±0.3 




R(2) 


7716.528 


0.302 


1.4±0.3 






7.4±0.3 


2.1±0.2 




1.8±0.2 




Q(2) 


7722.095 


0.400 


1.9±0.3 












2.9±0.3 




P(2) 


7725.819 


-0.299 








3.0±0.2 










R(4) 


7714.944 


0.223 


0.9±0.3 


0.4±0.2 






1.2±0.2 




1.6±0.3 




Q(4) 


7724.219 


0.400 


1.9±0.3 


1.1±0.3 




9.9±0.3 


2.6±0.3 








P(4) 


7731.663 


-0.076 










0.6±0.2 








R(6) 


7714.575 


0.189 


0.7±0.3 








0.4±0.2 




0.9±0.3 




Q(6) 


7727.557 


0.400 


0.8±0.3 


0.5±0.2 




7.1±0.4 


1.8±0.3 




2.4±0.3 




P(6) 


7738.737 


-0.015 








2.5±0.3 










R(8) 


7715.415 


0.169 


0.7±0.2 






1.4±0.2 


0.4±0.2 




0.4±0.2 




Q(8) 


7732.117 


0.400 


0.4±0.2 








0.5±0.2 




1.4±0.3 




Q(10) 


7737.904 


0.401 








1.7±0.2 











Note. — Equivalent widths are in mA; uncertainties are ler; limits are 3<r. Values for blended lines [(2-0) R(4)/R(8) and 
Q(8)/P(4)] are given in parentheses; (2-0) Q(10) is usually blended with a telluric line; (3-0) Q(2) is sometimes blended 
with a weak diffuse interstellar band. Assumed / — 0.00140 for the (2-0) band and 0.00065 for the (3-0) band. 



Tabic 9. C2 Column Densities 



Star 


Band 




.7 = 




2 






4 






6 






8 






10 




12 




14 


16 


18 




0-6 




Total a 


HD 24534 


F-X(l-O) 





13±0 


.03 





,71±0 


.07 





73±0 


.07 





55±0 


.05 





,37±0 


.05 





,27±0.05 





,14±0 


.05 








2 


.12±0 


.11 


3.07±0.36 




F-X(O-O) 





.17±0 


.02 





.68±0 


.05 





77±0 


.06 





43±0 


.05 





,16±0 


.03 





,05±0.03 





,12±0 


.04 








2 


.05±0. 


.09 


2.97±0.34 




F-X(0-0) b 





18±0 


.02 





.65±0 


.04 





,67±0 


.05 





.42±0 


.04 





12±0 


.03 




<0.05 





.08±0 


03 


0.11±0.03 


0.08±0.03 


0.05±0.02 


1 


.92±0. 


.08 


2.78±0.31 




D-X(O-O) 





15±0 


.02 





,67±0 


.05 





,69±0 


.04 





,55±0 


.04 





,37±0 


.04 





23±0.03 





,22±0 


03 


0.16±0.03 


0.04±0.03 




2 


.06±0 


.08 


2.99±0.33 




A-X(3-0) 





,31±0 


.05 


1 


11±0 


.10 





98±0 


11 





54±0 


.10 





33±0 


.11 


















2 


.94±0 


19 


4.26±0.52 




A-X(2-0) 





,14±0 


.04 





.73±0 


.08 





,77±0 


.09 





,35±0 


.00 





18±0 


.08 


















1 


.99±0. 


.16 


2.89±0.37 




Average 





,16±0 


.02 





,72±0 


.06 





,74±0 


.06 





,51±0 


.05 


(0 


33±0 


.04) 


(0 


22±0.03) 





,18±0 


.04 


0.16±0.03 


0.04±0.03 




2 


.13±0 


.10 


3.09±0.35 


HD 27778 


F-X(O-O) 





15±0 


.04 





,49±0 


.08 





,46±0 


.08 





25±0 


.08 





08±0 


.07 





08±0.07 














1 


.35±0. 


14 


2.00±0.29 




D-X(O-O) 





,15±0 


02 





.57±0 


.04 





.62±0 


.04 





.47±0 


.03 





,32±0 


03 





,19±0.02 





15±0 


.02 


0.11±0.02 


0.09±0.02 


0.06±0.02 


1 


.81±0 


.07 


2.68±0.29 




A-X(3-0) 





27±0 


05 





,45±0 


.12 





.54±0 


.10 





.23±0 


.10 
























1 


.49±0. 


.19 


2.21±0.36 




A-X(2-0) 





.12±0 


.02 





,65±0 


.04 





,68±0 


.04 





73±0 


.04 





,44±0 


.04 





09±0.08 














2 


.18±0. 


07 


3.23±0.34 




Average 





15±0 


.02 





.59±0 


.05 





.62±0 


.05 





.52±0 


.04 


(0 


,34±0 


.04) 


(0 


18±0.03) 





15±0 


.02 


0.11±0.02 


0.09±0.02 


0.06±0.02 


1 


.88±0 


.08 


2.78±0.31 


HD 147888 


D-X(O-O) 





,10±0 


.02 





.33±0 


.02 





.37±0 


.03 





.28±0 


.02 





18±0 


.02 





14±0.02 





.08±0 


.02 


0.05±0.02 






1 


.09±0. 


.04 


1.64±0.18 




A-X(2-0) 





,09±0 


.04 





,37±0 


.07 





,24±0 


.07 





,42±0 


.07 





18±0 


.06 


















1 


.12±0. 


13 


1.69±0.26 




Average 





10±0 


.02 





,33±0 


03 





,36±0 


03 





,30±0 


03 





18±0 


03 





14±0.02 





,08±0 


.02 


0.05±0.02 






1 


.09±0 


05 


1.64±0.18 


C Oph 


D-X(O-O) 





,08±0 


.02 





.29±0 


.02 





,33±0 


.02 





25±0 


.02 





19±0 


.02 





,14±0.02 





11±0 


.01 


0.08±0.02 


0.06±0.02 


0.06±0.01 





.94±0. 


.04 


1.62±0.17 




A-X(3-0) 





11±0 


03 





,30±0 


.05 





,28±0 


06 





,26±0 


.05 





16±0 


.04 





13±0.04 

















.95±0. 


.10 


1.63±0.23 




A-X(2-0) 





,08±0 


.01 





,24±0 


.02 





,34±0 


.05 





,26±0 


.04 





16±0 


.05 





14±0.03 





,10±0 


.02 











.92±0. 


07 


1.58±0.20 




Average 





08±0 


.01 





,26±0 


.02 





.32±0 


.03 





.25±0 


.03 





,18±0 


03 





14±0.02 





,11±0 


.02 


0.08±0.02 


0.06±0.02 


0.06±0.01 





.92±0. 


.05 


1.58±0.17 


HD 204827 


A-X(3-0) 


1 


,87±0 


15 


5 


75±0 


.15 


5 


19±0 


.15 


3 


,46±0 


.15 


1 


,59±0 


.15 





76±0.15 














16.27±C 


|.30 


20.01±1.99 




A-X(2-0) 


1 


64±0 


.15 


5 


.50±0 


.15 


4 


93±0 


.15 


3 


.29±0 


.15 


2 


00±0 


.15 





,79±0.15 





.64±0 


.15 








IE 


>.36±C 


i.30 


18.89±1.88 




Average 


1 


76±0 


.15 


5 


63±0 


.15 


5 


06±0 


.15 


3 


,38±0 


.15 


1 


,80±0 


.15 





78±0.15 





,64±0 


.15 








IE 


>.82±C 


i.30 


19.45±1.93 


HD 206267 


F-X(O-O) 





,52±0 


.05 


1 


55±0 


.10 


1 


30±0 


09 





,46±0 


.00 





13±0 


.07 





11±0.08 














3 


.83±0. 


.17 


4.90±0.54 




A Y ( Q n't 





.41±0 


.05 


1 


30i0 


. 10 


1 


25i0 


. 10 





.70i0 


.00 





.36i0 


.00 


















3 


.66±0 


.17 


4.b8±0.0o 




A-X(2-0) 





40±0 


.03 


1 


.37±0 


.06 


1 


32±0 


.06 





73±0 


.06 





59±0 


.05 









,11±0 


.05 








3 


.82±0 


.11 


4.89±0.52 




Average 





43±0 


.04 


1 


39±0 


.08 


1 


30±0 


.08 





66±0 


.07 


(0 


46±0 


.06) 


(0 


22±0.08) 





11±0 


.05 








3 


.78±0. 


.14 


4.84±0.52 


HD 207198 


D-X(O-O) 





,20±0 


.02 





83±0 


.04 





91±0 


.05 





.75±0 


.04 





,51±0 


.04 





,37±0.04 





.23±0 


.04 


0.17±0.04 


0.13±0.04 


0.13±0.04 


2 


.68±0. 


.08 


4.15±0.42 




A-X(2-0) 





25±0 


.04 





52±0 


.06 





.83±0 


06 





54±0 


.07 





54±0 


06 





41±0.07 





.17±0 


06 








2 


.16±0. 


.12 


3.32±0.37 




Average 





,21±0 


03 





73±0 


.05 





,88±0 


.05 





69±0 


.05 





52±0 


.05 





,38±0.04 





,21±0 


.05 


0.17±0.04 


0.13±0.04 


0.13±0.04 


2 


.50±0. 


09 


3.88±0.40 


HD 210121 


A-X(3-0) 





,34±0 


.08 


1 


.47±0 


.16 


1 


.16±0 


.16 


1 


21±0 


.16 





50±0 


.15 


















4 


.18±0 


29 


6.06±0.75 




A-X(2-0) 





,38±0 


.04 


1 


.24±0 


.07 


1 


37±0 


.07 





90±0 


.07 





,72±0 


.07 





50±0.07 





.18±0 


06 








3 


.88±0. 


13 


5.64±0.61 




Average 





36±0 


.05 


1 


28±0 


.08 


1 


34±0 


.08 





95±0 


.08 





68±0 


.08 





50±0.07 





18±0 


06 








3 


.94±0. 


.15 


5.71±0.63 


HD 210839 


A-X(2-0) 





,21±0 


.05 





40±0 


.10 





52±0 


.00 





35±0 


.08 
























1 


.48±0. 


.16 


3.42±0.51 



Table 9 — Continued 



Star Band J — 2 4 6 8 10 12 14 16 18 0-6 Total a 



Note. — Column densities are in units of 10 13 cm -2 . Uncertainties 
arc l-cr, and include contributions from photon noise and uncertainties 
in continuum placement and fc-value. Average values for J — 8 and 
10 in parentheses use 2x the nominal values from the F-X(O-O) band. 
Values for £ Oph are based on data from Lambert et al. (1995) [D-X(0- 
0)], van Dishocck & Black (1986b) [A-X(3-0)], and Danks & Lambert 
(1983) and Hobbs &c Campbell (1982) [(A-X(2-0)]. 

^Estimated total N(C2) (including unobserved higher J popula- 
tions) obtained by multiplying observed sum for levels 0—6 by ratio of 
theoretical sums for J — 0-20 and J — 0-6, for the best-fit n and Tk 
model. 

b Column densities Nj from fit to combined (3-epoch) spectrum, 
with only lines for J < 10 shifted from predicted wavelengths (see 
Fig. 11); these Nj are not used in the weighted average. 



Tabic 10. Total Column Densities 



Star 


H 2 




12 co 




13 co 


c 2 


12 CO/H 2 


12 CO/ 13 CO 


C2/H2 


HD 


24534 


20.92±0 


04 


16.15±0 


.07 


14.29±0.02 


13.49±0.05 


-4.77±0.08 


1.86±0.07 


-7.43 ± 0.06 


HD 


27778 


20.79±0. 


.06 


16.08±0. 


.03 


14.18±0.05 


13.44±0.05 


-4.71±0.07 


1.90±0.06 


-7.35 ± 0.08 


HD 


73882 


21.11±0 


08 


16.55±0. 


.20 


15.15±0.20 


13.60±0.06 


-4.56±0.25 


1.40±0.35 


-7.51 ± 0.10 


HD 


147888 


20.48±0. 


.08 


15.30±0. 


08 


<13.5 


13.22±0.05 


-5.18±0.12 


>1.8 


-7.26 ± 0.09 


HD 


185418 


20.76±0. 


.05 


14.75±0. 


03 




<12.85 


-6.01±0.06 




< -7.91 


HD 


192639 


20.69±0. 


05 


14.13±0 


13 




<12.85 


-6.56±0.13 




< -7.84 


HD 


206267 


20.86±0. 


04 


16.04±0. 


04 


14.39±0.03 


13.68±0.05 


-4.82±0.05 


1.66±0.05 


-7.18 ± 0.06 


HD 


207198 


20.83±0. 


04 


15.52±0. 


04 


[13.75±0.10] 


13.59±0.05 


-5.31±0.06 


[1.77±0.09] 


-7.23 ± 0.07 


HD 


210121 


20.75±0. 


12 


15.83±0. 


.08 


13.82±0.10 


13.76±0.05 


-4.92±0.15 


2.01±0.13 


-6.99 ± 0.13 


HD 


210839 


20.84±0. 


.04 


15.46±0 


.06 


[13.59±0.08] 


13.53±0.07 


-5.38±0.07 


[1.87±0.09] 


-7.31 ± 0.08 



Note. — Column densities arc log [N (cm 2 )]; uncertainties arc 1-cr; limits arc 3-<r. Values for 13 CO in square braces 
assume same AT /JVi as found for 12 CO. JV(H 2 ) arc from Rachford et al. 2002, except for HD 147888 (Cartlcdge et al. 2004). 
JV(CO) and JV(C 2 ) are from this paper, except for CO values toward HD 24534 (Sheffer et al. 2002a). 



Table 11. Correlation Coefficients and Fits 



Y 


X 


N d a 


r b 


Slopc c 


rms c 


Slope d 


rms d 


C 2 


H 2 


23 


0.678 


1.79 ± 0.35 


0.11 


3.33 ± 0.78 


0.11 


12 CO 


CN 


31 


0.900 


1.50 ± 0.11 


0.16 


1.51 ± 0.11 


0.16 


12 CO 


CH 


44 


0.854 


4.00 ± 0.32 


0.16 


3.81 ± 0.30 


0.17 


12 co 


C 2 


21 


0.450 


3.13 ± 0.78 


0.18 


3.70 ± 0.82 


0.17 


12 co 


c 3 


7 


0.586 


4.00 ± 2.08 


0.11 


4.76 ± 2.72 


0.11 


12 co 


HD 


19 


0.905 


1.82 ± 0.20 


0.35 


1.68 ± 0.17 


0.35 


12 co 


K I 


40 


0.733 


3.23 ± 0.42 


0.21 


3.23 ± 0.31 


0.21 


12 CO/H 2 


CN/CH 


28 


0.815 


1.54 ± 0.19 


0.17 


1.56 ± 0.16 


0.17 


13 CO/H 2 


CN/CH 


11 


0.868 


1.74 ± 0.38 


0.16 


1.65 ± 0.28 


0.16 


C 2 /H 2 


CN/CH 


26 


0.208 






0.05 ± 0.10 


0.17 


12 CO/ 13 CO 


12 CO 


17 


-0.647 


-0.34 ± 0.11 


0.23 


-0.37 ± 0.10 


0.22 


12 CO/ 13 CO 


CN/CH 


12 


-0.549 






-0.71 ± 0.23 


0.22 



a Number of detected values 
b Lincar correlation coefficient 
c Wcighted fit 
d Unweighted fit 



Tabic 12. CO Excitation Temperatures 



Star Name 12 CO 13 CO Ref 







Toi 


Tl2 


T23 


T02/T03 


Toi 


T12 


T23 




HD 22951 


40 Per 


3.0±0.3 














1 


HD 23180 


o Per 


3.6±0.3 














1 


HD 23625 




2.7±0.3 














1 


HD 24534 


X Per 


3.3 








3.7 






2 


HD 147933 


p Oph A 


2.7±0.1 


7.6±0.5 


8.4±0.5 




(7.7±2.5) 


(7.7±2.5) 


(7.7±2.5) 


3 


HD 149757 


C Oph 


3.4±0.4 


4.5±0.2 


5.9±0.5 




4.2±0.4 


4.4±0.2 


5.9±0.5 


4 


HD 148184 


X Oph 


3.0±0.3 


5.2±0.6 


7.5±3.7 




(4.4±2.9) 


(4.4±2.9) 


(4.4±2.9) 


3 


HD 203374A 




3.2 


3.5 


8.7 










5 


HD 24534 


X Per 


4.6±1.7 


5.3±0.7 


5.6±1.0 


5.1/5.4 


3.3±0.3 


4.0±0.2 




6 


HD 27778 


62 Tau 


5.2±0.7 


5.6±0.4 


6.3±0.6 


5.5/5.8 


4.3±0.8 






6 


HD 147888 


p Oph D 


11.1: 


7.9±2.3 


8.9±1.5 


8.6/8.7 








6 


HD 185418 




3.3±0.3 


4. Oil. 2 




3.7/... 








6 


HD 192639 




2.6±1.2 






.../... 








6 


HD 206267 




6.3±1.4 


5.6±0.5 


6.6±0.7 


5.8/6.1 


7.1±0.4 






6 


HD 207198 




3.8±0.5 


3.9±0.5 


<9.1 


3.9/... 








6 


HD 210121 




6.3±3.2 


7.6±2.4 


4.2±0.5 


7.2/5.3 


7.0±4.2 






6 


HD 210839 


A Ccp 


3.8±0.7 


4.2±0.3 


<5.7 


4.1/... 








6 


Diffuse (50 K) 




6.7 


4.9 


6.9 


5.3/5.9 


7.3 


5.1 


7.0 


7 


Diffuse (TE) 




5.4 


4.4 


6.2 


4.6/5.3 


5.7 


4.6 


6.4 


7 


Translucent (25 K) 


17.5 


18.4 


15.2 


18.2/16.6 


9.3 


5.8 


6.5 


7 


Translucent (r 


=0.5) 


12.2 


6.0 


7.2 




8.8 


5.4 


5.2 


7 


Translucent (TE) 


9.5 


8.7 


6.8 


8.9/7.7 


5.4 


4.3 


4.7 


7 



Note. — First set of sight lines is from other references; second set is from this paper; last set gives theoretical values 
from isothermal and thermal equilibrium (TE) models of Warin ct al. 1996. CO values in parentheses are from fits to all 
observed levels. 

References. — 1 = Wannicr ct al. 1997; 2 = Shcffcr ct al. 2002b [A-X (8-0)]; 3 = Fcderman et al. 2003; 4 = Lambert ct 
al. 1994; 5 = Shcffer ct al. 2003; 6 = this paper; 7 = Warin ct al. 1996 (theory) 



Tabic 13. C2 Excitation: Temperature and Density 



Star 


Name 


H 2 a 






Ca- 
























T m 


Tb 2 


To 4 


lk 


n b 




Rcfs d 


HD 23180 


Per 


48 




48±37 


(50) 


(200;f °) 


(265l 3 6 f) 


1 


HD 24398 


C Per 


58 




84±53 


(50) 


(150tf °) 


(215l 2 7 8 5 ) 


2 


HD 62542 




43±11 




43±16 


30±10 


3B0±*Sg 


500t 5 ™ 


3 


HD 80077 






40±26 


33± 8 


20± 5 


200l 2 ™ 




4 


HD 110432 




68± 5 


13± 4 


27± 6 




100 + 2 r 


140lf 5 


5 


HD 147084 


Sco 






57±29 


35±10 


150+™ 




5 


HD 147889 






52±22 


116±28 


(70) 


(350lt°° ) 




6 


HD 148184 


X Oph 


46± 6 


39±16 


57±12 


50±15 


150+1™ 


185lf 5 


6 


HD 154368 




51± 8 


27± 8 


38± 5 


20± 5 


1501™ 


2401™ 


6 


HD 169454 






24± 6 


37± 5 


20± 5 


200+™ 


[265lf 5 ] 


4 


HD 172028 






64±40 


48± 7 


40±10 


2501™ 




4 


HD 179406 


20 Aql 
Cyg OB2 5 


59± 7 




57±11 
49±20 


(55) 
35±15 


>250 
1501™ 


[>335] 


7 

8 




Cyg OB2 12 




60±24 


61± 9 


40±15 


2OOI50 




8 


HD 200775 






20± 9 


28± 7 


15± 5 


200 + 2 ° 


[265lf 5 ] 


9 


HD 207538 




73± 8 




62±14 


35±15 


150lf 5 


190lf 


10 



HD 


24534 


X Per 


57± 4 




72±13 


45±10 


2001™ 


' 3zo _S0 


11 


HD 


27778 


62 Tau 


55± 7 


65±43 


67±12 


50±10 


200+2^ 


280+3? 


11 


HD 


73882 




51± 6 


23± 5 


25± 4 


20± 5 


3501™° 


520l^ 5 5 


3 


HD 


147888 


p Oph D 


47±12 


38±20 


66±13 


55± 5 


200+25 


2151™ 


11 


HD 


149757 


C Oph 


54± 3 


36±12 


74±15 


50± 5 


150lf 5 


215±|» 


12 


HD 


185418 




110±15 












11 


HD 


192639 




90±10 












11 


HD 


204827 






35± 7 


51± 3 


40± 5 


3001™ 


[400lf 5 ] 


11 


HD 


206267 




65± 5 


36± 9 


51± 5 


35± 5 


2501™ 


0-1 r + 125 


11 


HD 


207198 




66± 5 


43±19 


81±14 


60±10 


200lf 


2451™ 


11 


HD 


210121 




51±11 


46±21 


64± 8 


45± 5 


2001™ 


3151™ 


11 


HD 


210839 


A Cep 


72± 6 


16± 6 


47±12 


(30) 


(loot™) 


(115l 5 3?) 


11 



Note. — C2 7k an d inferred n and nn values for upper set of sight lines are based on equivalent width 
data from various references; values for most of the sight lines in the lower set are based on data given in this 
paper. Values in parentheses are estimates based on the assumption that Tk < T01 from H2. All values for nu 
assume a collisional de-excitation cross-section a — 2 X 10~ 16 cm 2 and the average Galactic radiation field. 

References. — 1 = Hobbs 1981; 2 = Chaffee ct al. 1980; 3 = Grcdcl ct al. 1993; 4 = Grcdcl 1999; 5 = van 
Dishoeck & Black 1989; 6 = van Dishocck & de Zecuw 1984; 7 = Fcdcrman et al. 1994a: 8 = Grcdcl ct al. 
2001; 9 = Fcdcrman et al. 1997b; 10 = Galazutdinov ct al. 2006; 11 = this paper; 12 = Hobbs & Campbell 
1982, Danks & Lambert 1983, van Dishocck Sz Black 1986b, Lambert ct al. 1995 

a Savagc ct al. 1977; Rachford ct al. 2002 

b Density of collision partners (H I + H2). 

c Dcnsity of hydrogen nuclei (H I + 2H 2 ). For HD 169454, HD 179406, HD 200775, and HD 204827, /(H 2 ) 
— 0.5 is assumed. 

d Reference for C2 equivalent widths (temperatures and densities are those calculated in this paper). 
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Table 14. Component Structures for K I, CH, and CN 



Star 


V 

(km s" 1 ) 


b 

(km s" 1 ) 


N(K I) 
(10 10 cm" 2 ) 


JV(CH) 
(10 12 cm- 2 ) 


V 

(km s- 1 ) 


b 

(km s- 1 ) 


AT(CN) 
(10 12 cm- 2 ) 


HD 24534 


11.2 


0.8 


9.1 


5.0 










12.4 


0.6 


46.7 


6.0 


12.3 


0.8 


0.7 




14.5 


1.0 


63.5 


25.9 


14.7 


1.0 


6.6 




16.5 


0.8 


2.8 










HD 27778 


12.6 


0.6 


2.2 


0.3 










14.1 


0.4 


13.6 


4.0 










15.2 


0.7 


19.2 


12.8 


15.1 


1.0 


8.9 




17.4 


0.8 


14.2 


9.4 


17.3 


0.9 


5.1 




18.6 


0.4 


17.5 


2.0 










20.0 


0.7 


1.8 


1.3 








HD 73882 


16.7 


0.7 


0.8 












19.4 


0.7 


4.5 


8.9 










21.0 


0.7 


8.3 


16.4 


21.3 


0.8 


38.1 




22.6 


0.7 


4.4 


8.7 










24.4 


1.0 


1.2 












27.6 


0.9 


3.9 


0.8 








HD 147888 


-12.6 


0.5 


1.2 












- 9.7 


0.8 


24.6 


5.9 










-8.5 


0.8 


46.0 


10.1 


-8.5 


0.8 


2.1 




-7.3 


0.9 


23.7 


6.7 








HD 185418 


-10.9 


1.1 


54.0 


8.6 










-8.2 


0.5 


10.0 


1.1 










-6.5 


0.5 


12.0 


3.2 








HD 192639 


-18.5 


0.8 


10.0 


2.3 










-16.5 


1.1 


25.0 


4.7 










-14.0 


1.0 


6.0 


6.1 










-10.7 


1.2 


20.0 


6.8 










-8.4 


1.1 


14.0 


7.8 








HD 206267 


-21.7 


1.0 


1.1 


1.2 










-18.7 


0.7 


29.3 


4.5 


-17.9 


0.7 


0.4 




-17.1 


0.8 


49.1 


5.6 


-16.6 


0.9 


6.5 




-15.7 


0.8 


35.4 


2.6 










-13.6 


0.6 


26.8 


7.3 


-13.8 


0.7 


1.0 




-12.0 


1.0 


26.0 


4.9 










-9.7 


0.8 


31.1 


3.3 


-10.7 


0.7 


0.2 




-7.3 


1.0 


2.7 


0.9 








HD 207198 


-21.4 


1.0 


4.8 


0.6 










-19.4 


0.6 


29.8 


3.7 


-20.1 


0.4 


1.3 




-17.1 


1.0 


19.9 


3.5 










-15.3 


1.2 


62.3 


10.7 


-16.2 


1.2 


2.5 




-12.9 


0.9 


70.4 


8.5 


-13.7 


1.0 


0.6 




-10.2 


0.8 


19.1 


3.8 










-8.3 


0.7 


46.1 


3.6 










-7.0 


1.0 


4.8 


1.2 








HD 210121 


-16.5 


1.1 


21.0 


5.4 










-14.2 


1.2 


94.0 


24.2 


-14.3 


1.1 


16.0 



Table 14 — Continued 



Star 



(km s 



v 



-1 



b N(K I) JV(CH) v b N(CN) 

) (kms- 1 ) (10 10 cm- 2 ) (10 12 cm- 2 ) (kms" 1 ) (kms" 1 ) (10 12 cm" 2 ) 



HD 210839 



-34 
-32. 
-31 
-17 
-15 
-13 
-11. 
-9. 
-8. 
-6 
-5. 
-2 
-1 



.6 
.9 
.1 
.3 
.4 
.4 

6 
.4 
.6 
.2 
.7 
.5 







1. 
0. 

{) 
1. 

0. 

{) 
1. 


0. 

{) 



.8 
6 
.5 

6 
.7 
1 
6 
.8 

6 
.5 
.8 



11.2 
1.3 
1.9 
3.9 
41.5 
39.4 
12.3 
3.9 
1.8 
1.7 
7.6 
9.7 
0.7 



I. 7 

II. 2 
6.7 
1.2 



-15.3 
-13.5 



0.5 
0.5 



3.1 
0.7 



Note. — CH fitted with K I component structure; CN fitted independently. K I structures for HD 206267, 
HD 207198, and HD 210839 are from Welty & Hobbs 2001. CN structure for HD 73882 is from Grcdel et al. 1991 
(FWHM ~ 5 km s™ 1 ). CN structure for HD 147888 is from Pan ct al. 2004. CN structure for HD 210121 is from 
Welty & Fowler 1992 (FWHM ~ 3 km s _1 ). CH structure for HD 210839 is based on new fit to profile from Crane 
ct al. 1995, assumed aligned with strongest components in K I and CN. 
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Table 15. 


c 2 


Band Structures 








Line 


E(./)/k 




F-X bands 




D-X bands 




A-X bands 








A(l-O) 


A(0-0) 


S a 


A(0-0) 


S a 


A(3-0) 


A(2-0) 


S a 




(K) 


(A) 


(A) 




(A) 




(A) 


(A) 






00 


1314 170 




1 000 


9313 163 


1 000 


7719 329 


O 1 O 1 . uou 


1 000 




15 64 


1314 079 


1 "341 ^39 


400 


9^1 9 7fiQ 


600 


771 ^9& 


S7^s 040 


400 






1314 946 


J- 1.IUO 


500 






7799 DQ^ 


«7fi1 104 


500 






1314 35X 




100 


9^1 ^ 7A^ 


400 


779^ Q 


1 uu.uai 


100 


RC41 


52.11 


1314.015 


1341.458 


0.333 


2312.371 


0.556 


7714.944 


8751 685 


0.333 






1314 9Q4 


1341 759 

1041. 1 O £ 


500 






7724 219 


£763 751 


500 






1314 517 


1 "iAI QSS 

±o*±± .yoo 


167 


2314.126 


0.444 


77^1 fifi^ 


oil O.tOU 


167 




109 43 


1313 977 


1341 411 


308 


2311.967 


0.538 


7714 575 


8750 S4S 


308 


0(6) 




1314.368 


1341.822 


0.500 






7727.557 


8767.759 


0.500 


P(6) 




1314.703 


1342.176 


0.192 


2314.503 


0.462 


7738.737 


8782.308 


0.192 


R(8) 


187.57 


1313.967 


1341.388 


0.294 


2311.560 


0.529 


7715.415 


8751.486 


0.294 


Q(8) 




1314.469 


1341.917 


0.500 






7732.117 


8773.221 


0.500 


P(8) 




1314.915 


1342.388 


0.206 


2314.874 


0.471 


7747.037 


8792.648 


0.206 


R(10) 


286.53 


1313.983 


1341.390 


0.286 


2311.147 


0.524 


7717.470 


8753.578 


0.286 


Q(io) 




1314.597 


1342.038 


0.500 






7737.904 


8780.141 


0.500 


P(10) 




1315.155 


1342.626 


0.214 


2315.240 


0.476 


7756.583 


8804.499 


0.214 


R(12) 


406.27 


1314.027 


1341.418 


0.280 


2310.730 


0.520 


7720.748 


8757.127 


0.280 


Q(12) 




1314.752 


1342.183 


0.500 






7744.900 


8788.558 


0.500 


P(12) 




1315.422 


1342.889 


0.220 


2315.601 


0.480 


7767.371 


8817.826 


0.220 


R(14) 


546.78 


1314.098 


1341.471 


0.276 


2310.309 


0.517 




8762.145 


0.276 


Q(14) 




1314.934 


1342.353 


0.500 








8798.459 


0.500 


P(14) 




1315.715 


1343.178 


0.224 


2315.956 


0.483 




8832.679 


0.224 


R(16) 


708.05 


1314.196 


1341.550 


0.273 


2309.884 


0.515 




8768.627 


0.273 


Q(16) 




1315.143 


1342.549 


0.500 








8809.841 


0.500 


P(16) 




1316.036 


1343.491 


0.227 


2316.307 


0.485 




8849.071 


0.227 


R(18) 


890.02 


1314.321 


1341.653 


0.270 


2309.455 


0.514 








Q(18) 




1315.379 


1342.770 


0.500 












P(18) 




1316.384 


1343.830 


0.230 


2316.652 


0.486 








R(20) 


1092.68 


1314.473 


1341.781 


0.268 


2309.022 


0.512 








Q(20) 




1315.642 


1343.015 


0.500 












P(20) 




1316.758 


1344.193 


0.232 


2316.991 


0.488 








Adopted / 


0.06 


0.10 




0.0545 




0.00065 


0.00140 





Note. — Wavelengths arc adapted from Herzberg ct al. 1969 [F-X bands], Sorkhabi ct al. 1998 [D-X(O-O)], 
and van Dishocck & dc Zccuw 1984 [A-X(2-0)]. 

a Rclativc /-value for individual line. Values for higher-J A-X bands differ slightly from those tabulated by 
van Dishocck &c dc Zccuw 1984. 



- 76 - 



Table 16. Column Densities (H 2 , CH, C 2 , CN, CO) 



HD 




Name 


H 2 a 




CH b 




C 2 C 




CN d 




12 CO c 




13 CO f 




HD 


21278 




















12.87 ± 0.12 


1 






HD 


22951 


40 Per 


20.46 ± 0.18 


1 


13.08 ± 0.04 


1 


12.04 ± 0.12 


1 


11.85 


1 


14.22 


2 






HD 


23180 


o Per 


20.60 ± 0.15 


1 


13.33 ± 0.03 


2 


13.27 ± 0.12 


(2) 


12.17 ± 0.01 


2 


14.93 


2 






HD 


23408 


20 Tau 


19.75 ± 0.26 


2 


< 11.30 


2 






< 11.82 


3 


13.73 ± 0.05 


1 






HD 


23625 








13.08 ± 0.15 


1 


12.91 ± 0.05 


1 


11.85 


1 


14.26 


2 






HD 


23630 


T] Tau 


19.54 ± 0.18 


2 


< 12.34 


3 










< 12.34 


1 






HD 


24398 


C Per 


20.67 ± 0.18 


1 


13.34 ± 0.03 


2 


13.22 ± 0.07 


(3) 


12.51 ± 0.01 


4 


14.86 ± 0.17 


1 






HD 


24534 


X Per 


20.92 ± 0.04 


3 


13.57 ± 0.02 


4 


13.49 ± 0.05 


4 


12.87 ± 0.01 


5 


16.15 ± 0.07 


3 


14.29 ± 0.02 


1 


HD 


24760 


e Per 


19.52 ± 0.26 


1 


< 12.30 


3 










11.93 ± 0.08 


1 






HD 


24912 


£ Per 


20.53 ± 0.15 


1 


13.12 ± 0.09 


2 


12.90 ± 0.13 


5 


11.44 ± 0.16 


6 


13.49 ± 0.04 


3 






HD 


27778 


62 Tau 


20.79 ± 0.06 


3 


13.48 ± 0.04 


4 


13.44 ± 0.05 


4 


13.19 ± 0.02 


4 


16.08 ± 0.03 


4 


14.18 ± 0.05 


2 


HD 


28497 




14.82 


1 


< 12.08 


5 










< 12.67 


1 






HD 


30614 


a Cam 


20.34 ± 0.15 


1 


12.84 ± 0.10 


2 


< 12.75 


6 


< 11.82 


3 


14.49 ± 0.04 


1 






HD 


35149 


23 Ori 


18.30 ± 0.10 


4 


12.70 ± 0.03 


4 


< 12.70 


6 


< 11.73 


5 


< 13.00 


5 






HD 


36486 


S Ori 


14.74 ± 0.05 


5 














< 12.04 


1 






HD 


36861 


A Ori 


19.12 ± 0.20 


1 














12.59 ± 0.06 


1 






HD 


37043 


l Ori 


14.69 ± 0.20 


1 














12.43 ± 0.17 


1 






HD 


37128 


e Ori 


16.28 ± 0.20 


5 














< 12.30 


1 






HD 


37202 


C Tau 


< 17.67 


1 














< 12.32 


1 






HD 


37742 


C Ori 


15.86 


5 














< 12.72 


1 






HD 


37903 




20.92 ± 0.06 


3 


12.93 ± 0.05 


6 


< 12.85 


6 


11.85 ± 0.09 


7 


14.00 ± 0.04 


6 






HD 


40111 


139 Tau 


19.73 ± 0.20 


1 














< 12.32 


1 






HD 


47839 


S Mon A 


15.54 


1 


< 12.26 


7 


< 12.85 


6 


< 11.70 


3 


< 12.81 


1 






HD 


53138 


o 2 CMa 


















< 12.97 


1 






HD 


57060 


29 CMa 


15.78 ± 0.18 


1 


< 12.38 


8 










< 12.67 


1 






HD 


57061 


T CMa 


15.48 ± 0.18 


1 


< 12.26 


8 










< 12.61 


1 






HD 


66811 


C Pup 


14.45 


1 


< 11.78 


8 










< 12.11 


1 






HD 


68273 


7 2 Vol 


14.23 ± 0.20 


1 














< 12.46 


1 






HD 


73882 




21.11 ± 0.08 


3 


13.54 ± 0.04 


4 


13.60 ± 0.06 


(7) 


13.58 ± 0.05 


8 


16.55 ± 0.20 


4 


15.15 ± 0.20 


2 


HD 


87901 


a Leo 


















< 12.28 


1 






HD 


91316 


p Leo 


15.61 


1 


< 12.08 


7 


< 12.63 


6 


< 11.52 


3 


< 12.78 


1 






HD 


93030 


6 Car 


< 17.65 


1 


< 11.78 


8 










< 11.95 


1 






HD 


110432 




20.64 ± 0.04 


3 


13.19 ± 0.07 


9 


13.38 ± 0.07 


(8) 


12.20 ± 0.08 


8 


14.60 ± 0.40 


7 






HD 


112244 




20.14 ± 0.11 


1 


12.70 ± 0.08 


3 






< 11.52 


6 


13.14 ± 0.07 


1 






HD 


135742 


Lib 


















< 12.71 


1 






HD 


141637 


1 Sco 


19.23 ± 0.18 


1 


< 12.26 


4 










< 12.41 


1 






HD 


143018 


7T SCO 


19.32 ± 0.20 


1 


< 11.18 


4 


< 12.23 


5 






11.92 ± 0.11 


1 






HD 


143275 


8 Sco 


19.41 ± 0.20 


1 


12.34 ± 0.04 


2 


< 12.85 


6 


< 11.12 


6 


12.49 ± 0.13 


1 






HD 


144217 


/3 1 Sco 


19.83 ± 0.06 


1 


12.41 ± 0.04 


2 


< 12.23 


5 


< 11.12 


6 


13.63 ± 0.05 


1 






HD 


144470 


a; 1 Sco 


20.05 ± 0.11 


1 


12.51 ± 0.05 


3 


< 12.23 


5 






12.95 ± 0.04 


1 






HD 


145502 


v Sco 


19.89 ± 0.15 


1 


12.77 ± 0.03 


3 










12.76 ± 0.04 


1 






HD 


147165 


a Sco 


19.79 ± 0.15 


1 


12.52 ± 0.07 


2 






< 11.12 


6 


12.48 ± 0.03 


1 






HD 


147888 


p OphD 


20.48 ± 0.08 


6 


13.30 ± 0.02 


4 


13.22 ± 0.05 


4 


12.35 ± 0.02 


9 


15.30 ± 0.08 


4 


< 13.50 


2 


HD 


147933 


p OphA 


20.57 ± 0.15 


1 


13.33 


2 


13.54 ± 0.24 


6 


12.32 ± 0.04 


6 


15.28 ± 0.06 


8 


13.18 ± 0.10 


3 


HD 


148184 


X Oph 


20.63 ± 0.18 


1 


13.48 ± 0.02 


2 


13.46 ± 0.05 


0) 


12.09 ± 0.04 


6 


14.58 ± 0.12 


8 


12.51 ± 0.19 


3 


HD 


149038 


I_l Nor 


20.45 ± 0.15 


1 


13.02 ± 0.04 


8 






< 11.30 


6 


13.91 ± 0.05 


1 






HD 


149757 


C Oph 


20.64 ± 0.08 


1 


13.40 ± 0.02 


2 


13.20 ± 0.05 


(10) 


12.47 ± 0.01 


4 


15.40 ± 0.03 


9 


13.18 ± 0.03 


4 


HD 


157246 


7 Ara 


19.23 ± 0.23 


1 


< 11.78 


8 










12.26 ± 0.04 


1 






HD 


164353 


67 Oph 


20.26 ± 0.28 


1 


12.60 ± 0.07 


2 


< 12.70 


6 


< 11.22 


6 


13.04 ± 0.08 


1 






HD 


175191 


a Sgr 


< 14.00 


1 














< 12.26 


1 






HD 


179406 


20 Aql 


20.73 ± 0.07 


3 


13.27 ± 0.03 


2 


13.71 ± 0.05 


(11) 


12.55 ± 0.01 


7 


15.64 ± 0.01 


6 


14.28 ± 0.02 


5 


HD 


185418 




20.76 ± 0.05 


3 


13.12 ± 0.09 


4 


< 12.85 


6 


11.60 ± 0.11 


9 


14.75 ± 0.03 


4 






HD 


192639 




20.69 ± 0.05 


3 


13.45 ± 0.07 


4 


< 12.85 


6 


11.70 ± 0.11 


9 


14.13 ± 0.13 


4 






HD 


200120 


59 Cyg 


19.30 ± 0.18 


1 














12.18 ± 0.12 


1 






HD 


200775 








13.51 ± 0.03 


6 


12.99 ± 0.08 (12) 


13.08 ± 0.11 


7 


17.30 ± 0.02 


6 


16.04 ± 0.04 


5 
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Table 16 — Continued 



HD 


Name 


H 2 a 




CH b 




C 2 ° 




CN d 




12 co° 




13 co f 




HD 203374A 




20.70 ± 0.01 


7 


13.35 ± 0.04 


10 






12.32 


10 


15.41 ± 0.02 


10 






HD 206165 


9 Ccp 


20.78 ± 0.02 


7 


13.40 ± 0.04 


4 


13.56 ± 0.08 


6 


12.41 ± 0.05 


5 


15.18 ± 0.01 


10 






HD 206267A 




20.86 ± 0.04 


3 


13.48 ± 0.03 


4 


13.68 ± 0.05 


4 


12.89 ± 0.01 


5 


16.04 ± 0.04 


4 


14.39 ± 0.03 


2 


HD 206773 




20.47 ± 0.01 


7 


13.05 ± 0.04 


10 






11.65 


10 


14.20 ± 0.05 


10 






HD 207198 




20.83 ± 0.04 


3 


13.56 ± 0.02 


4 


13.59 ± 0.05 


4 


12.65 ± 0.03 


5 


15.52 ± 0.04 


4 


[13.75 ± 0.10] 


2 


HD 207308 




20.86 ± 0.01 


7 


13.51 ± 0.05 


10 






13.04 


10 


15.93 ± 0.05 


10 






HD 207538 




20.91 ± 0.05 


3 


13.63 ± 0.02 


4 


13.92 ± 0.05 (13) 


12.72 


10 


15.40 ± 0.04 


10 






HD 208266 








13.53 ± 0.04 


10 










16.05 ± 0.04 


10 






HD 208440 




20.34 ± 0.01 


7 


13.06 ± 0.04 


10 






< 11.95 


10 


14.20 ± 0.03 


10 






HD 209339 




20.25 ± 0.01 


7 


12.90 ± 0.05 


10 






< 11.95 


10 


13.95 ± 0.05 


10 






HD 209952 


a Gru 


















< 12.23 


1 






HD 210121 




20.75 ± 0.12 


3 


13.47 ± 0.03 


4 


13.76 ± 0.05 


4 


13.09 ± 0.07 


8 


15.83 ± 0.08 


4 


13.82 ± 0.10 


2 


HD 210839 


A Cop 


20.84 ± 0.04 


3 


13.32 ± 0.05 


2 


13.53 ± 0.07 


4 


12.58 ± 0.02 


5 


15.46 ± 0.06 


4 


[13.59 ± 0.08] 


2 


HD 216532 




21.10 ± 0.01 


7 


13.31 ± 0.04 


10 






12.49 


10 


15.15 ± 0.03 


10 






HD 216898 




21.05 ± 0.01 


7 


13.45 ± 0.05 


10 






11.90 


10 


15.04 ± 0.03 


10 






HD 217035A 




20.95 ± 0.01 


7 


13.23 ± 0.05 


10 






< 11.95 


10 


14.57 ± 0.03 


10 






HD 217312 




20.80 ± 0.01 


7 


13.41 ± 0.05 


10 






12.18 


10 


14.38 ± 0.05 


10 






HD 217675 


o And 


19.67 ± 0.18 


2 


< 12.11 


3 










12.75 ± 0.09 


1 






HD 218376 


1 Cas 


20.15 ± 0.18 


1 


12.90 ± 0.22 


3 


< 12.85 


6 


< 11.82 


3 


14.80 ± 0.09 


1 







a H 2 : 1 = Savage ot al. 1977; 2 = Jenkins et al. 1986; 3 = Rachford et al. 2002 or Rachford, priv. comm.; 4 = Welty et al. 
1999; 5 = Jenkins & Peimbert 1996 or Jenkins et al. 2000; 6 = Cartledge et al. 2004; 7 = Shcffer et al. 2003 or Pan et al. 
2005. 

b CH: 1 = Wannicr et al. 1999; 2 = Crane et al. 1995; 3 = Federman 1982; 4 = Welty et al. 1999 or Welty, Snow, & Morton, 
in prep.; 5 = Penprase 1993; 6 = Knauth et al. 2001; 7 = Thorburn et al. 2003; 8 = Danks et al. 1984; 9 = Crawford 1995; 
10 = Pan et al. 2004. 

C C 2 : 1 = Wannier et al. 1999: 2 = Hobbs 1981; 3 = Chaffee et al. 1980; 4 = this paper: 5 = Lambert et al. 1995; 6 = 
Thorburn et al. 2003; 7 = Grcdel et al. 1993; 8 = van Dishocck & Black 1989; 9 = van Dishoeck fc de Zeeuw 1984: 10 = 
Hobbs & Campbell 1982, Danks & Lambert 1983, van Dishoeck & Black 1986b, Lambert et al. 1995: 11 = Federman et al. 
1994a; 12 — Federman et al. 1997b; 13 — Galazutdinov et al. 2006. Values for references in parentheses are estimated in this 
paper from equivalent widths taken from those references. 

d CN: 1 = Wannier et al. 1999; 2 = Roth & Meyer 1995; 3 = Thorburn ct al. 2003; 4 = Roth & Meyer 1995; 5 = Welty, 
Snow, & Morton, in prep, or Welty et al. 1999; 6 = Federman 1984; 7 = Knauth et al. 2001; 8 = Gredel et al. 1991; 9 = this 
paper: 10 = Pan et al. 2004. 

ol2 CO: 1 = Crcnny & Federman 2004 or Federman et al. 1980; 2 = Wannier et al. 1999; 3 = Sheffer ct al. 2002a, 2002b; 
4 = this paper; 5 = Welty et al. 1999; 6 = Knauth et al. 2001; 7 = Rachford ct al. 2001; 8 = Federman et al 2003; 9 = 
Lambert ct al. 1994; 10 = Sheffer ct al. 2003 or Pan et al. 2005. 

fl3 CO: 1 = Sheffer ct al. 2002a, 2002b: 2 = this paper; 3 = Federman et al. 2003; 4 = Lambert ct al. 1994; 5 = Knauth 
ct al. 2001. Values in square braces assume that Nq/Ni is the same for 13 CO as for 12 CO. 
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Table 17. Derived Quantities 



HD 




Name 


E(B - 


vy A v a 


iV(H tot 


) b log[/(H 2 )] b T i(H 2 ) 


b T ox (C 2 ) c CN/CH 


( 12 CO/H 2 ) 


12 CO/ 13 CO 


HD 


22951 


40 Per 


0.27 





66 


21.23 


-0.47 


63±10 






-1 


23 


-6.24 




HD 


23180 


o Per 


0.31 





84 


21.20 


-0.30 


48 


.../48/... 




-1 


16 


-5.67 




HD 


24398 


C Per 


0.31 





83 


21.20 


-0.23 


58 


.../84/... 




-0 


83 


-5.81 




HD 


24534 


X Per 


0.59 


2 


05 


21.34 


-0.12 


57± 4 


.../72/45 




-0 


70 


-4.77 


1.86 


HD 


24912 


£ Per 


0.33 


1 


03 


21.29 


-0.46 


61 






-1 


68 


-7.04 




HD 


27778 


62 Tau 


0.37 


1 


01 


21.34 


-0.25 


55± 7 


.../67/50 




-0 


29 


-4.71 


1.90 


HD 


30614 


a Cam 


0.3 





88 


21.09 


-0.45 


80 




< 


-1 


02 


-5.85 




HD 


35149 


23 Ori 


0.11 






20.74 


-2.14 


100 




< 


-0 


97 


< -5.30 




HD 


37903 




0.35 


1 


28 


21.50 


-0.28 


68± 7 






-1 


08 


-6.92 




HD 


73882 




0.72 


2 


36 


21.59 


-0.18 


51± 6 


23/25/20 







04 


-4.56 


1.40 


HD 


110432 




0.40 


1 


32 


21.20 


-0.26 


68± 5 


13/27/... 




-0 


99 


-6.04 




HD 


112244 




0.29 






21.17 


-0.73 






< 


-1 


18 


-6.95 




HD 


143275 


S Sco 


0.17 






21.17 


-1.46 


56 




< 


-1 


22 


-6.92 




HD 


144217 


/3 1 Sco 


0.19 





65 


21.14 


-1.01 


89 




< 


-1 


29 


-6.20 




HD 


147165 


a Sco 


0.41 


1 


22 


21.36 


-1.27 


65 




< 


-1 


40 


-7.31 




HD 


147888 


p Oph D 


0.47 


1 


91 


(21.62) 


-0.84 


47±12 


38/66/55 




-0 


95 


-5.18 


> 1.60 


HD 


147933 


p Oph A 


0.48 


2 


68 


21.63 


-0.76 


47± 5 






-1 


01 


-5.29 


2.10 


HD 


148184 


X Oph 


0.52 


2 


34 


21.36 


-0.43 


46± 6 


39/57/50 




— 1 


39 


-6.05 


2.07 


HD 


149038 


H Nor 


0.37 






21.19 


-0.44 






< 


-1 


72 


-6.54 




HD 


149757 


C Oph 


0.32 


1 


06 


21.15 


-0.21 


54± 3 


36/74/50 




-0 


93 


-5.24 


2.22 


HD 


164353 


67 Oph 


0.11 





33 


21.14 


-0.58 


74 




< 


-1 


38 


-7.22 




HD 


179406 


20 Aql 


0.33 





94 






59± 7 


.../57/... 




-0 


72 


-5.09 


1.36 


HD 


185418 




0.51 


2 


03 


21.39 


-0.34 


110±15 






-1 


52 


-6.01 




HD 


192639 




0.66 


1 


87 


21.49 


-0.45 


90±10 






-1 


75 


-6.56 




HD 


200775 




0.63 


3 


64 








20/28/... 




-0 


43 




1.26 


HD 


203374A 




0.60 










76± 1 






-1 


01 


-5.29 




HD 


206165 


9 Ccp 


0.47 










68± 4 






-0 


92 


-5.60 




HD 


206267A 




0.53 


1 


37 


21.54 


-0.38 


65± 5 


36/51/35 




-0 


59 


-4.82 


1.66 


HD 


206773 




0.50 










94± 1 






-1 


39 


-6.27 




HD 


207198 




0.62 


1 


36 


21.55 


-0.42 


66± 5 


43/81/60 




-0 


91 


-5.31 


[1.77] 


HD 


207308 




0.50 










57± 1 






-0 


47 


-4.93 




HD 


207538 




0.64 






21.59 




73± 8 


.../62/35 




-0 


89 


-5.51 




HD 


208440 




0.34 










75± 2 




< 


-1 


11 


-6.14 




HD 


209339 




0.37 










90± 1 




< 


-0 


93 


-6.30 




HD 


210121 




0.38 





83 


21.19 


-0.14 


51±11 


46/64/45 




-0 


38 


-4.92 


2.01 


HD 


210839 


A Ccp 


0.57 


1 


57 


21.45 


-0.31 


72± 6 


16/47/... 




-0 


74 


-5.38 


[1.87] 


HD 


216532 




0.86 










83± 8 






-0 


81 


-5.95 




HD 


216898 




0.85 






21.05 




86± 6 






-1 


55 


-6.01 




HD 


217035A 




0.76 










75± 1 




< 


-1 


28 


-6.38 




HD 


217312 




0.66 










76± 3 






-1 


23 


-6.42 




HD 


218376 


1 Cas 


0.25 





54 


21.07 


-0.62 


84 




< 


-1 


08 


-5.35 





a Thorburn ct al. 2003; L. Hobbs, priv. comm.; B. Rachford, priv. comm. 

b Savage ct al. 1977; Rachford ct al. 2002; Sheffer et al. 2003; Pan ct al. 2005; Cartlcdgc ct al. 2004 
c The three values are T 02 , T 4, and T k : sec Tabic 13. 

d Values in square braces assume that Nq/Nx is the same for 13 CO as for 12 CO. 



